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la  symmetric  ipecimeae  the  crack  advances  iato  the  relatively  undamaged  region  between  two 
plastic  shear  xoaes.  A  crack  aear  a  weld  or  shoulder,  loaded  iato  the  plaatic  range,  may  have  oaly  a 
siagie  shear  bead,  aioag  which  the  crack  grows  iato  p restrained  aad  damaged  material  with  lees 
ductility  thaa  the  usual  symmetrical  configurations.  A  crack  ductility  caa  be  defined  as  the 
minimum  displacement  per  unit  crack  growth.  A  low  crack  ductility  requires  higher  stiffness  of  the 
surrounding  structure  for  fracture-stable  design.  Tests  of  six  alloys  showed  that,  for  the  low-  j. 

hardening  alloys,  the  crack  ductility  ia  the  asymmetric  case  is  less  thaa  a  third  that  of  the 
symmetric.  In  the  higher  hardening  alloys  the  crack  ductility  ia  the  asymmetric  case  is  smaller  by  a 
factor  of  1.2  at  most.  A  noteworthy  result  is  the  presence  of  a  Mod*  1  opeaiag  component  even  with 
asymmetry,  as  is  shown  by  the  far  field  displacement  vector  being  more  thaa  45*  from  the  transverse 
direction.  The  crack  direction  is  less  than  45*  indicating  the  effect  of  triaxiaiity  oa  cracking. 

A  macro-mechanical  model  for  crack  advance  by  combined  fracture  aad  sliding  off  along  two 
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slip  planes  gives  the  independent  physical  parameters  (crackiag  aad  two  shear  directions,  relative 

amounts  of  crackiag  aad  shearing)  ia  terms  of  the  observable  quantities  of  the  macroscopic  fracture 

(flank  angles,  flank  lengths,  back  sagle).  This  two  slip  plane  model  admits  a  Mode  I  opening  - 

component  and  describes,  based  oa  an  idealization  of  underlying  physical  mechanisms,  the 

development  of  deformation  ia  ductile  crack  growth  for  both  the  symmetric  and  asymmetric 

specimens.  A  finite  element  study  of  the  asymmetric  specimens  gave  a  crack  direction  within  two 

degrees  aad  a  far  field  displacement  vector  at  initiation  within  three  degrees  of  that  experimentally 

found.  Stress  aad  strain  fields  indicate  the  presence  of  a  Mode  I  component.  Early  growth,  studied  , 

by  successive  removal  of  the  most  damaged  element,  resulted  ia  crack  growth  rate  for  the  lower 

hardening  cane  about  twice  that  of  the  higher  hardening  one.  Finally,  a  logarithmic  tensile 

singularity  in  the  mean  normal  stress  was  found  for  rigid-plastic  flow  past  a  growing  crack  of  finite 

angle  with  rigid  flanks  under  combined  shear  aad  tension.  The  tensile  singularity  predicts  yielding  of 

the  crack  flanks. 
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CHAPTER  ONE 


INTRODUCTION 

In  symmetric  singly  grooved  unconstrained  tensile  specimens  the  crack 
advances  into  the  relatively  undamaged  region  between  two  symmetric  shear  zones. 
An  asymmetry,  introduced  through  a  weld  fillet  or  a  harder,  heat-affected  zone  or  a 
shoulder  on  one  side  of  the  crack  suppresses  one  of  the  two  slip  lines  that  would 
appear  in  a  symmetrical  specimen.  This  is  likely  to  give  asymmetric  cracking  along 
the  remaining  active  slip  line,  with  less  ductility  because  the  crack  is  advancing  into 
prestrained  and  pre-damaged  material.  A  reduced  ductility  requires  higher  stiffness 
of  the  surroundings  for  fracture-stable  design. 

Near  the  tip  of  the  growing  crack,  strain  hardening  will  cause  the  deformation 
field  to  fan  out.  For  power  law  creep  or  deformation  theory  plasticity  and  a 
stationary  crack,  the  asymptotic  stress  and  strain  distribution  may  be  found  from 
the  extended  by  Shih  [1]  HRR  [2,3]  solutions  for  the  general  mixed  mode  I  and  II 
case.  Notice,  however,  that  such  a  superposition  of  stationary  singularities  does  not 
take  into  account  the  hardening  of  the  mate  il  left  behind  the  growing  crack. 
Indeed,  the  stress  and  strain  fields  near  the  tips  of  growing  cracks  in  ductile 
materials  are  known  to  differ  from  the  stress  and  strain  state  around  stationary 
cracks  in  the  same  materials  as  is  shown  from  asymptotic  solutions  [4,5,6] 
supplemented  through  finite  element  calculations  [7,8]. 

McClintock  and  Slocum  [9]  developed  a  formulation  for  the  accumulation  of 
damage  directly  ahead  of  an  asymmetric  crack,  based  on  strain  increments  adapted 
from  Shih's  [1]  analysis.  The  crack  was  assumed  to  follow  the  center  of  the  45° 
shear  band.  It  was  found  that  the  crack  growth  per  unit  displacement  increases 


approximately  as  the  logarithm  of  the  total  crack  advance  per  inclusion  spacing  p 
and  varies  inversely  as  the  critical  fracture  strain  7  .  Little  effect  of  strain 
hardening  on  the  growth  rate  was  found. 

The  objective  of  the  current  study  is  to  investigate  through  experimental, 
analytical  and  numerical  work,  the  ductility  of  asymmetric,  fully  plastic, 
unconstrained  configurations.  First,  the  approximate  pure  Mode  II  incremental 
solution  [9]  is  extended  to  admit  a  crack  growing  at  an  angle  to  the  shear  band. 

This  deviation  from  the  shear  band  is  expected  from  the  higher  triaxiality.  Far  field 
displacement  is  assumed  again  to  be  parallel  to  the  shear  band.  Next,  tests  results 
on  symmetric  and  asymmetric  specimens  of  six  alloys  are  presented.  A  method  for 
quantifying  and  representing  the  ductility  is  suggested.  A  macromechanical  model  of 
crack  growth  by  combined  fracture  on  one  plane  and  sliding  off  along  two  others 
describes,  for  this  idealization  of  the  physical  mechanisms,  the  ductile  crack  growth 
for  both  the  asymmetric  and  symmetric  specimens.  To  account  for  the  effect  of  the 
finite  width  of  the  shear  band  and  study  the  stress  and  strain  fields  at  initiation,  a 
finite  element  investigation  is  undertaken.  Early  growth  is  also  studied  by  successive 
element  removal.  Finally,  a  stream  function  technique  is  used  to  investigate  whether 
rigid-plastic  strain  hardening  flow  past  a  growing  crack  of  finite  angle  with  rigid 
flanks  can  be  sustained.  The  last  chapter  contains  an  o  rerview  of  the  results  and 
summarizes  the  conclusions.  It  also  contains  recommendations  for  further  research. 
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ABSTRACT 


Welds,  shoulders,  or  other  asymmetries  may  eliminate  one  of  the  shear  bands 
of  symmetrically  cracked  parts  and  thus  give  crack  propagation  through  pre¬ 
damaged  material,  instead  of  through  the  relatively  unstrained  region  between  the 
two  plastic  shear  zones  of  the  symmetric  case.  Previous  work  is  extended  to  include 
sites  at  several  angles  ahead  of  the  crack.  Far  field  displacement  is  assumed  to  take 
place  parallel  to  the  sheai  band.  Strain  increments  are  approximated  from  the 
mixed  mode,  power-law  elastic  solution  for  a  stationary  crack  and  used  with  a 


g 


a 


r* 
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fracture  criterion  for  hole  growth  in  shear  bands  to  predict  the  direction  and  rate  of 
crack  growth.  The  crack  is  assumed  to  advance  to  advance  in  the  direction  that 
requires  the  minimum  far  field  displacement  to  reach  critical  damage.  For  a  shear 
band  at  45°  the  crack  progresses  at  21°-30°  from  the  transverse  (depending  on  strain 
hardening),  indicating  the  effect  of  higher  triaxiality.  The  crack  growth  rate  is  about 
6-15%  higher  than  if  the  directional  effects  are  neglected.  Lower  strain-hardening 
results  in  a  5%  higher  rate  of  crack  advance  per  unit  displacement,  a  higher  fracture 
strain,  and  the  final  crack  orientation  being  closer  to  the  45°  shear  band. 

INTRODUCTION 

Most  fracture  tests  use  symmetric  specimens.  The  crack  advances  into 
relatively  undamaged  material  between  two  shear  bands.  This  will  not  happen  if  one 
of  the  bands  is  eliminated  due  to  a  weld  bead,  or  a  harder  heat-affected  zone,  for 
example  (Fig.  1).  A  fatigue  crack  or  other  defect  near  such  an  asymmetry  will  tend 
to  advance  along  the  remaining  shear  band  through  highly  strained  material.  Lower 
ductility  is  thus  expected.  An  example  of  lowered  ductility  in  asymmetric  flow  is  the 
formation  of  a  shear  lip  at  the  end  of  an  ordinary  cup  and  cone  fracture  in  a  tensile 
test. 


McClintock  and  Slocum  (ij  developed  an  approximate  formulation  for  the 
accumulation  of  damage  ahead  of  the  crack  in  a  power-law  strain  hardening 
material,  by  using  the  strain  and  displacement  fields  derived  by  Shih  [2]  for  a 
stationary  mixed-mode  crack  and  the  McClintock,  Kaplan,  and  Berg  [3]  criterion  for 
fracture  by  hole  growth.  It  was  assumed  that  the  crack  advances  directly  along  the 
shear  band.  Preliminary  experiments,  however,  have  indicated  that  the  crack 
actually  advances  at  an  angle  from  the  shear  band  which  reflects  the  effect  of  the 


higher  triaxiality  on  one  side. 


In  the  following  we  modify  the  Pure  Mode  II  [1]  solution  by  considering  several 
sites  around  the  crack  tip.  The  far  field  displacement  is  again  assumed  to  be  parallel 
to  the  shear  band  and  the  strain  increments  to  follow  Shih’s  [2)  mixed  mode 
stationary  crack  fields.  The  accumulated  damage  from  initiation  and  prior  growth  is 
calculated  and  the  necessary  far-field  displacement  for  critical  damage  is  found  for 
each  site.  The  crack  is  assumed  to  advance  in  the  direction  requiring  the  least  far- 
field  displacement. 
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*1  r  Ii/n(n*MP) 


D^“+1)  fTjj(0,Mp,n) , 


£ij  = 


<ri  r  I1/n(n,MP) 


l/(n+l) 


Cy^.MP.n) , 


ui 


<rl  r  I1/n(n,MP) 


i/(»+l)  . 


u^.MP.n) 


The  J  integral  is  defined  (Xj  axis  along  crack)  as 


f  5ui 

J  =  /  Wdx„  -  T.— ids  . 

J  2  id\l 


(3) 


where  Tj  is  the  traction  vector,  Uj  is  the  displacement  vector,  and  W  is  the  work  per 
unit  volume.  For  a  single  shear  band,  J  can  be  evaluated  in  terms  of  the  shear 
strength  k,  the  far  field  relative  displacement  U  and  the  angle  <(>  between  crack  and 
shear  band  (Fig.  2)  by 
kU 

J  = - .  (4) 

COS  <f> 

The  dimensionless  functions  <7^(0, Mp,n),  £jj(0,Mp,n)  and  I^Jn.MP)  have  been 
numerically  determined  by  Shih  (2]  for  n=l/3  and  n=l/13.  The  dimensionless 
functions  u^.MP.n)  are  derived  in  the  Appendix  from  the  strain  functions  and  are 
shown  in  Fig.  3  for  n=l/13  and  Mp=0.5. 


Assume  that  currently  the  shear  band  forms  an  angle  <j>  with  a  recent  average 
crack  direction  to  be  defined  below.  With  the  crack  not  advancing  along  the  shear 
band,  there  is  no  longer  pure  Mode  II.  When  the  relative  far-field  displacement  is 
assumed  to  be  parallel  to  a  single  narrow  shear  band,  as  is  valid  in  the  non-hardening 
limit,  and  that  direction  is  assumed  to  be  the  same  as  the  local  relative  displacement 
across  the  flanks,  the  mixity  parameter  can  be  determined  from  the  angular 
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functions  (')  of  the  displacement  field  relative  to  0=-ir  since,  from  Fig.  2, 
ue  u,(jr,Mp,n) 

tan0  =  —  = - .  (5 

uf  ur(5T,Mp,n) 

Fig.  4  shows  the  resulting  variation  of  the  mixity  parameter  with  the  angle  <p  for 
n=l/13.  Thus,  the  angle  4>  determines  the  angular  stress  and  strain  functions  and 
hence  the  local  stress  and  strain  for  a  given  J  through  (3).  The  angular  functions 
turn  out  to  affect  the  fracture  criterion  through  the  triaxiality  and  the  shear  strain, 
as  follows.  The  mean  normal  stress  for  plane  incompressible  flow  is 

*rr+<TM 


The  triaxiality  <r/r  used  in  the  fracture  criterion  is  given  in  terms  of  the 
dimensionless  principal  shear  stress  f,  defined  by 

-  f-  2  ^  /  V*^2ll/2 
r=  <V  +  (-T— ) 

it 


The  angular  variation  of  this  triaxiality,  afr,  is  shown  in  Fig.  5  for  n=l/13.  Note 
that  the  triaxiality  is  highest  for  negative  values  of  6  for  all  cases  except  pure  Mode 
I.  This  is  the  primary  reason  for  exploring  the  directional  effects. 

Similarly,  the  dimensionless  principal  shear  strain  7  can  be  expressed  for  the 
plane  incompressible  case  as:  , 

i  =  4,2+(^)T/2  =  2^,, 2  +  V  m 

it 

Introducing  (4)  into  the  first  of  (3)  and  solving  for  the  displacement  gives  the  far- 
field  displacement  U  in  terms  of  the  principal  shear  strain  7  at  any  point  in  the  near 
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field: 

—  —  I.  /„(n,Mp)  rcos<£  (7/;r)ll+1  . 

A. 


(10) 


The  critical  displacement  for  crack  initiation  occurs  when  the  fracture  strain  is 
reached  at  the  point  (p,0)  where  p  is  a  fracture  pro*.  >  •'•»<»  size  (e.g.  the  mean 

inclusion  spacing).  The  fracture  strain  is  found  by  using  the  fracture  >°rion  of 
McClintock,  Kaplan  and  Berg  [3]  by  which  it  is  postulated  that  fracture  due  to 
micro-void  coalescence  occurs  when  the  "damage”,  t),  reaches  unity.  The  damage  is 
expressed  in  terms  of  a  hole  growth  ratio  Ft,  the  principal  shear  strain  7  and  the 
triaxiality  <r/r. 


= - [lnv/l+V*  + 

InF. 


(l-n)<r 
-sinh( - )]  . 


t  2(l-n)  t 

The  above  damage  equation  is  associated  with  growth  of  cylindrical  holes. 
Alternatively  one  might  use,  for  example,  the  eqs  for  growth  of  spherical  holes  in 
nonhardening  material  (Rice  and  Tracey,  (5)). 


(11) 


A  Newton-Raphson  technique  is  used  to  solve  (10)  and  (11)  at  r =p  for  the 
critical  far-field  displacements  for  initiation  in  a  number  of  directions.  The  actual 
initial  direction  is  that  which  minimizes  the  required  displacement.  Once  the 
initiation  displacement  Uj  for  the  critical  strain  at  the  point  (p,  0.)  is  known,  the 
strain  at  all  other  points  can  be  found  by  re-arranging  (10): 


kU 

7=(_1 _ ]l/(o+l)^. 

nr^n.MPjcos^ 


(12) 


t 


» 


2.  Crack  growth.  After  initial  growth  by  Ac,  further  growth  requires  reaching  the 
critical  damage  at  some  new  site  p  from  the  current  crack  tip.  The  damage  at  each 
site  is  that  from  crack  initiation  plus  those  for  any  following  crack  growth 
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increments.  The  strain  to  bring  the  damage  to  unity  is  found  by  differerentiating 
(11)  to  find  the  damage  increments  in  terms  of  the  strain  increment  and  the  strain 
itself: 


1.7  1 

d  tj  = -  - -  + 


(l-n)<r 

-sinh - ]  d^ 


(13) 


lnFt  I+72  2(l-n)  r 

In  the  absence  of  an  incremental  strain-hardening  solution  for  a  growing  crack,  we 
follow  McClintock  and  Slocum  [1]  and  approximate  the  strain  increment  in  terms  of 
the  far-field  displacement  increment  by  differentiating  and  rearranging  (10)  (this  is 
strictly  valid  only  for  non-linear  elasticity): 
k7 


= 


(7/7)“dU 


(14) 


(n+l)<r1I(n,Mp)rcos^ 

The  damage  at  any  point  in  front  of  the  growing  crack  is  given  by  the  sum  of  the 
damage  due  to  crack  initiation,  as  found  from  (12)  and  (11),  and  all  of  the  damage 
increments  from  prior  crack  growth,  as  found  from  (13)  and  (14),  with  r  taken  to  be 
the  distance  from  the  prior  crack  tip  to  -he  roint  in  question,  and  <f>  the  prior  angle 
between  crack  and  shear  band. 


The  necessary  increment  in  damage  for  fracture  is  6ij=l-ri.  The  corresponding 
strain  increment  can  be  found  from  (13): 


&1  =  (l-»j)lnFt  / 


1+7  2(l-n) 


-sinh 


(l-n)<r 


(15) 


The  necessary  increment  in  far-field  displacement  to  cause  this  strain  increment  can 
then  be  found  from  (14): 

(^I^tn.MP^cos^ 


SU  =  (n+1) 


7k 


(7/7)"^  . 


(16) 


The  crack  will  advance  in  the  direction  <j>  which  requires  the  least  far  field 
displacement  to  reach  critical  damage,  not  necessarily  toward  the  most  severely 


i.e.  the  angle  between  crack  and  shear  band  is  smaller  and  the  triaxiality  is  smaller. 

(5)  For  the  low  strain  hardening  n=l/13,  increasing  gave  smaller  initiation 
displacements  and  strains  and  larger  initial  crack  growth  rates  (dc/du). 

(6)  Strains  and  triaxialities  during  growth  are  relatively  insensitive  to  the  initial 
angle  between  crack  and  shear  band. 

(7)  For  both  strain  hardening  exponents  and  all  the  angles  <^0,  the  final  average  angle 
between  crack  and  shear  band  0avg  after  growth  by  c/p=100  was  between  23°  and 
32°  from  the  shear  band.  The  Mode  I  mixity  corresponding  to  the  final  crack 
orientation  was  also  within  a  correspondingly  narrow  range  for  each  of  the  strain 
hardening  exponents. 

(8)  These  rigid-plastic  results  do  not  predict  instability  (infinite  crack  advance  per 
unit  far  field  displacement).  Instability  could,  however,  arise  from  the  compliance  of 
the  surrounding  structure. 

SUMMARY  OF  CONCLUSIONS 

A  part  containing  a  crack  near  a  weld  or  a  shoulder,  loaded  into  the  plastic 
range,  can  give  an  asymmetric  shear  band  extending  from  the  crack  tip.  The 
resulting  crack  propagation  into  previously  damaged  material  gives  less  ductility 
than  the  typical  symmetric  case.  A  previous  incremental  solution  for  crack  growth 
using  Shih’s  asymptotic  fields  for  a  stationary  crack  in  nonlinear  elastic  material  is 
extended  to  account  for  the  effect  of  triaxial  stress  in  advancing  a  crack  at  an  angle 
to  the  shear  band.  Far  field  displacement  is  assumed  to  take  place  along  the  shear 
band.  Cracking  is  assumed  to  occur  at  the  site  around  the  crack  tip  that  needs  the 
least  far  field  displacement  for  critical  damage.  For  a  45°  shear  band,  it  is  found 
that  the  crack  does  not  advance  along  the  shear  band  but  at  an  angle  of  about  21° 
from  the  transverse  under  a  higher  triaxiality.  The  crack  growth  rate  is  higher  by 


about  6-15%  (the  larger  increase  with  less  hardening)  than  if  the  directional  effects 
are  neglected  and  the  crack  is  assumed  to  progress  along  the  shear  band.  A  higher 
strain  hardening  decreases  slightly  (about  5%)  the  crack  growth  rate  and  the  final 
angle  from  the  transverse  of  the  growing  crack.  Strains  and  triaxiaiity  during 
growth  are  not  sensitive  to  the  initial  angle  between  the  crack  and  the  shear  band. 
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TABLE  1 

Comparison  of  numerical  and  approximate  pure  Mode  II  solutions 

Numerical  Approx.  Mode  H 
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Appendix  -  Displacement  Functions 
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The  displacement  functions  Uj  are  determined  from  the  strain  functions  for  the 
plane  strain  considered  here.  The  radial  displacement  uf  may  be  found  from  the 
radial  strain 
<9u 


€rr  = 


dr 


(17) 


so 


Ur  =  £ errdr  +  Ur(0,i9)  . 


(18) 


For  zero  rigid-body  translation  at  r=0,  ur(O,0)=O.  Eliminating  e  with  (3)  and 
integrating  gives 

u  J  ,,  n+1 

~=[  -  L/(°+1)( - ] *„  ■  (19) 

r  (T1rIl/n(n,MP)  “  * 

Using  the  displacement  equation  (3)  gives  the  radial  displacement  function  uf  relative 
to  that  at  0=-ir. 


n+1 

ur(0,Mp,n)  =  — — (€rr(0,Mp,n)  -  in(-v, Mp,n))  . 
n 

The  tangential  displacement  function  u^  is  determined  from 
1  1  dwa 


(20) 


f-aa  =  -U.  +  - 


r  dO 


(21) 


as 


ne  =  f  K<fur 

J-n 


)d  9  +  f(r)  . 


(22) 


Noting  that  ( ^  —  -crr  for  plane  strain  incompressibility  and  using  (3)  with  (20)  gives 
the  tangential  displacement: 


r  ajrlj^fn.MP)  n  in 


etTd6  +  f(r)/r  . 


(23) 
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With  respect  to  the  displacement  at  0=-n,  f(r)=0.  By  using  (3)  we  can  thus  find  the 

dimensionless  tangential  displacement  function  relative  to  the  displacement  at  9=-ir, 

in  terms  of  the  dimensionless  strain  function  e  : 

rr 

2n+l  / 

ufl(0,Mp,n)  = - /  £rr(0,Mp,n)d0  .  (24) 

n  JLn 

As  an  example,  the  displacement  functions  for  n=l/13,  Mp=0.50,  determined 
numerically  by  (2 1),  (20)  are  given  in  Fig.  3. 


mgular  variation  of  the  radial  and  tangential  near  tip 
dimensionless  displacement  functions  for  plane  strain 
with  n-1/13  and  Up=0  5 
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EXPERIMENTAL  STUDY 

TABLE  OF  SYMBOLS 

minimum  extension  rate  (eq.  10) 
crack  ductility  (eq.  3) 
modulus  of  elasticity 
load  factor  (eq.  12) 

J-integral 

shear  yield 

initial  ligament 

strain  hardening  exponent 

nominal  load 

tensile  strength 

tearing  modulus 

eq.  8 

eq.  9 

initiation  displacement 

growth  displacement 

total  displacement  vector 

growth  displacement  vector 

specimen  width 

yield  strength 

mean  inclusion  spacing 

flow  stress  at  unit  strain 

displacement  vector  angle  from  transverse 

fracture  strain 

crack  opening  angle 

crack  direction  from  transverse. 


ABSTRACT 

Most  fracture  tests  use  symmetric  specimens,  with  the  crack  advancing  into 
the  relatively  undamaged  region  between  two  plastic  shear  zones.  However,  a  crack 
near  a  weld  or  shoulder,  loaded  into  the  plastic  range,  may  have  only  a  single  shear 
band,  along  which  the  crack  grows  into  prestrained  and  damaged  material  with  less 
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ductility  than  the  usual  symmetrical  configurations.  An  experimental  study  on  six 
alloys  shows  that  while  the  crack  initiation  displacements  are  similar,  the  growth 
displacement  is  much  less  for  the  asymmetrical  specimens,  especially  with  less 
hardening.  Indeed,  for  the  low-hardening  alloys  (n^O.l)  the  crack  growth  ductility, 
defined  as  the  minimum  displacement  per  unit  ligament  reduction,  is  less  in  the 
asymmetric  case  than  the  symmetric  by  a  factor  of  three,  In  the  higher  hardening 
alloys  the  crack  growth  ductility  is  less  in  the  asymmetric  case  by  a  factor  of  1.2  at 
most.  Triaxiality  on  one  side  of  the  asymmetric  shear  crack  diverts  it  from  45°  to 
38°-41°  from  the  transverse  direction,  the  larger  angles  with  smaller  strain 
hardening.  In  addition,  the  far  field  displacement  vector  is  more  axial  than  the  45° 
line,  at  51°  to  63°  from  transverse,  suggesting  a  Mode  I  component  even  with 
asymmetry. 

INTRODUCTION 

For  fracture-stable  structures  it  is  important  not  only  that  fully  plastic 
conditions  be  attained  before  fracture,  but  also  that  the  load  does  not  fall  off  too 
rapidly  during  crack  growth.  Flow  fields  such  as  Fig.  1,  in  which  the  far-field 
deformation  consists  of  a  single  shear  band,  may  arise  in  practice  due  to  the 
constraint  of  weld  material.  These  specimens  may  exhibit  less  ductility  than  the 
symmetric  ones,  because  the  crack  is  advancing  into  pre-strained  and  damaged 
material,  rather  than  into  the  new  material  encountered  by  a  crack  advancing 
between  two  symmetrical  shear  bands.  Being  able  to  predict  such  increased  crack 
growth  can  have  useful  applications  in  the  design,  inspection  and  maintenance  of 
pressure  vessels  and  ships. 

Near  the  tip  of  the  growing  crack,  strain  hardening  will  cause  the  deformation 


field  to  fan  out.  For  power  law  creep  or  deformation  theory  plasticity,  the  stress  and 
strain  in  the  neighborhood  of  a  stationary  crack  may  be  found  from  Shih’s  [1]  mixed 
mode  solutions.  More  realistically,  a  corresponding,  fully-plastic,  incremental 
plasticity  solution  should  be  obtained  for  a  growing  crack,  taking  into  account  the 
hardening  of  the  material  left  behind  the  growing  crack.  McClintock  and  Slocum  {2] 
developed  a  formulation  for  the  accumulation  of  damage  directly  ahead  of  an 
asymmetric  crack,  based  on  strain  increments  adapted  from  Shih’s  [l]  analysis  for 
stationary  cracks  in  a  power  law  material.  The  crack  was  assumed  to  follow  the 
center  of  the  45°  shear  band.  It  was  found  that  the  crack  growth  per  unit 
displacement  increases  approximately  as  the  logarithm  of  the  total  crack  advance  per 
inclusion  spacing  p  and  varies  inversely  as  the  critical  fracture  strain  7  .  To  correct 
for  triaxial  effects,  several  sites  around  the  current  crack  tip  were  considered  in 
chapter  2.  The  damage  at  each  site  due  to  crack  initiation  and  prior  growth  was 
determined  and  then  the  necessary  increment  in  far  field  displacement  was  found  for 
each  site.  The  crack  was  assumed  to  advance  in  the  direction  requiring  the  least 
displacement.  This  numerical  investigation  resulted  in  growth  directions  not  along 
the  45°  shear  band  but  at  a  smaller  angle  from  the  transverse  depending  on  the 
hardening  and  the  initial  crack-shear  band  angle,  and  lower  ductility  by  6-15°^ 

(larger  decrease  with  less  hardening)  than  with  growth  along  the  shear  band. 

A  test  with  pure  shear  (Mode  II)  loading  was  carried  out  by  Chant  et  al.  [3]  of 
high  hardening  carbon  manganese  steel  (B.S.  I501-151-430A,  Y.S.=329  MN/mr,  T.S. 
=  490  MN/m2).  Small  specimens  were  subjected  to  both  Mode  II  and  Mode  I 
testing  but  the  ductility,  measured  by  dJ/da,  was  practically  the  same  although  the 
microscopic  features  for  the  pure  shear  specimens  are  different  than  those  observed 
in  the  Mode  I  specimens.  The  objective  of  (he  current  study  is  to  present 
experimental  evidence  on  the  ductility  of  asymmetric  crack  configurations. 


EXPERIMENTAL  PROCEDURES 


Material.  Tests  were  performed  on  six  alloys  with  the  mechanical  properties 
listed  in  Table  1.  True  stress-true  plastic  strain  curves  (Figs.  2a, b)  for  these  alloys 
were  obtained  using  standard  6.35  mm.  dia.  specimens  with  25.4  mm.  gage  length.  It 
is  convenient  to  represent  a  curve  of  equivalent  stress  <t  vs.  equivalent  plastic  strain 
t*  by: 

°  =  *1(fo+rP)D’  (!) 

where  <Xj,  eQ  and  n  are  three  constants  that  were  determined  from  the  flow  strengths 
at  yield  point,  e=0.125  and  £=0.250,  and  are  given  in  Table  2.  The  lower  hardening 
alloys  are  the  1018  cold  Finished  steel,  HY-80  and  HY-100  steels  (n~0.10)  and  the 
higher  hardening  alloys  are  the  1018  normalized  and  A36  hot  rolled  steels  (n~0.24). 
The  5086-Hlll  aluminum  is  between  these  two  groups. 

Test  Method.  From  12.7  mm.  dia.  round  bars  of  each  alloy,  seven  specimens 
were  First  machined  as  shown  in  Fig.  3a,  with  side  grooves  to  ensure  a  straight 
fatigue  pre-crack  approximately  1.3  mm  deep.  For  the  four  asymmetric  specimens 
(Fig.  3b),  further  side  grooves  were  machined  at  40°  from  the  transverse  direction. 
This  corresponded  to  the  crack  direction  found  in  preliminary  tests  and  served  to 
reduce  3-dimensional  effects.  For  the  three  symmetric  specimens,  since  the  crack 
grows  by  alternating  shear  at  ±45°,  orthogonal  triangles  were  machined,  as  shown  in 
Fig.  3c. 

Stability  of  the  tests  turned  out  to  be  an  important  consideration  due  to  the 
high  crack  growth  rate  expected  in  the  asymmetric  case.  Thus  short  specimens,  stiff 
adapters,  and  locknuts  were  used.  The  tensile  tests  were  performed  on  an  MTS  50 
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metric  ton  testing  machine  with  resulting  compliances  of  2.3xl0'6,  4.SxlO‘6,  1.08x1  O'6 
mm/N  for  the  steel  specimens,  the  adapters,  and  the  machine  respectively.  The 
axial  and  transverse  displacements  across  the  notch  were  plotted  continuously. 

A  typical  plot  of  load  vs.  axial  displacement  is  shown  in  Fig.  4;  the 
breakthrough  point  is  when  the  fracture  first  breaks  through  the  back  surface,  with 
some  shear  lips  remaining  on  the  sides.  The  displacement  during  crack  initiation  and 
growth,  u;  and  ug,  are  found  from  the  drawings  of  the  crack  path.  The  topographies 
of  the  fracture  surfaces  were  thus  subsequently  plotted  using  a  metallurgical 
microscope  with  a  travelling  stage.  The  horizontal  and  vertical  coordinates  of  the 
travelling  stage  are  recorded  with  two  linear  potentiometers;  several  points  are 
obtained  to  give  an  impression  of  the  surface  profile  of  the  broken  specimens.  A 
typical  microscope  plot,  as  in  Fig.  5,  consists  ol  the  60°  notch,  the  fatigue  crack 
(with  some  amount  of  deformation,  v^-v^),  an  initiation  zone  which  shows  some 
blunting,  and  a  growth  zone.  The  initiation  displacement  is  v^-v^.  These  quantities 
were  also  checked  against  the  data  from  the  load-extension  curves.  In  addition, 
fracture  profiles  were  used  to  obtain  the  angular  quantities  such  as  the  crack  opening 
angle,  u ;,  the  lower  and  upper  flank  angles,  9 j  and  9 u,  and  the  orientation  of  the  total 
displacement  vector,  <t>,  in  the  asymmetric  case. 

RESULTS 

Initiation  Displacement.  Stable  tests  were  obtained  except  for  the  lower 
hardening  alloys,  which  were  unstable  for  less  than  209o  of  the  falling  part  of  the 
load-displacement  curve.  The  results  of  the  tests  are  summarized  in  Table  3.  An 
idealized  initiation  displacement,  u^/g,  can  defined  as  the  extension  between  the 
initial  elastic  loading  and  the  steepest  unloading  parts  of  the  load-displacement  curve 


at  maximum  load,  normalized  with  the  initial  ligament  lQ  (Fig.  4).  This  quantity, 
given  in  the  first  row  of  the  table,  is  a  convenient  measure  of  initiation  and  can  be 
compared  with  the  the  initiation  displacement  Uj,  measured  from  the  fracture  surface 
profiles  after  complete  separation.  The  normalized  crack  initiation  displacement 
Uj//0,  does  not  appreciably  differ  between  the  asymmetric  and  symmetric 
configurations.  It  is,  however,  dependent  on  the  strain  hardening,  being  for  the 
higher  hardening  alloys  two  to  four  times  that  of  the  lower  hardening  ones. 


Ductility.  For  a  measure  of  the  crack  growth  resistance,  the  crack  ductility, 

D  ,  is  defined  as  the  minimum  displacement,  du  ,  per  unit  projected  ligament 
reduction,  dl.  Thinning  of  the  ligament  from  the  far  side  in  fully  plastic  flow  makes 
the  reduction  in  ligament  rather  than  crack  advance  more  appropriate  for  describing 
load  drop.  The  displacement  duc  is  associated  with  the  crack  opening  stretch  and 
consists  of  the  gauge  displacement  du  and  the  elastic  unloading  duun,  (Fig.  4): 

duc  =  du  +  duun,  .  (2) 


The  ligament  reduction  can  be  approximated  from  the  relative  load  drop  and  thus 
we  can  define: 


D  = 

g  VdP/p _ J  min  \di  Jmi 


1vdP/Pmax-  min  "*•  '  min 

From  Fig.  6,  this  is  also  related  to  the  crack  opening  angle  by 


(3) 


Dg  ~  (COA)/cos20c  ,  (4) 

where  0.  is  the  crack  orientation.  In  addition,  the  above  defined  quantity  is  the 
normalized  compliance  requirement  for  fracture-stable  design: 


Compliance  of  surrounding  <  Dg/Q/Pmax  . 


(5) 
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The  crack  ductility,  given  in  the  second  row,  is  smaller  for  the  asymmetric  case  by  a 
factor  ranging  from  3.4  for  the  lower  hardening  HY100  to  1.1  for  the  higher 
hardening  A36  hot  rolled  steel.  Notice  that  in  the  lower  hardening  alloys  1018  CF, 
HY-80,  HY-100  steel  the  factor  by  which  the  ductility  is  reduced  is  larger  than  three 
which  shows  also  that  these  alloys  have  much  larger  stiffness  requirements  for 
stability.  A  comparison  of  Dg  among  the  alloys  reveals  that,  in  the  asymmetric  case, 
the  crack  growth  rate  in  the  lower  hardening  alloys  is  about  2  times  larger  than  in 
the  higher  hardening  alloys.  In  the  symmetric  case,  on  the  contrary,  the  crack 
growth  rate  is  practically  insensitive  to  strain  hardening. 


P  ^ 


r 


The  third  row  is  a  parameter  analogous  to  the  "tearing  modulus”  T  of  Paris  et 
al.  [4]  defined  in  terms  of  the  yield  or  tensile  strength  er0,  the  modulus  E  and  the  J- 
integral  by: 


E  dJ 
<tq2  dc 


(6) 


To  approximate  the  J-integral,  consider  the  simple  case  of  the  far-field  displacement 
taking  place  along  a  single  shear  band  [2]  and  express  it  in  terms  of  the  shear 
strength,  k,  and  the  displacement  along  the  band  \x\fil, 


J  =  kux/i", 


(7) 


and  thus  define  a  parameter,  T  ,  analogous  to  the  tearing  modulus  T,  which  allows 
comparing  the  ductility  of  alloys  of  different  strength.  In  terms  of  the  tensile 
strength  T.S.  ~  kx/^,  by: 


asym 


=WD‘- 


(8) 


In  the  symmetric  case  the  expression  for  the  J-integral,  J=2ku  [5]  leads  to  an 
analogous  to  (8)  expression. 


T  m  is  about  3  times  larger  for  the  higher  hardening  alloys,  as  is  shown  in  Table 
3. 

The  load-displacement  curve  of  Fig.  4  can  be  described  in  terms  of  the  initial 
elastic  compliance,  the  idealized  initiation  displacement  Uj1,  and  the  minimum  gauge 
displacement  per  unit  crack  advance  (steepest  slope  of  the  falling  part),  D  t,  given 
by 


This  definition  includes  the  effect  of  the  compliance  in  the  shoulders  and  is  thus 
smaller  than  Dg.  Results  for  25  mm  gauge  length  are  given  in  the  fourth  row  of 
Table  3. 

Growth  Displacement.  The  growth  displacement  until  the  fracture  breaks 
through  the  back  of  the  specimen,  ug,  can  be  found  from  the  fracture  surface 
profiles.  The  normalized  displacements  during  crack  growth,  ug//Q  is  more  than  3 
times  larger  in  the  symmetric  than  the  asymmetric  specimens  for  the  lower 
hardening  alloys  but  only  about  18%  larger  for  the  higher  hardening  A36  HR  steel. 
It  is  also  larger  by  about  a  factor  of  two  in  the  higher  hardening  relative  to  that  of 
the  low  hardening  alloys. 

The  far-field  displacement  vector  angle  from  the  transverse  in  the  asymmetric 
case,  defined  from  the  slope  of  the  transverse-axial  displacement  curves,  is  found  to 
be  greater  than  45°  and  larger  initially  in  the  lower  hardening  alloys.  As  the  crack 
grows  the  displacement  vector  becomes  less  axial  (Fig.  7a).  The  final  orientation  0  , 
measured  after  fracture  from  the  microscope  plots  (Fig.  5)  is  between  53°  and  63°, 
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larger  for  the  higher  hardening  ease.  The  fact  that  the  axial  component  of  the 
displacement  is  larger  than  the  transverse  one  suggests  a  Mode  I  mixity  of  the  local 
plastic  flow. 

Crack  direction.  In  the  symmetric  specimens  the  crack  runs  within  ten  degrees 
of  horizontal  except  for  shear  lips  near  the  ends  of  the  cracks.  In  the  lower 
hardening  alloys,  even  with  the  symmetric  geometry,  the  fracture  turned  often  into 
the  asymmetric  mode,  the  fracture  advancing  close  to  the  -15°  slip  plane  or,  in  some 
cases,  half  of  the  cross  section  following  the  one  and  the  other  half  following  the 
other  slip  plane.  In  the  asymmetric  specimens  the  crack  progresses  at  an  angle  of 
about  38-11°  from  the  transverse.  This  smaller  than  45°  angle  was  expected  from 
the  higher  triaxiality.  In  the  lower  strain  hardening  alloys  the  crack  grows  closer  to 
the  45°  band,  at  40°-41°  from  the  transverse  and  in  the  higher  hardening  alloys  at 
38°-39°  (Table  3).  Finally  more  blunting  occurred  with  the  higher  hardening  alloys 
and  in  the  symmetric  case. 

Load.  To  summarize  the  load  performance,  a  parameter  dealing  with  the 
maximum  load  will  be  defined.  The  nominal  load-carrying  ability  is  simply  the 
tensile  strength  multiplied  by  the  net  area  at  the  end  of  the  fatigue  crack  (and 
corrected  by  the  plane  strain  factor).  In  terms  of  the  initial  ligament  lQ,  the  width  w, 
and  the  tensile  strength  T.S., 

P-on,  =  <,w(T.S.X2/n/3)  '  (><> 

A  load  factor  F,  can  be  defined  in  terms  of  the  actual  maximum  load  P  as: 

L  max 

F,  =  P  /P  (12) 

L  max'  nom  '  ' 

Table  3  also  contains  the  load  factors.  They  are  in  general  larger  in  the  symmetric 


case,  where  the  overall  deformation  is  bigger,  and  in  the  higher  hardening  alloys. 


The  normalized  load-normalized  displacement  (and  transverse-axial 
displacement  for  the  asymmetric  case)  curves,  obtained  for  the  lower  hardening 
HY-100  steel  and  higher  hardening  A36  HR  steel  for  both  the  asymmetric  and 
symmetric  case  are  shown  in  Figs.  7a,  7b,  8a  and  8b.  For  the  HY-100  steel,  notice 
the  sharp  increase  in  the  slope  of  the  falling  part  of  the  load-displacement  curve  of 
the  asymmetric  case  relative  to  that  of  the  symmetric;  this  is  not  the  case  in  the  A36 
HR  steel.  The  microscope  plots  for  these  alloys  are  shown  in  Figs.  9a,  9b  10a  and 
10b.  In  the  HY-100  steel,  notice  the  large  reduction  in  the  crack  opening  angle  of 
the  asymmetric  case  relative  to  that  of  the  symmetric,  whereas  in  the  A36  steel  the 
difference  in  the  crack  opening  angle  between  the  two  geometries  is  not  appreciable. 

Size  effects.  To  investigate  size  effects,  tests  were  performed  in  38.1  mm.  dia. 
specimens  of  5086-Hlll  aluminum  and  the  results  were  compared  with  those  from 
the  12.7  mm.  specimens.  Table  4  summarizes  the  results  Comparing  with  the  data 
given  in  Table  3  for  the  smaller  5086-Hlll  specimens,  we  conclude  that  the  ductility 
and  the  normalized  growth  displacement  is  only  4%  smaller  for  the  larger  specimens 
and  the  load  factor  is  slightly  larger.  Notice  that  the  size  effects  that  were  predicted 
in  [2]  are  associated  with  a  transient  behavior  (increasing  crack  advance  per  unit  far 
field  displacement). 

Marking  the  crack  front.  In  the  large  38.1  mm  dia.  5086-Hlll  aluminum 
specimens,  the  crack  front  was  marked  by  imposing  unloading-loading  cycles  at 
selected  points  during  crack  advance.  The  spacing  of  these  fatigue  marks  was 
measured  with  a  stereo  microscope  at  about  50x.  The  corresponding  displacements 
were  then  obtained  from  the  load-displacement  curves.  In  this  manner,  points  on  the 
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c-u  curve  can  be  accurately  determined.  Figure  11a  and  11b  show  the  load- 
displacement  curves  and  in  Figure  12  the  corresponding  crack  growth-displacement 
data. 


Comparing  with  theoretical  formulations.  An  approximate  formulation  for  the 
accumulation  of  damage  ahead  of  an  asymmetric  crack,  based  on  strain  increments 
following  a  power  law  relationship  was  presented  by  McClintock  and  Slocum  (2). 

The  crack  was  assumed  to  progress  along  the  45°  shear  band  (pure  Mode  II)  with  the 
far  field  displacement  along  the  shear  band.  The  initiation  displacement  was 
expressed  by: 


2c 


(13) 


where  7c  is  the  fracture  strain,  k  is  the  shear  strength,  c  =  0.88  for  the  assumed 
pure  Mode  II  and  1^  =  0.72-0.83  for  n  =  0. 1-0.2.  The  initiation  displacement  is 
thus  of  the  order  of  the  inclusion  spacing  (0.010  mm),  much  smaller  than  the  one 
found  experimentally.  This  discrepancy  is  due  to  the  blunting  that  occurs  during 
crack  initiation.  For  a  quasi-steady  growth,  the  crack  advance  per  unit  displacement 
was  practically  insensitive  to  the  strain  hardening  exponent  n  and  was  found: 
d(u/u.)  n+1 

- -  = - .  (14) 

d(c/p)  ln[(c-Cj)/p  +  exp(n-fl)] 

The  above  formula,  for  a  mean  inclusion  spacing  p  =  .01  mm  and  growth  by  the 
ligament  length  of  lQ  =  2.54  mm,  gives  du/dc  «  0.200,  which  is  closer  to  the  test 
data  for  the  higher  hardening  alloys.  Equation  (14)  underestimates  the  crack  growth 
rate  in  the  lower  hardening  alloys  by  a  factor  of  two.  For  the  size  effects,  predicted 
in  (2],  use  the  mean  inclusion  spacing  of  about  10  microns  and  find  the  ratio  of  the 
crack  growth  rate  for  the  large  38.1  mm  specimens  (initial  ligamemt  /Q=7.62  mm)  to 
that  of  the  small  12.7  mm  (initial  ligament  lQ=  2.54  mm)  ones  as 
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(dc/duW  =  1d[(c~C|)/ P] large  =  ^  =  {  2Q 
(dc/du)smaji  ln[(c-Ci)/p]smaU  ln254 

Thus  the  resulting  from  the  integration  of  stationary  crack  fields  (2]  increasing  crack 
advance  per  unit  displacement  (associated  with  the  strain  distribution  flattening  out 
in  front  of  the  crack  at  a  decreasing  rate)  leads  to  larger  size  effects  than  those 
experimentally  observed.  Notice,  however,  that  a  solution  based  on  a  superposition 
of  stationary  singularities  does  not  take  into  account  the  hardening  of  the  material 
left  behind  the  growing  crack.  More  realistically,  a  corresponding  fully-plastic, 
incremental  plasticity  solution  should  be  obtained  for  a  growing  mixed  mode  crack. 

To  study  the  directional  effects,  an  incremental  solution  was  developed  in 
chapter  2.  The  far  field  displacement  was  again  assumed  to  be  along  the  shear  band. 
At  the  initiation  and  at  each  growth  step  several  sites  around  the  current  crack  tip 
were  considered  and  the  crack  was  assumed  to  advance  to  the  direction  requiring  the 
least  far-field  displacement  to  reach  critical  damage.  The  program  predicted  that  a 
smaller  strain-hardening  coefficient  would  cause  the  crack  to  grow  closer  to  the  shear 
band  and  this  was  confirmed  from  the  experimental  results.  It  also  gave  6-15% 
higher  crack  growth  rates  (the  larger  increase  with  less  hardening)  than  the  pure 
Mode  II  [2]  solution  and  thus  closer  to  the  experimental  findings.  The  effect  of  strain 
hardening  was  very  small,  although  it  was  correctly  found  that  a  lower  strain 
hardening  increases  the  crack  growth  rate.  The  45°  shear  band  gave  however  a 
crack  angle  of  21°  from  the  transverse.  The  experimentally  found  angle  of 
approximately  40°  from  the  transverse  can  be  obtained  by  assuming  a  65°  shear 
band.  Notice  however  that  the  displacement  vector  angles  (Table  3)  suggest  that  we 
cannot  assume  the  far-field  displacement  taking  place  along  a  45°  shear  band  as  this 
model  did. 


CONCLUSIONS 


Id  asymmetrical  configurations  with  only  a  single  shear  band,  (which  can  occur 
with  cracks  near  welds  for  example),  the  crack  progresses  into  prestrained  material 
instead  of  the  new  material  between  the  two  shear  bands  of  the  symmetric  case. 
Experiments  on  six  alloys  have  shown  that  the  resulting  reduction  in  ductility  is 
primarily  dependent  on  the  strain  hardening  exponent.  In  the  lower  hardening  alloys 
the  crack  ductility,  defined  as  the  minimum  displacement  per  unit  ligament 
reduction,  in  the  asymmetric  case- is  less  that  a  third  that  of  the  symmetric  one  but 
in  the  higher  hardening  alloys  the  reduction  is  no  more  than  20%.  The  high  crack 
growth  rate  of  the  asymmetric  configuration  leads  also  to  correspondingly  higher 
stiffness  requirements  for  fracture-stable  design.  The  initiation  displacement  is  not 
much  different  and  a  fair  amount  of  blunting  was  observed  during  initiation  for  both 
cases.  The  crack  growth  direction  is  38°-41°  from  the  transverse  (instead  of  45°)  as 
expected  from  triaxiality,  the  higher  angles  with  the  smaller  strain  hardening.  The 
displacement  vector  is  at  about  51°-63°  from  the  transverse.  Angles  greater  than  45° 
suggest  a  Mode  I  component,  even  with  asymmetry. 
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TABLE  1 


Ambient  temperature  mechanical  properties  of  the  six  alloys  tested. 


Yield 

strength 

Tensile 

Strength 

(T.S.) 

Hardness 

HBN 

Reduction 

in  area 

MN/m2 

Alloy 

MN/m2 

kgf/mm2 

Percent 

1018  steel 

cold  finished 
0.15-0.20*  C, 
0.60-0.90*  Mn 

586 

600 

157 

49.3 

1018  steel 

normalized 
at  1700°F 
in  argon  flow 

321 

355 

101 

61.7 

A36  steel 

hot  rolled 

0.29*  C  max, 
0.60-0.90*  Mn 

281 

348 

105 

61.1 

HY80  steel 

0.18*  C.  2-3.25* 
0.10-0.40*  Mn, 
0.15-0.35*  Si 

587 

Ni. 

692 

175 

69.9 

HY100  steel 

693 

772 

195 

68.6 

0.20*  C.  2.25-3.50*  Ni, 
0.10-0.40*  Mn. 
0.15-0.35*  Si 


5086-H1 1 1 

aluminum 

4*  Mg,  0.4*  Mn. 


210 


264 


70 


45.8 


TABLE  2 


Stress-strain  equation  parameters . 


Alloy 

1018  steel 

cold  finished 

796 

0.05152 

1018  steel 

normalized 

818 

0.01716 

A36  steel 

hot  rolled 

697 

0 . 02628 

HY80  steel 

1107 

0.00702 

HY100  steel 

1180 

0.00488 

5086-Hill 

aluminum 

589 

0.00554 

! 

i 


I 


I 


n 

0.10 

0.23 

0.24 

0.12 

0.10 

0.19 


C£> 

f 

L_ 
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TABLE  3  -  TEST  RESULTS  (Ligament  lQ= 2.54  mm) 

Alloy  1018  CF  HY80  HY100  1018  norm.  A36  HR  5086-Hill 

INITIATION 

Idealized  initiation  displacement,  u^/Zq  (Fig.  4) 


m 

Sym  0.072 

0.108 

0.083 

0.348 

0.172 

0.179 

Asym  0.073 

0.110 

0.100 

0.252 

0.206 

0.161 

DUCTILITY  MEASURES 

Crack  Growth  Ductility,  D  . 

eq.  (3).  ~(du  /d/)min 

Sym  0 . 233 

0.320 

0.354 

0.258 

0.192 

0.166 

Asym  0 . 072 

0.096 

0.105 

0.215 

0.181 

0.108 

Modified  Tearing  Modulus,  T* 
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TABLE  4  -  RESULTS  FROM  LARGE  SPECIMENS  (Ligament  /  =7 . 62  mm) 


5086-Hill  Aluminum 


INITIATION 

Idealized  initiation  displacement, 

Sym  0.081 

Asym  0 . 072 


DUCTILITY  MEASURES 
Crack  Ductility,  D 
Sym 
Asym 


0.165 

0.105 


Modified  Tearing,  Modulus.  T" 

Sym 

Asym 


43.4 

13.8 


Min.  extension  rate.  Dext.  for  25  mm  gauge  length 
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Asym  0 . 080 

DISPLACEMENTS  from  fracture  profiles 
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Growth  Displ.,  u  //Q. 
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Displacement  vector  angle,  <f> 
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ue  Stress-  True  Plastic  Strain  for  HY-80,  HY-100  and 

.36  Hot  Rolled  Steel .  Dashed  lines  indicate  Bridgman  necking 

orrection  from  the  semiempirical  relations. 


ure  2b  True  Stress-  True  Plastic  Strain  for  1018  Cold  Finished. 

1018  Normalized  Steel  and  5086-Hill  Aluminum  Dashed  lines 
indicate  Bridgman  necking  correction  from  the 
semiempincal  relations. 


Figure  3a  Machining  for  precracking  of 


Second  Machining  (after  fatigue  precracking)  for  the 
Asynnetric  Specimens;  a  is  the  fatigue  crack 
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Syuetric  Specimens 
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Schematic  of  the  Fracture  Surface  Profil 
for  the  Asymmetric  and  Symmetric  cases 
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Figure  7a  Test  Data  for  the  HY-100  Steel  Asymmetric  Specimens 
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Figure  8a  Test  Data  for  the  A36  Hot  Rolled  Steel 
Asymmetric  Specimens 


est  Data  for  the  A36  Hot  Rolled  Steel 
Symmetric  Specimens. 
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Profile  for  the  HY-100  Steel 
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Fracture  Surface  Profile  for  the  HY-100  Steel 
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Figure  11a  Test  Data  for  the  38  1  mm  dia.  5086-Hill 

asymmetric  specimens  showing  the  unloading-loading 
points  for  marking  the  crack  front 
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Figure  lib  Test  Data  for  the  38  1  am  dia.  5086-Hill 

symmetric  specimens  showing  the  unloading-loading 
points  for  marking  the  crack  front. 
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Figure  12  Crack  advance-displacement  data  for  the  specimens  of 
Figs  lla.  lib  The  fatigue  marks  provided  the 
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APPENDIX 

FRACTOGRAPHIC  OBSERVATIONS  IN  ASYMMETRIC  AND  SYMMETRIC 
FULLY  PLASTIC  SPECIMENS. 

Observations  of  ductile  fracture  suggest  that  it  results  from  a  multi-step 
process  initiated  by  the  cracking  of  inclusions  or  the  separation  of  inclusion-metal 
interfaces,  followed  by  void  growth  and  coalescence.  The  coalescence  has  been 
observed  to  occur  on  a  plane  of  high  shear  stress,  giving  elongated  dimples  form,  or 
on  a  plane  normal  to  the  direction  of  maximum  tensile  stress,  giving  equiaxed 
dimples  [1].  Furthermore,  fracture  surfaces  have  been  studied  to  identify  and 
classify  the  characteristic  surface  markings  that  are  produced  by  the  deformation 
mechanisms  [2]. 

Tests  on  symmetric  and  asymmetric  specimens  were  performed  on  six  alloys  for 
which  X-ray  spectrography  gave  the  predominant  inclusions:  1018  cold  finished  steel 
with  Si-bearing  inclusions,  1018  steel  normalized  at  1700°,  A36  hot  rolled  steel  with 
MnS  inclusions,  HY80  steel  with  Al-bearing  inclusions,  HY100  steel  with  MnS 
inclusions  and  5086-Hlll  aluminum  with  F e-bearing  inclusions.  These  alloys  can  be 
separated  into  the  lower  hardening  ones  (1018  cold  finished,  HY80  and  HY100  steel) 
and  the  higher  hardening  ones  (A36  hot  rolled,  1018  normalized  steel).  In  this  work 
the  microscopic  features  of  the  fracture  surface  for  the  two  geometries  are 
quantitatively  compared. 

In  general,  for  a  given  crack  tip  opening  displacement,  the  amount  of  crack 
extension  can  be  separated  into  two  components:  a  sliding  off  component  and  a 
fracture  component.  To  quantify  the  ductility,  as  observed  from  the  fractographs, 
an  “apparent  crack  ductility”,  DAC,  observed  fractographically,  can  be  defined  as 
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the  ratio  of  that  part  of  the  projected  crack  area  exposed  by  pure  plastic  flow  to  the 
total  projected  area,  including  that  exposed  by  fracture.  For  instance,  with  n 
parabolic  dimple  markings  per  unit  area,  each  having  tip  radius  r,  the  apparent 
crack  ductility  may  be  found  by  assuming  that  the  area  jit2  of  each  parabola  opens 
up  before  arrival  of  the  crack  front,  and  the  balance  of  the  surface  is  formed  by 
sliding  off.  Then  DAC  =  1  -  njrr2.  Due  to  the  difficulty  in  measuring  these 
quantities,  only  rough  approximations  for  DAC  can  be  obtained.  Table  1  shows  these 
approximate  findings  (estimated  from  the  lower  flank  fractographs,  surface  normal  to 
the  beam)  for  the  lower  hardening  HY100  steel  and  the  higher  hardening  1018 
normalized  steel  from  the  asymmetric  and  symmetric  specimens  which  are  more 
ductile.  These  results  are  another  manifestation  of  the  fact  that  higher  hardening 
alloys  are  more  ductile  than  the  lower  hardening  ones  in  the  asymmetric 
configuration  but  almost  equally  ductile  in  the  symmetric  one. 

TABLE  1  -  Apparent  Crack  Ductility  DAC 

HY-100  steel  Asymmetric  0.51 

HY-100  steel  Symmetric  0.04 

1018  normalized  steel  Asymmetric  0.68 

1018  normalized  steel  Symmetric  0.67 

Fig.  1  shows  micrographs  of  the  upper  and  lower  flanks  for  5086-Hill 
aluminum  with  different  degrees  of  void  formation  and  shearing.  F racture  is  more 
“shear  type”  in  the  lower  flank,  indicating  a  larger  sliding  off  component  in  the 
crack  extension.  This  suggests  a  macro-mechanical  model  for  crack  growth  by 
combined  void  growth  and  sliding  off,  where  the  lower  flank  slides  off  along  the 
upper  slip  plane  and  the  upper  flank  slides  off  along  the  lower  slip  plane  by  a  smaller 
amount.  Thus  the  lower  flank  consists  of  a  larger  amount  of  “sheared”  material 
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than  the  upper. 

A  qualitative  understanding  of  the  differences  in  ductility  from  the 
fractographs  can  be  obtained  by  comparing  in  Fig.  2,  the  micrograph  for  the  less 
ductile  HY100  asymmetric  specimen  with  the  corresponding  one  for  the  more  ductile 
higher  hardening  1018  normalized  steel  (larger  and  less  elongated  voids). 

To  compare  the  symmetric  and  asymmetric  cases,  consider  Fig.  3  which  shows 
micrographs  of  the  low-hardening  1018  cold  finished  asymmetric  and  symmetric 
specimens.  This  alloy  shows  a  substantial  reduction  in  ductility  in  the  asymmetric 
configuration.  In  the  asymmetric  case  the  fracture  is  more  the  “shear  type”  with 
voids  elongated  in  the  direction  of  crack  growth;  in  the  symmetric  case  the  fracture 
is  more  the  “normal  type”  with  equiaxed  voids.  In  the  high-hardening  A36  hot 
rolled  steel,  with  small  differences  in  the  ductility  between  the  asymmetric  and  the 
symmetric  cases,  the  corresponding  micrographs  (Fig.  4)  are  not  much  different:  the 
fracture  in  the  asymmetric  case  is  almost  as  much  the  “normal  type”  as  in  the 
symmetric  case. 

“Zig-zagging”  of  the  fracture  surface  is  a  characteristic  of  some  symmetric 
specimens,  where  two  slip  planes  are  active  and  the  crack  grows  by  alternating 
shear.  Fig.  5  shows  this  for  the  5086-Hlll  aluminum.  The  wavy  (zig-zag)  region 
followed  the  fatigue  precrack.  In  the  end  the  fracture  turned  into  a  shear  lip. 
Symmetric  specimens  in  the  lower  hardening  alloys  often  turned  into  asymmetric 
ones,  following  only  one  slip  plane.  In  some  instances,  half  of  the  specimen  followed 
the  +45°  slip  plane  and  half  the  -45°  plane. 


In  conclusion,  fractographic  observation  of  deformation  during  crack  extension 


in  the  asymmetric  specimens  suggests  a  mechanism  by  fracture  followed  by  a 
different  amount  of  sliding  off  in  the  two  flanks.  The  usual  symmetric  case  suggests 
alternating  shear  and  fracture  and  in  some  cases  the  macroscopic  surface  is 
characterized  by  zig-zagging.  Noteworthy  features  in  the  asymmetric  specimens  are 
the  “shear  type”  fracture,  more  evident  in  the  lower  hardening  alloys  and  a  larger 
amount  of  sliding  off  in  the  lower  flank.  The  symmetric  specimens,  with  the  larger 
ductility,  show  in  turn  the  “normal  type”  fracture  with  more  equiaxed  voids  than 
the  corresponding  asymmetric  specimens. 
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Figure  1  Fracture  surface  of  5086-Hill  aluminum  asymmetric  specimen 
showing  the  difference  between  the  two  flanks 
(a)  Upper  flank,  (b)  Lower  flank  with  more  shearing 
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Figure  3  Fracture  surface  of  1018  cold  finished  steel  (lower 
hardening  alloy),  (a)  Asymmetric,  (b)  Symmetric  Distinctly 
more  "shear  type"  fracture  in  the  asymmetric,  less  ductile  c 


Figure  4  Fracture  surface  of  A36  hot  rolled  steel  (higher  hardening 
alloy),  (a)  Asymmetric.  (b)  Symmetric,  without  appreciable 
difference  Both  cases  have  almost  the  same  ductility 
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:igure  5 


Fracture  surface  of  5086-Hill  aluminum  symmetric  specimen 
seeding  the  macroscopic  roughness  of  the  specimen  After  the 
fatigue  crack  a  ‘  '  *avy  ’  ‘  region  [shown  better  m  (b)]  and  finai’y 


a  shear  lip  at  the  end 
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APPENDIX 

ON  FRACTURE  CHARACTERIZATION 

In  this  section  the  concepts  that  are  commonly  used  in  characterizing  fracture 
are  discussed. 

Growth  Resistance.  Representing  ductile  crack  propagation  has  been  based  on  the 

introduction  of  d(COD)/da  [l,2,3j  and  the  tearing  modulus  T  or  dJ/da  concept  [4]. 

In  general,  past  work  has  implied  that  the  results  are  all  charactrerized  by  a  common 

triaxiality  (or  are  independent  of  it)  whereas  in  fact  both  d(COD)/da  and  T  should 

depend  on  triaxiality,  since  both  cleavage  and  hole  growth  do.  Analogous  to  the 

above  measures  of  growth  resistance  are  the  previously  defined  crack  ductility 

Dg=ducp/d/  or  gauge  displacement  per  unit  reduction  in  ligament  Dext=du/d/,  and 

the  modified  tearing  modulus  T*  to  include  asymmetric  cracks.  Table  1  gives  some 

values  of  d(COD)/da  and  T  for  some  common  tests.  They  can  be  compared  with  the 

much  lower  values  of  about  0.010  for  D„  and  0.060  for  D  .  and  of  about  15  for  T* 

g  ext 

found  in  the  asymmetric  low  hardening  tests.  Symmetric  tests,  on  the  other  hand, 
show  values  of  Dg  and  Dext  close  to  those  in  Table  1.  It  may  also  be  that  an 
asymmetry,  introduced  in  the  bending  specimens  could  result  in  values  of  d(COD)/da 
below  those  of  the  symmetric  bending  specimens. 

J-controlled  growth.  In  large-scale  yielding  the  HRR  singularity  is  embedded  in  a 
plastic  zone  that  extends  throughout  the  remaining  ligament.  J-control  depends  on 
material  and  crack  geometry  (McMeeking  and  Parks  [7]).  The  finite  element  study 
of  the  asymmetric  specimens  showed  that  stress  and  strain  fields  are  consistent  with 
the  HRR  singularity  at  initiation.  The  COA  (crack  opening  angle)  concept  has  also 
been  alternatively  used  to  characterize  growth.  Shih  [1]  in  his  experimental  study 
found  that  the  COA  appears  to  be  constant  over  a  larger  range  of  growth  than  the 
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tearing  modulus.  Instead  of  dJ/da,  Shih  (1]  suggested  characterizing  growth  by  COA. 
A  condition  in  terms  of  a  tearing  modulus  based  on  the  COA  (=d6/da): 

dSE 

Ts  =  - - >>  1 

da<r0 

In  other  words  the  COA,  d5/da,  must  be  large  compared  to  the  yield  strain,  cJE. 

In  similar  fashion,  the  tearing  modulus  T  should  be  much  larger  than  unity.  How 
large  T  or  must  be  for  a  J-controlled  or  COA  controlled  growth  is  yet  to  be 
explored.  COA’s  of  more  than  10°-20°,  were  reported  in  the  COA-controlled  tests  in 
[1]. 

Stability.  Stability  depends  again  on  triaxiality  and  geometry.  Paris  et  al.  [5] 
developpd  instability  relations  for  fully-plastic  (nonhardening)  conditions  including 
some  common  test  piece  configurations.  For  example,  in  the  double  edge  cracked 
strip  in  tension  the  imposed  constraint  leads  to  a  critical  value  of  T  for  instability  six 
times  that  in  the  center  cracked  strip  in  tension.  In  the  expression  for  the  3-point 
bending  case  the  remaining  ligament  size,  /,  comes  into  the  instability  criterion,  so 
that  if  /  is  small  enough  in  the  first  place  the  situation  remains  stable  throughout. 

In  conclusion,  single-test  characterization  of  crack  propagation  can  apply  only  if 
crack  extension  occurs  in  a  certain  mode  and  configuration.  Instead  of  a  single 
parameter  representation  like  d(COD)/da,  a  set  of  d(COD)/da,  each  referring  to  a 
certain  mode  and  triaxiality,  could  conceivably  describe  adequately  the  material 
resistance  in  crack  propagation.  For  instance,  asymmetric  (mixed  mode  I  and  TI) 
fully  plastic  configurations  in  low  hardening  alloys  have  been  found  less  ductile  than 
the  corresponding  symmetric  singly  grooved  unconstrained  tensile  specimens. 
Extended  work  could  invlolve  studying  the  effect  of  triaxiality  by  performing 
constrained  asymmetric  tests.  For  example,  tensile  testing  on  doubly-grooved 
specimens  with  the  asymmetry  introduced  through  varying  notch  angles  and 
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positions;  or  wedge-splitting  of  a  doubly  grooved  specimen;  or  ductile  fracture  under 
asymmetric  bending  with  the  asymmetry  introduced  not  only  by  specimen  geometry 
but  also  through  shear  loading. 
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TABLE  1 

d (COD) /da  and  T  for  some  common  tests 


Material 

Test 

d (COD) /da 

A533B  steel 

compact 

0.205 

Y.S.  =  443  MPa 

tension  4T 

T.S.  =  574  MPa 

ASTM  A471  rotor  steel 

3-pomt  bend 

Y.S  =931  MN/m2 

8  x  1  x  0  5  in 

T.S.  =  1022  MN/m2 

Notched  and  Fatigue 

precracked,  a/w 

=0  502 

Free-cutting 

3-pomt  bend 

0  300 

Mild  steel 

fatigue  precracked 

BS  4360  Grade  50 

steel 

Y  S  =  359  MN/m2 

3-point  single 

0.250 

T.S.  =  526  MN/m2 

edge  notch  bend 

Asymmetric  and  Symmetric  tensile  tests. 

1018  CF  steel 

Asymmetric 

0  072  (0 

Y  S  =  586  MN/m2 

T.S  =  612  MN/m2 

Symmetric 

0  233  (0 

HY80  steel 

Asymmetric 

0  096  (0 

Y.S.  =  587  MN/m2 

T  S  =  708  MN/m2 

Symmetric 

0  320  (0 

* 


046) 

199) 

060) 

285) 


values  in  parentheses  are  based  on  gauge  extension 


CHAPTER  FOUR 


SHEAR-BAND  CHARACTERIZATION  OF  MIXED 
MODE  I  AND  II  FULLY  PLASTIC  CRACK  GROWTH 

TABLE  OF  SYMBOLS 

f  amount  of  fracture 

s;  amount  of  slip  along  lower  slip  plane 

s^  amount  of  slip  along  upper  slip  plane 

9f  fracture  plane 

9s[  lower  slip  plane 

upper  slip  plane 

\  fracture  parameter  (=f/su> 

£  shearing  parameter  (=s^/su) 

/j  projected  lower  flank  length 

/  projected  upper  flank  length 

lQ  initial  ligament 

9 u  upper  flank  angle  from  transverse 

9l  lower  flank  angle  from  transverse 

jj  crack  opening  angle 

uy  axial  displacement 

<t>  displacement  vector  angle  from  transverse 

0  upper  back  angle 

0.  lower  back  angle 

7u  strain  for  upper  shearing 

*j[  strain  for  lower  shearing 

ABSTRACT 

Asymmetric  fully  plastic  specimens  give  higher  crack  growth  rates  and  thus 
smaller  deformation  to  fracture  than  the  corresponding  symmetric  specimens.  A 
macro-mechanical  model  of  crack  growth  by  combined  fracture  on  one  plane  and 
sliding  off  along  two  others  describes,  for  this  idealization  of  the  physical 
mechanisms,  the  ductile  crack  growth  for  the  general  mixed  mode  I,  n  case.  The 
analysis  allows  determining  the  independent  physical  parameters  (shear  and  cracking 
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directions,  relative  amounts  of  cracking  and  shearing)  in  terms  of  the  observable 
quantities  of  the  macroscopic  fracture  (flank  angles,  flank  lengths,  back  angle). 

INTRODUCTION 

In  studying  fracture  there  is  a  need  for  understanding  and  characterizing  the 
deformation  and  crack  growth  in  the  fully  plastic  range  for  both  the  usual  symmetric 
case  and  for  the  asymmetric  case  shown  in  Fig.  1.  The  asymmetric  configuration 
may  occur  near  welds  due  to  the  constraint  of  a  heat-affected  zone  or  due  to  some 
geometric  asymmetry,  such  as  near-by  shoulders.  These  cracks  exhibit  less  ductility 
than  symmetric  ones,  because  the  crack  is  advancing  into  prestrained  and  damaged 
material  rather  than  into  the  new  material  encountered  by  a  crack  advancing 
between  two  symmetric  shear  bands. 

The  nonhardening  rigid  plastic  flow  field  of  Fig.  lb  consists  of  a  single  slip  line 
at  45°.  Strain  hardening,  however,  causes  the  deformation  field  to  fan  out.  It  also 
leads  to  adding  a  Mode  I  component,  as  suggested  by  the  direction  of  the  far  field 
displacement  being  more  axial  than  45°.  To  account  for  the  presence  of  the  Mode  I 
component  and  the  spreading  out  of  the  deformation  in  the  more  general  asymmetric 
case  we  assume  two  slip  planes  at  arbitrary  angles. 

Crack  growth  is  a  mixture  of  sliding  off  and  fracture.  In  the  general  case  it 
may  be  idealized  by  assuming  cycles  of  first  sliding  off  on  the  upper  slip  plane,  then 
on  the  lower,  and  finally  fracture  on  possibly  a  third  (Fig.  2).  The  combination  of 
cracking  and  sliding  off  gives  the  two  new  surfaces  of  the  macro  fracture.  These 
define  the  crack  opening  angle  and  the  crack  direction.  In  the  symmetric  case  the 
two  slip  planes  and  the  fracture  are  symmetric.  These  ideas  will  now  be  developed 
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quantitatively,  giving  a  description  of  the  mixed  mode  ductile  crack  growth  based  on 
an  idealization  of  the  underlying  physical  mechanisms.  The  single-band  pure  Mode  II 
asymmetric  and  the  pure  Mode  I  symmetric  behavior  can  be  obtained  as  limiting 
cases. 

ANALYSIS 

Consider  lower  and  upper  slip  planes  at  angles  9g{  and  9gu  (Fig.  2).  The  upper 
crack  flank  is  formed  by  sliding  off  along  the  lower  slip  plane  through  a  distance  s^  at 
9gl,  combined  with  fracture  over  a  distance  f  at  an  angle  0 f.  The  lower  flank  is 
formed  by  sliding  off  along  the  upper  slip  plane  at  9gxj  and  the  fracture  f  at  0f.  A 
cracking  parameter  x=f/su  aQd  a  shearing  parameter  f=s^/su  can  be  defined.  As 
independent  physical  variables  consider  the  cracking  and  shearing  parameters  x  and 
the  fracture  angle  9{  and  the  slip  angles  *S/>  °su-  The  limiting  case  of  Mode  I,  with 
two  symmetric  slip  lines  corresponds  to  0^=0,  9gu—-9gl,  s/=s|J  and  the  limiting  Mode 
II  case  of  slip  on  a  single  plane  corresponds  to  Sj=0. 

Observable  quantities  that  allow  solving  for  the  above  physical  variables  turn 
out  to  be  the  angles  between  the  faces  of  the  crack  and  the  transverse  direction  0u, 

9p  the  transverse  components  of  the  crack  flank  lengths  after  complete  separation, 
normalized  with  the  initial  ligament,  IJIq,  and  the  angle  that  the  deformed 
upper  back  surface  makes  to  the  load  axis,  0  (Fig.  3).  Other  dependent  variables  of 
interest  are  the  crack  opening  angle  (COA),  the  total  axial  displacement  per  initial 
ligament  u y//Q,  the  orientation  of  the  displacement  vector  <p,  and  the  angle  that  the 
deformed  lower  back  surface  makes  to  the  load  axis  0^  (Fig.  3).  These  can  be 
deduced  from  the  analysis  and  observed  from  the  tests,  except  for  the  lower  back 
angle  which  is  suppressed  by  the  shoulder. 


The  orientations  of  the  crack  flanks  9 and  9(  from  the  transverse  direction  can 
be  found  from  Fig.  1: 

f  sin0f  4-  su  sin^su  X  ^  +  sin0su 

al  =  tan  - =  tan  - ,  (1) 

f  cos0f  +  s  cos 9  \  cos0.  +  cos 0 

,  f  sin0.  +  s,  sin0  ,  ,  x  sin0.  +  £sin0  , 

0u  =  tan*1 - - - - - —  =  tan*1 - - - —  .  (2) 

f  cos0f  4-  Sj  cos0&l  x  cos 9f  4-  £cos 9&i 

The  crack  opening  angle  is 

COA  =  u  =  9r  0a  ,  (3) 

Fig.  4  shows  the  crack  opening  angle  as  a  function  of  x  and  f  with  the  slip  angle 
difference  0SU-#S/  as  a  parameter.  An  increasing  cracking  ratio  x=f/su  or  shearing 
ratio  f=Sj/su  leads  to  a  decreasing  COA  (notice  that  the  COA  is  more  sensitive  to  x 
than  £).  A  larger  slip  angle  difference,  °sn0sl’  leads  to  a  larger  COA.  Large  slip 
angle  differences  represent  spreading  out  of  the  deformation  and  can  simulate  the 
effect  of  a  high  strain  hardening  exponent  which  has  been  found  experimentally  to 
result  in  a  bigger  crack  opening  angle. 

The  original  ligament  thickness,  lQ,  projected  onto  the  transverse  direction,  r 
reduced  to  zero  by  the  cracking  f  and  by  the  sliding  off  and  su  (Fig.  2): 

/0  =  f  cos0f  +  su  cos0su  +  s/  cos 0?/  .  (4) 

The  corresponding  axial  extension  is 

uy  =  su  sin0su  -  s,  sin0s, .  (5) 

From  (4)  and  (5)  the  deformation  ratio,  defined  as  the  total  axial  displacement  per 
initial  ligament,  which  is  a  measure  of  the  ductility,  is  found  to  be 
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uy  sin0gu  -  £  sin^s/ 

—  — 

/0  x  cos 6f  +  cos 9SU  +  £  cosdsl 

The  deformation  ratio  behaves  as  the  COA,  being  higher  for  a  lower  x  or  f  and 
smaller  for  a  larger  slip  angle  difference  0su-0sf  Fig.  5  shows  an  example  of  variation 
of  uy//Q. 

An  expression  for  the  flank  lengths  is  desirable  because  the  final  projected 
flank  lengths  per  initial  ligament  can  be  measured.  The  projected  upper  and  lower 
flank  lengths,  /u  and  are  given  by 


/u  =  f  cos0f  +  sl  cosO&l . 


/j  =  f  cos0f  +  su  cosflsu  . 


Using  (4)  and  substituting  the  expression;  for  the  cracking  and  shearing  ratios  x  and 


f  gives 

L  X  cos0,  +  £  COS0  . 

—  = - - - - - ,  •  (9) 

/Q  x  cos0f  +  £  cos08l  +  .„os0su 

/,  x  cos0f  +  COS0 

-  - - [ - - - .  (10) 

/0  X  cos 0f  +  £  cos0sl  +  cos0gu 

The  lower  flank  ratio  IJIQ  decreases  with  an  increasing  shearing  parameter  £  in 
contrast  to  the  upper  flank  ratio  /u//Q  which  increases  (Fig.  6).  The  amount  by 
which  the  upper  back  surface  is  drawn  in,  projected  along  the  transverse,  is  (Fig.  2) 


t  =  s  cos0 
u  u  su 


A  “thinning”  ratio  for  the  upper  surface,  t  JlQ,  can  be  defined  in  terms  of  the 
independent  parameters  by  noting  that 


v  •"«  —  J % 


For  the  lower  surface  similarly, 


V^o  — 1  ‘  'A  • 


Of  interest  is  also  the  orientation  of  the  displacement  vector  from  the 
transve.-se,  <f>.  In  terms  o'  the  flank  lengths  and  angles,  it  can  be  found  from 
i  (VW  tan 5.  -  (L/Iq)  tan^ 


( W  -  iW 

The  back  angle,  defined  as  the  angle  that  the  deformed  back  surface  makes  to 
the  load  axis,  can  be  observed  macroscopically.  For  the  upper  slip  line  from  Fig.  2 

ds  cos0s 


(df  sin(0su-0f)  +  ds/  sin(0su-08/)]/cosdsu  +  dsu  sin0su 
Substituting  the  expressions  for  the  cracking  and  shearing  ratios  we  can  write 


x(cos0f  tan0gu  -  sin0f)  +  £(cos0g,  tan0gu  -  sin0s,)  +  sin0 
Similarly  for  the  lower  surface 
ds,  cos0, 


[df  sin(0.r0f)  +  dsu  s'm(dsr6sJ]/cos9s[  +  ds,  sin0s 


\(cos0f  tan0s/  -  sin0f)  +  (cos0gu  tan0s/  -  sin0su)  +  £  sin0g/ 

There  are  five  independent  macroscopically  observable  parameters:  the  flank 
angles  0u,  0,,  the  projected  flank  lengths  per  initial  ligament,  /  //Q,  /,//Q,  and  the  back 
angle  /?u>  Equations  (1)  for  <?u,  (2)  for  9(,  (9)  for  IJIQ,  (10)  for  /,//Q,  and  (16)  for  /?u 
give  the  corresponding  physical  variables  £,  0f)  9  0gu,  as  described  in  the 
Appendix. 


-  "  *"  **.  •.  1*.  **.  •*,  v,.-, 


Finally,  an  expression  for  the  shear  strain  can  also  be  found.  It  can  be 
expressed  in  terms  of  the  slip  and  the  normal  separation  between  corresponding  slip 
planes.  For  the  upper  shear  band,  su, 

\  =  su  /  [f  sin(*su-0f)  +  sl  sin(0sM]  = 

=  V  [x  sin(^su-0f)  +  £  sin(0su-y]  .  (19) 

Similarly,  for  the  lower  shearing,  s{, 

7/  =  s,  /  [f  sin(0s/-0f)  +  su  sin(^r0su)]  = 

=  £  /  [x  sin(0sf0f)  +  sin(^r0su)]  .  (20) 

APPLICATION  TO  TESTS 

Tests  were  performed  on  12.7  mm  dia.  round  bars  of  six  alloys:  1018  cold 
finished,  1018  normalized,  A36  hot  rolled,  HY-80  and  HY-100  steel  and  5086-Hlll 
aluminum  in  both  the  asymmetric  and  symmetric  configurations.  The  alloys  tested 
can  be  grouped  into  the  lower  hardening  ones  (1018  cold  finished,  HY-80  and 
HY-100)  and  the  higher  hardening  ones  (A36  hot  rolled,  1018  normalized).  The  lower 
hardening  alloys  exhibited  a  significantly  lower  ductility  in  the  asymmetric 
configuration  than  the  symmetric;  the  higher  hardening  alloys  showed  only  a  small 
reduction.  The  profiles  of  the  fracture  surface  and  the  deformed  back  surface  were 
plotted  with  a  travelling  stage  microscope  to  obtain  the  flank  lengths,  the  flank 
angles,  the  back  angle,  the  displacement  to  separation  and  the  orientation  of  the 
displacement  vector. 

To  apply  the  above  model  to  the  tests,  the  projected  crack  length  ratios  /u//Q, 
/j//Q  for  the  upper  and  lower  flanks  and  the  flank  angles  0 u,  0i  were  measured  from 


the  profiles  of  the  fracture  surface.  In  addition  the  back  angle  for  the  upper  surface 
0  was  measured  from  the  microscope  plots  of  the  back  surface.  The  projected 
length  ratios  depend  on  the  strain  hardening  exponent,  being  smaller  for  a  higher 
strain  hardening.  These  quantities  were  used  in  equations  (1),  (2),  (9),  (10),  (16)  to 
yield  the  cracking  parameter  x=f/su,  the  shearing  parameter  £=s^/su,  the  slip 
angles  ds(  and  0gu  and  the  fracture  angle  The  axial  displacement  uy//Q  and  the 
orientation  of  the  displacement  vector  <t>  can  also  be  obtained  and  compared  with  the 
test  data.  Results  are  shown  in  the  tables  at  the  end  of  this  chapter.  For  the 
asymmetric  specimens  the  shearing  ratio  £  is  found  to  be  about  0.5  indicating 
shearing  in  lower  flank  about  twice  that  in  upper  flank.  SEM  fractographs  have 
confirmed  that  the  lower  flank  shows  indeed  more  “shear  type”  fracture  than  the 
upper  one  (chapter  three).  The  slip  angle  difference  0S/-0SU  is  a  measure  of  the 
spreading  out  of  deformation  and  is  found  to  be  4°-6°  in  the  high  hardening  alloys  as 
opposed  to  l°-2°  in  the  low  hardening  ones.  The  cracking  ratio  x  is  a  measure  of  the 
relative  amount  of  fracture  and  sliding  off  and  allows  defining  an  comparing  with  the 
‘  ‘apparent  crack  ductility”,  DAC,  as  the  sliding  off  to  total  area.  Thus, 

In  upper  flank  DAC  u  =  s//(f+s<)  =  l/(x/£+l) 

In  lower  flank  DAC  (  =  su  /(f+su)  =  1/(X+1) 

A  smaller  cracking  ratio  means  higher  ductility.  In  the  low  hardening  alloys  the 
cracking  ratio  x  is  smaller  in  the  symmetric  case,  whereas  in  the  higher  hardening 
alloys  it  is  about  the  same  in  both  the  symmetric  and  asymmetric  configurations. 

The  cracking  ratio  x  in  the  higher  hardening  asymmetric  specimens  is  also  smaller 
that  of  the  low  hardening  asymmetric  ones. 


The  Mode  I  symmetric  case  corresponds  to  the  limit  of  f?s^=-0su,  f?f=0,  and 


£=s^/su=l.  Results  are  also  shown  in  the  tables  at  the  end  of  the  chapter.  The 
projected  flank  ratio  and  the  flank  orientation  ^=-0^  were  measured  to 

give  the  cracking  ratio  \  and  the  slip  angle  0s.  As  dependent  variables,  the 
displacement  to  separation  uy//Q,  and  the  back  angle  /3n=-^(  can  be  found.  The 
displacement  is  more  than  twice  that  of  the  asymmetric  case  in  the  lower  hardening 
HY-100  steel.  Fig.  7  shows  the  variation  of  the  deformation  ratio  u  /lQ  and  the 
crack  opening  angle  u>  vs.  the  cracking  parameter  x  with  0s=9sn=-6sl  as  a 
parameter. 

CONCLUSIONS 

In  asymmetrical,  singly  grooved,  fully  plastic  tensile  specimens  the  crack 
progresses  into  pre-strained  material.  This  results  in  less  ductility  than  in 
symmetrical' specimens  where  the  crack  grows  into  new  material  between  two  shear 
bands.  A  macro-mechanical  model  for  crack  advance  by  sliding  off  along  two  slip 
planes  and  fracture  in  the  asymmetric  specimens  gives  the  independent  parameters 
(shear  and  cracking  directions,  relative  amounts  of  cracking  and  shearing)  in  terms  of 
the  observable  quantities  of  the  macroscopic  fracture  (flank  angles,  flank  lengths, 
back  angle).  This  two  slip  plane  model  accounts  for  the  presence  of  a  Mode  I 
component  (far  field  displacement  more  axial  than  45°)  that  was  experimentally 
confirmed  in  the  asymmetric  case.  Higher  hardening  alloys  are  found  to  exhibit  more 
thinning  of  the  ligament  (hence  smaller  projected  length  ratios),  a  larger  slip  angle 
difference,  indicating  more  fanning  out  of  the  deformation  and  a  bigger  sliding  off 
component.  The  analysis,  based  on  an  idealization  of  underlying  physical 
mechanisms,  describes  the  deformation  that  leads  to  a  larger  crack  opening  angle 
and  displacement  to  separation  in  the  higher  hardening  asymmetric  specimens 
relative  to  the  lower  hardening  ones  as  well  as  the  symmetric  specimens  relative  to 


the  asymmetric  ones. 
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APPENDIX 


The  problem  is  to  determine  X,  Z,  0f,  9st,  0SU  from  the  observable  quantities  6^ 
eu,  IJIq,  /u//q,  Pn-  The  relevant  equations  are  summarized: 

X  sin0f  +  sin0su 

6(  =  tan'1 - ,  (21) 

X  cos 0j  -(-  cos0gu 
X  sin0f  +  fsin#  . 

=  tan'1 - ,  (22) 

X  cos0f  -f  £cosOs[ 

/,  x  cos0.  +  cosfl 

-  = - 1 - - - ,  (23) 

L  x  cos0f  +  £  cos  6.  +  cos0 

£  X  cos^  +  t  cosl , 

-  = - [ - - - ,  (24) 

/Q  x  cos0f  +  Z  cosflg i  +  cos 0su 

i  cos0 

/?„  =  tan'1 - .  (25) 

X(cos0f  tan0gu  -  sin0f)  +  f(cos0s,  tan0su  -  sin0g/)  +  sin0gu 

For  convenience  define  the  upper  thinning  ratio  from  (24) 

cos* 

KK  =  *  -  ‘A  =  x  cos9f  +  { — ~  cos^  •  ,26> 

Dividing  Eqs  (21)-(25)  by  cos0gu  leaves  them  in  terms  of  four  parameters,  (A,  B,  C,  D 
defined  in  the  following)  plus  tan*  that  can  be  solved  from  the  observed  variables. 
The  first  is  found  by  dividing  (23)  by  (26): 

X  cos*f 

A  9  — —  =  yt  -  1  .  (27) 

cos(,s» 

Dividing  (23)  by  (24)  and  introducing  A  from  (27)  gives 
Z  cos  9. 

B  = - ~  =  (ljll)( A+l)-A  .  (28) 

cos*su 

Introducing  (27)  into  (21)  and  rearranging, 

X  sin*f 

C  = - =  (tan*,)(A+l)  -  tanfl  n  .  (29) 

Now  introduce  (27),  (28),  and  (29)  into  (22) 


- f kf — M  —  -  **■  -  ■<  -  j*  i  .  J  9  9  •  9*^9  *  9*  9_*  »  "  iu'  •— »  iO.  v’  r 


(30 


D  = - -  =  (tan0  )(A+B)  -  C  . 

C0S^su 

From  (25)  and  the  above  definitions  of  A,B,C,D, 


tan/?u  =  1  /  ( (A+B+l)  tan0gu  -  (C+D)  ]  . 


From  (30)  C+D  of  (31)  is  given  in  terms  of  observed  variables.  Solve  for  0 


tan^su  = 


1  +  (A+B)tan/?u  tan0u 


ou  (A+B+l)  tan/?u 

Having  found  0  ,  find  C  from  (29)  and  D  from  (30).  Now  find  9{  from  (27)  and  (29) 


tan0f  =  C/A  ,  (33 

and  0sl  from  (28)  and  (30): 

tan0g/  =  D/B  .  (34 

Then  \  is  determined  from  (27)  by  using  the  already  determined  values  of  0gu  in  (32) 
and  of  0f  in  (33) 

X  =  A  cos0su/cos0f ,  (35 


and  similarly  £  is  found  from  (28)  by  using  the  values  of  0gu  from  (32)  and  0g/  from 


£  =  B  COS0  /COS0  . 


(36 


- 
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TABLE  1 

Deformation  of  singly-grooved 
asymmetrical  specimens 


Alloy 


Observations 


HY-100  steel 
Clow  hardening,  n=0.10) 


Projected  upper  flank  0  820 

ratio,  l  /lQ 

Projected  lower  flank  0  900 

ratio,  /j//g 

Upper  flank  angle.  39® 

*u 

Lower  flank  angle,  41® 

*1 

Upper  back  angle,  /?u  14® 

Corresponding  slip  and  fracture  parameters 


Slip  angle 
Slip  angle  0gu 
Cracking  angle  6 j. 
Cracking  parameter  \ 


1018  normalized 
(high  hardening,  n=0  24) 


0.750 


Shearing  parameter  £ 


Dependent  variables 


2  912 


0  536 


deduced 


Growth  Displ.  ratio,  0  118 

u/'o 

Growth  displ  vector  56^ 

angle  from  transverse 


gauge 


0  115 


0  870 


1.518 


0  445 


deduced 


0  238 


gauge 


0.230 


Back  angle. 
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Deformation  of  singly-grooved 
asymmetrical  specimens 


Alloy  1018  CF  steel  HY-80  Steel 


Observations 

Projected  upper  flank 
ratio.  /u//Q 

0  890 

0  850 

Projected  lower  flank 
ratio.  /j//Q 

0  960 

0  930 

Upper  flank  angle. 

40° 

39° 

Lower  flank  angle. 

"i 

41° 

41° 

Upper  back  angle,  /? u 

13° 

to 

o 

Corresponding  slip  and 

fracture  parameters 

Slip  angle  0^ 

50° 

52° 

Slip  angle  0gu 

51° 

54° 

Cracking  angle  0 j. 

39° 

o 

CO 

CO 

Cracking  parameter  x 

6  335 

3  822 

Shearing  parameter  £ 

0  354 

0  440 

Dependent  variables 

deduced 

gauge 

deduced 

Growth  Displ  ratio. 

0.088 

0  084 

0  120 

V'o 

Growth  displ  vector 

51° 

52° 

56° 

Growth  displ  vector 
angle  from  transverse 
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Deformation  of  singly-grooved 
symmetrical  specimens 
Assumed  <?s=^su="^si.  £=s/su=1  •  0f0° 

Alloy  1018  CF  steel  HY-80  steel 


Observations 

Projected  flank  0.820 

ratio.  /u//0=/,//0 

Flank  angle.  9^ 

Corresponding  parameters 

Slip  angle  0g  36^ 

Cracking  parameter.  2  883 

X 


0  800 


2  204 


Dependent  variables 


Growth  Displ.  ratio, 

u/'o 

Back  angle. 


deduced 

observed 

deduced 

observed 

0  260 

0.262 

0  359 

0.362 

12° 

12° 

o 

o 

12° 

Deformation  of  singly-grooved 
asymmetrical  specimens 


Alloy :  A36 

hot  rolled  steel 

5086-Hill  aluminum 

Observations 

Projected  upper  flank 
ratio.  /U//Q 

0  770 

0  810 

Projected  lower  flank 
ratio.  /j//q 

0.890 

0  900 

Upper  flank  angle. 

CO 

C T) 

O 

41° 

Lower  flank  angle. 

*1 

41° 

39° 

Upper  back  angle.  0 

43° 

16° 

Corresponding  slip  and 

fracture  parameters 

Slip  angle  0 gj 

53° 

51° 

Slip  angle  0su 

57° 

53° 

Cracking  angle  0 j 

32° 

37° 

Cracking  parameter  x 

1  834 

2  821 

Shearing  parameter  £ 

0  425 

0  507 

Dependent  variables 

deduced  gauge 

deduced  gauge 

Growth  Displ  ratio. 

0  214  0  216 

0  126  0.13 

Vo 

Growth  displ  vector 

61°  60° 

55° 

Growth  displ  vector 
angle  from  transverse 
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Deformation  of  singly-grooved 
symmetrical  specimens 
Assumed  0s=0sl=-*sl.  ^8,/s^l.  Of  0°. 


Alloy :  A36  hot 

rolled  steel 

50S5-H111 

aluminum 

•  "* 
\~4 

Observations 

*.  < 

Projected  flank 
ratio.  /u//0=/,//0 

0  780 

0.760 

Flank  angle. 

10° 

10° 

Corresponding,  parameters 

Slip  angle  0& 

32° 

29° 

Cracking  parameter. 

X 

2  158 

1 

892 

Dependent  variables 

deduced 

observed 

deduced 

observed 

Growth  Displ.  ratio. 

V'o 

0  275 

0  254 

0  263 

0  278 

Back  angle. 

16° 

15° 

19° 

16° 

— 

-  — r 


vv> 


*  ^  •  K  *•  *'  *-%.*-* 

*  *  ?  ^  a  *  t  - 


crecK  opening  angle. 


Shearing  paran*r*r  £ 


Figure  4  Crack  opening  angle  as  a  function  of  the  cracking  parameter 
X  for  £=s/su=0  5.  and  the  shearing  parameter  £ 
for  x=f/su=2  0  In  both  cases  0f=  38°  and  0f/=4O° 


I 


V. 

•fc. 


f 


Shearing  paranefer  £ 


Figure  5  Deformation  ratio  as  a  function  of  the  cracking  parameter 
X  for  (=S;/su=0  5.  and  the  shearing  parameter  £ 
for  x=f/su=2  0  In  both  cases  0f=38°  and  0f/=4O° 


croc*  opening  angle  , 


0  12  3  + 

CracKlny  paranwler  X 


Figure  7  Deformation  ratio  and  crack  opening  angle  as  a  function 
of  the  cracking  parameter  x  for  tb®  limiting  Mode  I 
symmetric  case.  0{=O .  ^s/s^l.  0f=0su=-0ti 


CHAPTER  FIVE 


FINITE  ELEMENT  INVESTIGATION  OF  PLANE 
STRAIN  ASYMMETRIC  FULLY  PLASTIC  SPECIMENS 


TABLE  OF  SYMBOLS 

n  strain  hardening  exponent 

o,1  flow  stress  at  unit  strain 

r]  damage  (eq.  3) 

cr  mean  normal  stress 

7  principal . shear  strain 

r  principal  shear  stress 

(T  equivalent  stress 

eq  n 

f  equivalent  strain 

eq  n 

9  critical  orientation 

c 

p  mean  inclusion  spacing 

<j>  displacement  vector  angle  from  transverse 

axial  component  of  far  field  displacement  vector 
U  transverse  component  of  far  field  displacement  vector 

M^  mixity  parameter  (eq.  4) 

uy  axial  component  of  relative 

crack  tip  displacement  (=ux+_ux~) 

u^  transverse  component  of  relative 

crack  tip  displacenment  (=ux+-ux  ) 

ABSTRACT 


Crack  initiation  and  early  growth  in  asymmetric,  fully  plastic,  plane  strain 
configurations  in  power-law  hardening  materials  is  investigated  numerically  via  the 
finite  element  method.  In  such  asymmetric  configurations  a  single  shear  band  is 
present  instead  of  the  two  shear  bands  of  the  symmetric  case.  Results  for  two  strain 
hardening  exponents,  n=0.12  and  n=0.24,  indicate  that  cracking  occurs  at  an  angle 
of  39°-43°  from  the  transverse,  smaller  than  45°  due  to  the  higher  triaxiality.  The 
direction  of  cracking  is  closer  to  45°  for  lower  strain  hardening  exponents  and  is 
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within  2°  of  those  experimentally  found.  The  stress  and  strain  field  is  consistent 
with  the  power  law  singularity  of  the  HRR  fields.  The  far  field  displacement  vector 
is  not  along  the  shear  band  but  at  about  68°-70°  from  the  transverse  at  initiation, 
indicating  the  presence  of  a  Mode  I  component.  Early  growth,  studied  by  successive 
removal  of  elements  reaching  unit  damage,  results  in  crack  growth  per  unit 
displacement  for  the  lower  hardening  case  about  twice  that  of  the  higher  hardening 
one. 


INTRODUCTION 

Asymmetric  plane  strain  specimens  have  been  used  to  study  crack  growth 
along  a  single  shear  band.  Such  cases  may  occur  when  a  weld  fillet  or  a  harder  heat- 
affected  zone  on  one  side  of  the  crack  suppresses  the  other  shear  zone  that  wouid 
appear  in  a  symmetric  specimen.  Based  on  Shih’s  extension  to  mixed  mode  [1]  of  the 
HRR  [2,3]  fields,  McClintock  and  Slocum  [4]  developed  an  approximate  formulation 
for  the  accumulation  of  damage  directly  ahead  of  an  asymmetric  crack.  The  crack 
was  assumed  to  follow  the  center  of  a  45°  shear  band  and  the  far-field  displacement 
was  assumed  to  be  parallel  to  the  shear  band.  It  was  found  that  the  initiation 
displacement  was  of  the  order  of  the  fracture  process  zone  size  p.  To  study  the 
directional  effects,  several  sites  around  the  current  crack  tip  were  considered  in 
chapter  two  and  the  crack  was  assumed  to  advance  to  the  direction  requiring  the 
least  far  field  displacement  for  critical  damage.  The  far  field  displacement  vector 
was  again  assumed  to  be  at  45°  from  the  transverse.  Both  these  solutions  found  only 
little  effect  of  strain  hardening  on  the  crack  growth  rate.  However,  tests  have  show  n 
that  the  far  field  displacement  vector  is  not  at  45°  but  more  axial,  at  an  angle  of 
about  60°  from  the  transverse.  In  addition,  a  lower  strain  hardening  exponent  n  was 
found  to  increase  the  crack  growth  rate  dramatically.  Strain  hardening  causes  the 
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deformation  field  to  fan  out.  The  effect  of  the  finite  width  of  the  shear  band  can  be 
captured  with  a  finite  element  investigation.  In  the  following  the  finite  element 
method  is  used  to  study  crack  initiation  and  early  growth  in  fully  plastic  plane  strain 
asymmetric  specimens. 

TECHNIQUE 

The  finite  element  grid  used  is  indicated  in  Fig.  1,  with  the  details  of  the 
refined  mesh  for  the  first  circle  around  the  crack  tip  shown  in  Fig.  2.  An  increased 
element  concentration  near  the  45°  line  is  used  to  account  for  the  high  strain 
gradients  there.  Angular  spacings  of  3.75°  for  four  sectors,  7.5°  for  two  sectors,  15° 
for  four  sectors  and  30°  for  nine  sectors  are  used.  Minimum  radial  size  for  the  3.75° 
elements  is  p=0.01  mm,  the  approximate  value  for  the  mean  inclusion  spacing.  The 
radial  size  ratio  was  s=1.155  for  the  3.75°  sectors  becoming  s2  for  the  7.5°  sectors,  s4 
for  the  15°  and  s8  for  the  30°  sectors.  The  net  ligament  of  the  specimen  is  lQ— 2.55 
mm.  8-node  plane  strain  isoparametric  elements  are  used.  The  mesh  consisted  of  a 
total  of  207  elements  with  722  nodes  and  1444  degrees  of  freedom.  The  nodes  at  the 
bottom  were  on  rollers  with  the  center  node  pinned.  An  axial  displacement  with 
zero  shear  traction  was  applied  at  the  nodes  of  the  upper  end.  The  analysis  was 
carried  on  a  Data  General  computer  available  at  M.I.T.  and  the  general  purpose 
finite  element  code  ABAQUS  [5j  was  used. 

The  mesh  was  checked  by  comparing  the  theoretical  strain  distribution  for  the 
elastic  and  the  low  hardening  n= 1  / 13  HRR  [2,3]  fields  with  the  linear  variation  of 
strains  within  the  8-node  elements.  The  radial  variation  in  strain  showed  a 
maximum  deviation  of  15%  from  the  elastic  solution  for  the  first  element  around  the 
tip.  For  the  HRR  n=l/13  solution  the  deviation  was  33%,  dropping  to  5.6%  for  the 


second  element.  The  angular  distribution  in  showed  a  maximum  deviation  of  8% 
from  the  HRR  for  the  30°  sectors.  In  addition,  a  circular  portion  of  the  finite 
element  mesh  with  16  radially  elements  at  the  finest  sectors  was  tested  by  imposing 
Mode  I  HRR  displacement  boundary  conditions.  The  HRR  singularity  in  ctQ  was 
reproduced  with  no  more  than  5%  deviation  in  all  elements  except  the  first  one, 
where  the  maximum  deviation  was  14%  at  the  first  integration  point. 


The  material  is  modeled  as  isotropic  power-law  hardening:  the  stress  <x  is  given 
in  terms  of  the  plastic  strain  tp,  the  flow  stress  at  unit  strain  <Tj,  a  strain  hardening 
exponent  n,  and  a  pre-strain  fQ  by 

<r  =  .  (1) 

Two  cases  were  considered,  n=0.24,  ^=826  MN/m2,  yield  strength  Y=333  MN/m2 
and  n=0.12,  <^=909  MN/m2,  Y=435  MN/m2. 


The  fracture  criterion  of  McClintock,  Kaplan  and  Berg  [6]  is  used,  by  which  it 
is  postulated  that  fracture  due  to  micro-void  coalescence  occurs  when  a  quantity  ij, 
named  “damage”,  reaches  a  critical  value  of  unity.  The  damage  is  expressed  in 
terms  of  a  hole  growth  ratio  Ft,  the  principal  shear  strain  7,  and  the  triaxiality 
(defined  as  the  ratio  of  the  mean  normal  stress  <7  to  principal  shear  strain  r).  In 
terms  of  the  equivalent  stress  and  strain, 

T=<req/^’  T -*„>/».  (2) 


(3) 


Ill 


unity  at  a  point  (p,  9C)  where  p  is  the  fracture  process  zone  size  and  9q  is  the  critical 
orientation.  To  study  the  first  steps  of  crack  growth,  successive  elements  were 
removed  as  they  reached  a  damage  of  unity. 

RESULTS  AND  DISCUSSION 

The  axial  displacement,  at  the  upper  end,  Uy,  was  gradually  increased  and  the 
damage  from  (3)  was  calculated  at  each  site  around  the  tip.  Cracking  occurs  when 
the  fracture  criterion  of  17=  1  is  first  satisfied. 


The  initiation  conditions  (critical  orientation  from  the  transverse,  0  ,  critical 
strain  7  ,  critical  triaxiality  <t/t,  far-field  displacement  u J p)  are  shown  in  Table  1. 
The  crack  tip  initiation  displacement  is  of  the  order  of  the  mean  inclusion  spacing  as 
was  also  found  by  McCIintock  and  Slocum  (4)  and  in  chapter  2.  The  critical 
orientation  of  39-43°  from  the  transverse  and  the  far  field  displacement  vector 
orientation  of  about  68-70°  at  initiation  can  be  compared  with  the  values  of  38-41° 
for  the  crack  direction  and  58°-69°  for  the  displacement  vector  at  initiation  from 
tests.  The  lower  hardening  n=0.12  case  results  in  fracture  closer  to  the  shear  band, 
as  found  experimentally. 


For  a  crack  at  0=0°,  a  Mode  I  mixity  parameter  Mp  was  introduced  by  Shih 
[1],  defined  in  terms  of  the  near  field  stresses  by 


Mp  =  -  tan*1 
jr 


lim 


<Wr>  0=°) 


(4) 


|f-0  <rrj9(r,  0=0) 

The  mixity  parameter  varies  from  0  for  pure  Mode  II  to  1  for  pure  Mode  I.  This 
parameter  can  be  referred  to  either  the  initial  crack  direction  9=0°  or  the  final 
(critical)  one  9=9c,  giving  values  as  shown  in  Table  1.  Notice  that  the  above 
definition  of  the  Mode  I  mixity  is  with  respect  to  both  shearing  and  crack  advance  at 
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0=0,  for  both  the  limiting  cases  Mp=0  or  1.  In  the  problem  at  hand,  pure  Mode  I  is 
crack  advance  along  the  line  0=0  (corresponding  to  the  symmetric  case)  and  pure 
Mode  II  would  be  relative  deformation  and  crack  advance  along  the  45°  shear  band. 
Alternatively,  for  experiments  and  finite  element  studies,  a  definitition  of  a  Mode  I 
mixity  in  terms  of  the  displacement  field  is  helpful: 

2.  -i  ..  u«(r-,r) • 

Mj  =  -  tan  lim - .  (5) 

JT  *+0  ur(r,7T)  -  uf(r,-7r) 

Values  for  this  parameter  are  also  shown  in  Table  1.  Notice  that  for  the  non¬ 
hardening  rigid  plastic  pure  Mode  II  limit  with  a  single  slip  line  at  0=15°,  M,*=0.5 
for  the  crack  at  0=0°  but  Mj*=0  for  the  crack  at  0=45°. 


Fig.  3a  shows  the  angular  variation  of  the  <rrtf  stress  component.  The  curve  is 
consistent  with  the  Shih’s  [1]  curves  and  has  a  maximum  near  0=65°.  This 
compares  with  the  case  Mp=l,  n=l/l3  which  has  a  maximum  at  an  angle  near  97° 
and  the  Mp=l,  n=l/3  case  with  a  maximum  near  88°  whereas  the  case  of  Mp=0 
has  a  maximum  <rtB  at  0=0°.  Furthermore,  the  maximum  for  n=l/13,  Mp=0.82  is 
near  55°  and  for  n=l/3,  Mp=0.79  is  near  40°  [Ij.  The  0-variation  of  etB  is  shown  in 
Figs.  3b,  3c.  Of  the  two  peaks  in  e^,  the  one  for  positive  0  is  the  dominant  the  other 
peak  tending  to  vanish  during  growth  when  the  Mode  I  mixity  is  reduced  (Fig.  3d). 

The  radial  variation  of  the  equivalent  strain  for  n=0.12  (along  the  critical 
orientation)  is  shown  in  Fig.  4.  The  asymptotic  solution  for  power-law  hardening 
materials  yield  singularities  in  the  stress  and  strain  of  the  form  r'n^n+1^  and 
r-I/(n+l),  respectively.  The  agreement  between  the  theoretical  curve  and  the  finite 
element  results  is  within  b%. 

Fig.  5  shows  the  axial  displacement  of  the  upper  flank  relative  to  the  lower 


v./y.v, 


•  .  *  w  •  -  ■  -  *  .  •  _  n  -1  ,  *  4  *■  «  *  .*  .*  <%  ' ■*  -  ' •  *  *  *  -.  *  »  *  •  •  *  ■  »  •  »  •  .  ’  «*  *  *  -  • 


wn 


flank  at.  the  initiation  point.  The  components  of  the  relative  displacement  of  the 

crack  tip  u  (x=0,  y=0)  and  u  (x=0,  y=0)=CTOD  are  included  in  Table  1.  A 

x  y 

higher  CTOD  occurs  in  the  higher  hardening  case.  Figs.  6a  and  6b  show  the  angular 
variation  of  the  near  tip  displacement  field  for  the  two  hardening  exponents  n=0.24 
and  n=0.12,  along  with  the  nonhardening  limit.  The  far-field  displacement 
(displacement  at  the  upper  boundary)  components  U  and  U  at  the  initiation  point 
are  also  included  in  Table  1  together  with  the  far  field  displacement  vector 
orientation  from  the  transverse,  <i>.  The  value  of  <?=68°  to  70°,  instead  of  45°, 
indicates  that  we  cannot  consider  the  far  field  displacement  taking  place  parallel  to 
the  shear  band,  as  assumed  by  McClintock  and  Slocum  [4].  We  can  observe  that  the 
displacement  vector  at  initiation  is  more  axial  for  the  lower  hardening  case  with 
larger  Mp.  The  higher  triaxiality  for  angles  smaller  than  45°  is  the  main  reason  for 
the  cracking  direction  deviating  towards  the  transverse.  The  triaxiality  is  smaller 
for  n=0.24  because  of  the  smaller  Mode  I  mixity  Mp. 

Tests  have  shown  that,  in  the  asymmetric  case,  the  lower  hardening  alloys 
exhibit  a  maximum  crack  growth  rate  more  than  twice  that  of  the  higher  hardening 
alloys.  The  finite  element  mesh  of  Figure  l  was  used  to  study  the  early  growth. 
Crack  was  grown  by  successive  removal  of  the  most  heavily  damaged  element.  After 
initiation  and  removal  of  the  most  damaged  element,  the  far  field  displacement  is 
further  increased  until  critical  damage  17=1  occurs  in  the  next  row  of  elements.  At 
this  point  the  next  step  of  crack  growth  takes  place  by  removing  the  critical 
element.  After  growth  by  four  steps  {\.9%  of  the  ligament)  it  was  found  that  the 
average  displacement  per  unit  projected  crack  advance  Am/ Al  is  about  88°^  smaller 
for  the  lower  hardening  n=0.12  case  than  that  of  the  higher  hardening  n=0.24  case 
(Table  1).  Another  noteworthy  result  is  that  the  far  field  displacement  vector  U 
becomes  less  axial  as  the  crack  grows.  For  the  case  n=0. 12,  at  the  end  of  the  fourth 


step,  the  angle  of  the  displacement  vector  from  the  transverse  is  0=67.6°  instead  of 
the  initiation  value  of  0=69.5°.  Decreasing  0  angles  with  crack  growth  have  been 
experimentally  observed  (chapter  three).  During  these  steps  the  critical  elements 
were  at  the  same  angular  sector  and  no  appreciable  accelation  of  the  crack  was 
observed. 

CONCLUSIONS 

A  finite  element  investigation  of  fully  plastic  asymmetric  specimens  with  a 
single  slip  band,  as  might  be  encountered  near  a  weld,  has  provided  the  stress,  strain 
and  displacement  fields  around  the  tip.  Results  indicate  the  presence  of  a  large 
Mode  I  component  with  the  far  field  displacement  vector  at  initiation  not  along  the 
45°  shear  band  but  at  an  angle  about  67°  from  the  transverse.  The  initiation 
conditions  were  found  by  using  the  fracture  criterion  for  hole  growth  by  McClintock 
Kaplan  and  Berg  [6].  The  critical  direction  was  at  39-43°,  less  than  45°  from  the 
transverse,  increasing  for  a  lower  strain  hardening  exponent.  Displacement  to  crack 
initiation  is  about  twice  the  fracture  process  zone  size.  Stress  and  strain  fields  are 
consistent  with  the  solutions  for  the  mixed  mode  extended  IIRR  fields.  Early 
growth,  studied  by  successive  removal  of  the  most  damaged  element,  resulted  in 
crack  growth  rate  for  the  lower  hardening  case  about  twice  that  of  the  higher 
hardening  one.  The  angle  of  the  far-field  displacement  vector  from  the  transverse 
was  found  to  be  decreasing  with  crack  growth. 
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TABLE  1 

Results  of  the  finite  element  study 


•  1 

I" 

Initiation 

Conditions 

n  =  0  12 

Critical  angle 
from  transverse 

43.1° 

Far  field  displ 

VJp 

2.1 

components 

V* 

0.782 

Displacement 
vector-angle 
from  transverse 

<t> 

69.5° 

Displacements 

ux(x=0,  j=0)/p 

0  134 

at  crack  tip 

uy(x=0,  J=0)/p 

0  518 

i 

Principal  Shear 
Strain 

0.246 

•V 

p 

Triaxiality 

Mixity  parameter 

(t/t 

2  18 

.  • 

Mode  I  Mixity 

MP(Shih) 

defined  by  Shih 

(rel  to  0=0°) 

0.936 

(based  on  stresses)  (rel  to  0=0C) 

0  717 

Displ  based 

Mode  I  Mixity 

V 

0.752 

Early  Growth 


r*; 


0  24 

39  4° 

1.8 
0  737 

67  7° 

0  192 
0  564 

0  327 

1  995  ' 

0  927 
0.710 

0  815 


Far  field  displ  Au/A/ 

per  projected  crack 
advance  (4  steps) 


0  075 


0  143 


=  0.870  at  r//o  =  O.II3  ( 1st  integr.  point )  from  FE 


Radial  variation  of  the  equivalent  plastic  strain  vs  the 
theoretical  one  (HRR  singularity)  at  the  critical  angle 


CHAPTER  SIX 


ON  THE  FULLY  PLASTIC  FLOW  PAST  A  GROWING  ASYMMETRIC 
CRACK  AND  ITS  RELATION  TO  MACHINING  MECHANICS 


TABLE  OF  SYMBOLS 


H 

k 


ne  . 

V  er.e 


v, 

V 

11 

•i 

e 

u 


9 

8 


hardening  coefficient  (eq.  9) 
eq.  22 

stream  function 
mean  normal  stress 
displacement  rates 
strain  rates 
stress  deviators 
equivalent  strain 
equivalent  stress 

rigid  body  velocity  at  lower  flank 

rigid  body  velocity  at  upper  flank 

lower  boundary  of  deforming  region  (Fig.  2a) 

upper  boundary  of  deforming  region  (Fig.  2a) 

crack  opening  angle 

strain  at  upper  boundar 

“slip  angle"  (eq.  49) 


ABSTRACT 

A  tensile  logarithmic  singularity  in  the  mean  normal  stress  is  found  for  steady 
flow  of  rigid-plastic,  linearly  strain-hardening  material,  with  rigid  material  flowing 
past  straight  flanks.  For  cracks,  this  indicates  that  the  flanks  of  the  crack  tend  to 
deform.  For  the  machining  case  it  explains  the  tendency  for  precracking  ahead  of 
the  tool  which  contributes  to  a  built-up  edge,  or  the  formation  of  a  discontinuous 
chip.  Finally,  an  approximate  analysis  of  the  quasi-steady  integral  of  the  stationary 
crack  solution  shows  a  tendency  of  the  crack  flanks  to  form  a  cusp.  The  strains  for 
a  cusp  field  would  be  dominated  by  the  elastic-plastic  field  which  shows  instead  a 
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vertical  tangent  at  the  crack  tip. 

INTRODUCTION 

Fully  plastic  flow  before  fracture  is  desirable  even  in  structures  containing 
cracks  Such  ductility  is  reduced  if  plastic  flow  is  limited  to  one  shear  band,  for 
example,  by  a  weld  (Fig.  1).  In  such  asymmetric  mixed  mode  I  and  II  configurations, 
the  crack  accelerates  as  it  advances  into  pre-strained  and  damaged  material. 

Further  evidence  for  a  lowered  ductility  in  asymmetric  cracking  is  the  tendency  to 
form  a  shear  lip  at  the  end  of  a  cup  and  cone  fracture  in  a  tensile  test.  Non- 
hardening  plasticity  gives  a  shear  band  of  infinitesimal  thickness.  Strain  hardening, 
however,  causes  the  deformation  Held  to  fan  out,  leaving  a  finite  strain  except 
possibly  at  the  crack  tip. 

In  orthogonal  machining  the  geometry  is  similar,  with  the  cutting  tool 
progressing  steadily  below  the  plastic  zone.  Here,  again  due  to  strain  hardening, 
Christopherson  et  al.  (2j  found  that  the  plastic  zone  fans  out  over  10°-30°,  as 
opposed  to  the  single  plane  required  by  the  perfectly  plastic  solid. 

STRESS  SINGULARITY  WITH  RIGID  FLANKS 

Postulate  a  steady  flow  past  a  crack  in  rigid-plastic,  linearly  strain  hardening 
material.  The  mechanics  of  the  problem  should  determine  whether  or  not  the  crack 
tip  has  a  finite  angle.  Start  by  assuming  a  crack  of  finite  angle  u>  and  rigid  body 
velocity  of  the  material  flowing  along  the  flanks.  To  satisfy  incompressibility  assume 
a  stream  function  tp  in  polar  coordinates  r  and  9  .  Seek  the  form  of  the  stream 
function  in  the  immediate  vicinity  of  the  crack  tip  where  the  velocities  should  be 
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dr  r  dO  r  r 

Thus  the  only  component  of  strain  is  shear.  The  equivalent  strain  rate  is 

-  t  k  2-2  „  ,  Wl  Fff(0)  +  F(6) 

‘=W’  +e>  +2(r]r  -jr, — • 

The  stress  deviators  Sjj  are  found  from  the  stress-strain  relations  and  the 
equivalent  stress  cr  : 

3sii- 

€  ■  = - 1  €  . 


Since  ef  —  i&  —  0  ,  from  (4), 

sf  =  s^  =  0  and  s  ^  =  <r/\/Z  . 


Assume  the  material  is  rigid-plastic,  linearly  strain  hardening: 


I 


The  accumulated  equivalent  strain  Tis  calculated  by  integration  along  a  streamline, 
where  the  time  increment  is  expressed  in  terms  of  that  required  for  an  element  to 
traverse  an  increment  of  angle: 


’=  A-dt  =  f-rd9=  -  f — 
ioo  J  u-  -6  v 


J* oo  J  Vq  b  F 

Thus  the  equivalent  strain  is  independent  of  radius.  The  same  holds  for  the 
equivalent  stress  o,  by  (9),  and  also  for  the  shear  stress,  by  (8), 

sr0  =  sr<^)  • 


Now  turn  to  the  equilibrium  equations.  In  terms  of  the  mean  normal  stress  o, 
do  ds,  1  ds  s -s- 

— +  — + - —  +  — —  =  0  .  (12 

dr  dr  r  dff  r 

ds„-  1  da  1  5s-  2s  - 

—  +  -  —  + - 2  +  -*-0.  (13 

5r  r  dO  r  dO  r 

Introducing  (8)  to  eliminate  sfI  stf,  (11)  to  eliminate  5srfi/5r,  and  cross-differentiation 
to  eliminate  o  leads  to: 


d2srfl/d02  =  0  ,  from  which  ds^/dfl  =  const.  =  slS  g  . 


Now,  (12)  simplifies  with  (8),  and  can  be  integrated 
do  l  ds - 

—  +  -  “  °  '  <7  =  "Sr8  6  Nr/R)  +  W)  .  (15) 

dr  r  dO 

To  find  C(0),  differentiate  (13)  with  respect  to  9,  and  again  note  dsieJdr=0  from  (11) 
and  S0=O  from  (8).  Equating  the  result  to  the  second  partial  of  (15)  gives 
d2o  ds  -  d2C(0) 

^  =  -V  =  5^;CW  =  -s'#/  +  c‘s  +  c-  (16] 


Define  <r(R,0)  as  the  mean  normal  stress  at  9=0  and  a  convenient  radius 
R.  Equation  (15)  then  becomes: 


1 
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<r(r,0)  -  <A R,0)  =  -  srtf  tf(ln(r/R)  +  fi)  +  Cj0  .  (17) 

Thus,  the  assumption  of  rigid  flanks  subtending  a  finite  angle  would  require  that  the 
mean  normal  stress  at  the  crack  tip  (r-*0)  have  a  logarithmic  singularity.  Let  us 
now  complete  the  study  of  the  field  specified  by  the  stream  function  (1)  by  applying 
the  boundary  conditions  and  deriving  the  streamlines. 

Two  possible  flow  fields  are  consistent  with  the  constant  rate  of  shear  stress 
from  (14)  and  the  hardening  of  the  material  (increase  in  equivalent  stress  a  from  (8) 
as  it  flows  along  the  streamline).  The  first  field,  shown  in  Fig.  2a,  is  for  sjS  tf>0. 
From  (17)  this  model  gives  a  tensile  logarithmic  singularity  in  the  mean  normal 
stress  as  r— *0,  and  thus  the  field  will  be  called  “tensile”.  The  second  field,  shown  in 
Fig.  2b,  is  for  s^  8<0.  Here  the  singularity  in  the  mean  normal  stress  is  compressive 
and,  accordingly,  this  field  will  be  called  “compressive”.  A  compressive  singularity, 
however,  would  require  strains  of  order  unity  or  more  for  fracture.  Since  such  large 
strains  are  not  actually  observed  (1),  the  “compressive”  field  is  not  plausible  for  the 
growing  crack. 

Two  other  conceivable  fields  can  be  excluded.  A  single  band  being  split  by  the 
crack  (Fig.  2c)  would  have  shear  stresses  of  the  same  sign,  but  increasing  in 
magnitude  both  above  and  below  the  line  of  advance  due  to  increasing  strains  along 
a  streamline.  This  change  in  sign  of  s^  6  would  give  tensile  and  compressive 
singularities  adjacent  to  each  other,  and  a  discontinuity  in  normal  stress.  If  the 
shear  in  a  band  being  split  by  the  crack  were  to  change  sign,  on  the  other  hand, 
there  would  be  an  intermediate  region  below  yield,  and  the  band  would  separate  into 
two,  corresponding  to  those  of  Figs.  2a, b.  In  the  limit,  the  Mode  I  field  would  be 
approached. 


*  ~  .  '»  -N  '•  . 
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Thus  only  the  “tensile”  Held  of  Fig.  2a  remains.  From  (8)  and  (9)  for  positive 
shearing, 

dsf0  1  do  1  dr 

d6  s/5  66  v/5-  d0  ^ 

(18)  and  (14)  give,  with  sftf  g  =  const., 

1  He 

dO  H  v/3  d  O' 

Differentiating  (10)  gives: 
dr  F"  +  F 

—  = - .  (2( 

d  0  y/5? 

Introducing  (19)  into  (20)  gives  finally 

F"  +  k2F  =  0  ,  (21 


where 


k2  =  1  +  ilrM  =  1  + 


The  solution  of  (21)  is 

F(0)  =  A  cosk0  +  B  sink#  .  (23] 

Referring  to  Fig.  2a  we  denote  by  Vj,  Vu  the  (rigid  body)  velocities  at  the  lower  and 
upper  boundaries  of  the  deforming  region,  which  are  at  angles  0^  and  0u  respectively. 
Then  the  boundary  conditions  are: 
at  the  lower  boundary, 


uf  =  -V^cos 0( ,  and  by  (2),  F'(0j)  =  -VposOl , 
Uj  =  V^sin0p  and  by  (3),  F(0j)  =  -V^sinfy  ; 


similarly,  at  the  upper  boundary, 


F'(«J  =  -V„cos(»„  +  u)  , 


F(«.)  =  -Vasill(*„  +  «)  . 


Solving  (23)-  (27)  for  cj  in  terms  of  fy,  #u  gives 


uj  =  tan 


-l 


■ P - tankfy 
.  k(l  4-  Ptankfy) . 


where 


k  tan0.  +  tank#,, 

p  _ - I - —  . 

1  -  k  tan^  tank^u 

Substituting  back  into  the  boundary  conditions  gives 
Vu  k  sintfj  cosk#u  +  cos0l  sink#u 
Vj  k  sin(0u+oj)  coskfl^  +  cos(0u+w)  sinkfy 


(26) 

(27) 

(28) 

(29) 

(30) 


A  (Vu/V,)  sin(0u+w)  sinkfy  -  sink0u  sinfy 
V,  sink(0u-^) 

B  sinfy  cosk^u  -  (Vu/V^)  sin(0u+w)  coskfy 
yt  ~  sink (*.-*,)  ‘ 

Assume  now  a  critical  strain  7U  at  the  upper  boundary.  Then 

d\e/M  =  7U/(V*/)  ’ 


(31) 


(32) 


(33) 


and 

k2  =  1  +  .  (34) 

K-», 

The  streamlines  for  a  particular  example,  and  the  equation  for  the  rotation  of 
elements  are  given  in  the  Appendix. 
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According  to  Hill  [7j,  the  infinite  mean  normal  stress  by  (17)  cannot  be 
sustained  at  the  rigid  flank  and  this  will  lead  to  plastic  yielding.  For  further  insight, 
turn  to  the  approximate  superposition  of  singularities  for  stationary  cracks. 

SUPERPOSITION  OF  STATIONARY  SINGULARITIES 

Shih  [3j  solved  the  mixed  Mode  I  and  II  singular  fields  for  the  stationary  crack 
field,  extending  the  Mode  I  field  of  Hutchinson  [4j,  Rice  and  Rosengren  [5|  (HRR).  In 
terms  of  a  stress-strain  law  of  the  form  a  =  <7^°  ,  for  a  far  field  defined  by  the  path 
independent  integral  J  with  Mode  I  mixity  parameter  Mp  and  the  scalar  function 
Ii/n(Mp),  the  displacement  and  strain  components  at  r,0  for  the  fully  plastic 
parameters  of  interest  here  are  (see  e.g.  McClintock  [5]): 


Uj(r  ,0)  =  r 


rJ 


l/(n+l) 


Uj(0,  1/n,  Mp)  , 


l/(n+l) 


1/n,  Mp) 


(35) 


(36) 


L*1  I1/n(MP)  rJ 

Superposition  does  not  strictly  apply  to  (35)  for  two  reasons:  it  does  not  take  the 
convection  of  hardened  material  into  account  and  it  is  a  non-linear  relation  between 
displacement  and  J.  Qualitative  insights  may  be  obtained,  however,  by  assuming, 
following  (35),  that  the  displacement  increments  vary  with  radius  according  to 


ii.  a  rn/(n+1) , 


(37) 


and  that,  correspondingly,  the  strain  rates  vary  as 


iy  oc  r-1/<“+1> 


(38) 


For  the  non-hardening  material,  n=0,  (37)  correctly  indicates  displacement  rates 


independent  of  radius,  which  for  a  growing  crack,  integrates  to  displacements 
increasing  linearly  away  from  the  tip.  For  a  power  law  material  the  displacement 
rate  increases  as  a  fractional  root  of  the  distance  away  from  the  tip,  and  its  integral 
gives  displacements  of  the  order 

u g  =  J  ngdT  oC  r1+“/<n+1)  ,  (39) 

which  indicates  a  cusp. 

Correspondingly,  integration  of  the  strain  rates  from  (38)  with  respect  to  the 
distance  as  the  material  sweeps  by  the  tip  of  the  crack  gives  strains  varying  as 

€  oc  r"/<n+1)  .  (40) 

Thus,  due  to  convection,  the  strains  increase  continuously  behind  the  crack  tip, 
whereas  the  stationary  Shih  solution  gave  strains  that  decrease.  That  is,  at  any 
instant  during  the  integration,  the  material  behind  the  crack  is  actually  harder  than 
assumed  for  the  displacement  and  strain  rates  of  (37)  and  (38).  Therefore,  the  above 
superposition  exaggerates  any  cusp.  Furthermore,  the  nonhardening  solution  with 
flank  yielding,  for  bending  and  tensile  doubly  grooved  specimens  gives  linear 
displacement  increments  which,  when  integrated,  would  show  an  increasing  crack 
opening  angle  near  the  crack  tip  [6].  For  example,  the  normal  component  of  the 
displacement  field  at  the  flank  for  the  doubly  grooved  specimens  and  the  resulting 
from  integration  flank  shape  is  shown  in  Fig.  3. 

For  further  insight,  the  relative  dominance  of  the  field  for  a  cusp  will  be 
considered.  Rigid-plastic  flow  past  a  cusp-like  crack  would  not  exhibit  singularities 
in  the  strain  but,  instead,  in  higher  order  terms  like  strain  rates,  as  follows.  For  no 
crack  tip  opening  angle,  the  displacement  rates  are  of  the  form: 
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u;  =  0(rm)  +  rigid  body  motion  . 


(41) 


The  stream  function  can  thus  be  regarded  as  being  a  superposition  of  steady  state 
rigid  body  translation  of  the  material  past  the  crack  tip  and  a  strain  rate  singularity: 


ip  =  rm+lF(0)  -  crs'mff  . 


Then  the  velocities  are: 

1  dtp 

u,  = - =  rmF'(0)  -  c  cos9  , 

r  dr 

dip 

us  — - =  -  (m-H)rmF(0)  -f  c  sintf  . 

dr 


(42) 


(43) 

(44) 


Differentiating  the  displacement  rates  would  give  strain  rates,  and  hence  the 
equivalent  strain  rate,  of  the  order  0(rm*1).  The  accumulated  equivalent  strain  is 
found  by  integrating  along  a  streamline,  with  the  time  interval  to  traverse  a  given 
angular  increment  along  a  streamline  expressed  in  terms  of  the  tangential  velocity: 


f  =  /  —  rd0  . 


(45) 


00 


For  small  r,  uff  is  of  order  r°  from  (44)  while  t  is  of  order  r"1"1  and  so  the  integrand  in 
(45)  is  of  order  rm.  Thus 


=0(rm)  , 


(46) 


and  vanishes  for  small  r  (unless  m=0,  which  turns  out  to  be  the  nonhardening  case) 
and  the  strains  are  nonsingular.  Notice  that  no  stress-strain  relationship  has  been 
used  yet,  which  means  that  for  any  rigid-plastic  law  the  strains  are  nonsingular  for  a 
zero  crack  opening  angle.  This  field  would  thus  be  dominated  by  any  field  that 
exhibits  any  nonzero  crack  opening  angle.  For  example  the  elastic-perfectly  plastic 
field  shows  logarithmic  singularities  in  the  strains  but  drf/ar  is  unbounded  as  r— >0, 
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giving  a  vertical  tangent  at  the  crack  tip  [10].  The  strains  for  that  field  are  of  the 
order: 

7P  =  0[ln(R/r)]  +  0(r°) . 

which  goes  to  infinity  for  r— *0,  thus  dominates  the  local  strain.  The  large-scale  view 
of  a  fracture,  however,  may  look  like  a  cusp  and  since  the  strains  for  the  cusp  field 
are  0(rm),  increasing  with  r,  they  may  become  larger  at  a  sufficiently  distant  point. 

The  tendency  shown  above  for  flank  yielding  with  any  finite  opening  angle 
leads  to  the  need  for  an  exact  solution  of  the  growing  crack,  where  its  shape  is 
unknown  and  the  flank,  a  part  of  the  deforming  boundary,  is  traction-free.  Finally, 
connecting  the  steady-state  continuum  mechanics  solution  to  the  micromechanics  of 
hole  growth  would  require  a  transition  to  non-steady  analysis. 

DISCUSSION 

In  machining,  a  shear  band  with  an  undetermined  rigid-plastic  boundary  breaks 
through  to  a  free  surface.  The  problem  is  similar  to  mixed  mode  crack  growth, 
except  that  the  deformation  is  larger.  Christopberson  et  al  [2]  tried  to  assess  the 
effect  of  work  hardening  in  the  mechanics  of  orthogonal  machining.  By  modifying 
the  slip  line  equations  and  estimating  roughly  the  magnitude  of  the  added  term,  they 
pointed  out  that,  due  to  hardening,  the  hydrostatic  stress  changes  from  compressive 
at  the  free  surface  to  tensile  near  the  tool  point.  What  they  found  was  essentially 
the  qualitative  effect  of  the  logarithmic  singularity  derived  above  for  fully  plastic 
flow.  In  fact,  we  can  also  deduce  that,  for  a  certain  change  in  the  flow  strength 
between  the  chip  and  the  parent  material,  if  the  deforming  region  is  narrower,  the 
angular  change  of  the  shear  stress  (i.e.  s^  in  (14))  is  bigger  and,  consequently,  the 
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singularity  stronger,  in  accordance  with  their  observation  that  the  work-hardening 
effect  becomes  more  pronounced  as  the  plastic  zone  gets  narrower. 

It  is  worth  considering  now  the  region  of  dominance  of  the  logarithmic 
singularity  in  mean  normal  stress  that  would  characterize  the  flow  past  rigid  flanks. 
Using  typical  data  d<r/df  =  H  sa  Y  for  1020  steel  and  dc/d 0  0.8  ,  gives  from  (8) 

d<r  dr  dr 

s'“=^3?7r^Hrr138Y  ,47 

From  the  fully  plastic  flow  field  of  Prandtl  for  tension  of  grooved  plane  strain 
specimens  (see  e.g.  McClintock  (8|)  a  2.8Y  and  assuming  that  R  is  the  radius  at 
which  <r(r,0)  changes  sign,  gives  from  (17),  r/R  S3  0.1.  The  distance  R  is  within  the 
macroscopic  scale  as  is  evident  from  the  approximate  study  for  machining  field  done 
by  Christopherson  et  al.  [9j.  According  to  their  slip  line  theory  modified  to  include 
hardening  the  change  in  the  mean  normal  stress  Aa  from  the  free  surface  to  a  point 
in  the  band  is  roughly  estimated  in  terms  of  the  flow  strengths  in  work-piece  and 
chip,  Yw  and  Yc,  the  distance  s  from  the  free  surface,  and  the  local  width  t  of  the 
slipband: 

Aa  (Yc-Yw)  s/t\/3  .  (48 

In  machining  mild  steel,  Yc  may  be  40%  more  than  Yw,  so  Yc-Yw  «  0.4k.  Since  at 
the  free  surface  a  :=»  -Y/\/3",  the  mean  normal  stress  becomes  positive  at  about 
s/t=3,  which  for  10°  angular  width  happens  at  a  radius  R  approximately  1/3  the 
total  shear  band  length.  Thus  the  singularity  in  the  mean  normal  stress  dominates 
in  a  significant  region. 

Now  that  tension  has  been  shown  to  exist  near  the  tool  point,  it  is  possible 
that  brittle  (or  ductile)  fracture  may  occur  at  a  particular  history  of  stress  and  this 
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could  give  rise  to  the  characteristic  fracture  running  ahead  of  the  tool  point  and  the 
formation  of  a  built-up  edge  or  a  discontinuous  chip  [9].  In  particular,  according  to 
the  “tensile”  field,  the  maximum  strain  occurs  at  the  boundary  with  the  chip  (upper 
boundary  of  the  shear  band),  where  cracking  could  occur.  It  should  be  noted, 
however,  that  nonsteady  effects  have  not  been  considered. 

CONCLUSIONS 

A  logarithmic  tensile  singularity  in  the  mean  normal  stress  has  been  found  for 
rigid-plastic  flow  past  a  growing  crack  of  finite  angle  with  rigid  flanks  under 
combined  shear  and  tension.  Applied  to  the  machining  problem,  this  result  helps  to 
explain  the  formation  of  a  discontinuous  chip  or  the  precracking  ahead  of  the  tool. 

The  tensile  singularity  predicts  yielding  of  the  crack  flanks.  Approximate 
solutions  for  flank  yielding  give  contradictory  indications.  A  tendency  to  form  a 
cusp  has  been  found  from  an  order-of-magnitude  analysis  on  the  quasi-steady 
integral  of  the  extended  near  tip  HRR  singular  field.  That  result  indicates 
decreasing  strain  behind  the  crack  tip  and  hence  overestimates  any  cusp.  From  non¬ 
hardening  solutions  with  flank  yielding  for  bending  and  tension  of  doubly  grooved 
specimens  there  is  an  increasing  crack  opening  angle  near  the  tip.  Furthermore,  the 
strains  for  any  cusp  would  be  dominated  by  the  elastic-plastic  singularity  which 
gives  a  blunt  tip.  Thus  there  is  a  need  for  an  exact  solution  of  a  crack  growing  with 
deforming  flanks  into  strain-hardening  material. 
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APPENDIX 


Taking  an  example  from  machining  (Fig.  4)  for  7U  =  1.3,  0t  —  40°,  0U  = 
we  find  by  using  equations  (28)  -  (34): 

w  =  59°,  Vu/V,  =  0.73,  d7rfl/d0  =  7.44 


and,  for  9  in  radians, 


F (0)/V,  =  0.59  cos(2.90)  -  0.459  sin(2.90)  . 


A  streamline,  resulting  from  (1)  for  this  particular  example,  has  been  sketched  in  Fig. 
4.  The  velocity  triangle  shown  in  Fig.  2a  defines  a  “slip”  angle  0g: 

(V(l/V/)sin(j 

9 '  =  sin*1  - - - -j?  .  (49) 

*  L  [<V./V,)2  +  1  2(VU/V,)C0 Sa,]>/2  J 

For  this  particular  example  9&  =  45.06°. 

A  second  example  of  a  growing  crack  with  7u  =  0.25,  9{  =  30°,  0u  =  40°,  gives 

u  =  5°,  9&  =  35.2°,  Vu/V,  =  0.88,  d'1Jd9  =  1.43 

and 


F(6)/V,  =  0.0623  cos(  1.550)  -  0.744  sin(  1.550)  . 


Finally,  the  rotation  of  the  material  element  relative  to  that  of  the  stress  field  is 
important  in  hole  growth  and  thus  is  worth  considering.  The  rotation  of  the  element 
is 


m  2  V3r  r  r  at  ) 
and  from  (2)  and  (3) 


(50) 


<t>m  =  -  HF  +  F")  . 

2r 

while  that  of  the  stress  Field  is 


giving  a  relative  rotation 

^rel  =  -  (F  *  F")  • 
r  2r 


(51) 


(52) 


(53) 


For  F(0)  given  by  (23)  and  since  k>l  by  (22),  it  is  found  that  rotation  and  shear 


Fig.  2a.  The  flow  field  for  tension  in  the  band. 

machining  case  is  illustrated;  otherwise 
is  the  crack  opening  angle. 
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CHAPTER  SEVEN 

THE  ASYMMETRIC  (MIXED  MODE  I  AND  II) 

FULLY  PLASTIC  FRACTURE  -  OVERVIEW 

TABLE  OF  SYMBOLS 

crack  opening  angle 
crack  ductility  (eq.  1) 
apparent  crack  ductility  (upper  flank) 
apparent  crack  ductility  (lower  flank) 
modulus  of  elasticity 

amount  of  fracture  in  shear  band  model 

J-integral 

shear  yield 

initial  ligament 

projected  lower  flank  length 

projected  upper  flank  length 

mixity  parameter  (eq.  4) 

strain  hardening  exponent 

load 

amount  of  slip  along  lower  plane  in  shear  band  model 

amount  of  slip  along  upper  plane  in  shear  band  model 

tensile  strength 

tearing  modulus 

eq.  2 

eq.  2 

idealized  initiation  displacement 
initiation  displacement 
growth  displacement 
total  displacement  vector 
growth  displacement  vector 
mean  inclusion  spacing 

angle  of  total  displacement  vector  from  transverse 
fracture  strain 

crack  direction  from  transverse, 
fracture  plane 
upper  and  lower  slip  plane 
fracture  parameter  (=f/su) 
shearing  parameter  (=s//su) 

upper  and  lower  flank  angle  from  transverse 


r 


h 


upper  back  angle 


SUMMARY 

In  symmetric  singly  grooved  tensile  specimens  the  crack  advances  into  the 
relatively  undamaged  region  between  two  plastic  shear  zones.  A  crack  near  a  weld 
or  shoulder,  loaded  into  the  plastic  range,  may  have  only  a  single  shear  band,  along 
which  the  crack  grows  into  prestrained  and  damaged  material  with  less  ductility 
than  the  symmetrical  unconstrained  configurations.  In  this  chapter,  work  that  deals 
with  the  effect  of  asymmetry  in  crack  propagation  of  unconstrained  fully  plastic 
singly  grooved  tensile  specimens  is  summarized.  A  crack  growth  ductility  is  defined 
as  the  minimum  displacement  per  unit  crack  growth.  Tests  of  six  alloys  showed 
that,  for  the  low-hardening  alloys,  the  crack  ductility  in  the  asymmetric  case  is  less 
than  a  third  that  of  the  symmetric.  In  the  higher  hardening  alloys  the  crack 
ductility  in  the  asymmetric  case  is  smaller  by  a  factor  of  1.2  at  most.  A  noteworthy 
result  is  the  presence  of  a  Mode  I  opening  component  even  with  asymmetry,  as  is 
shown  by  the  far  field  displacement  vector  being  more  than  45°  from  the  transverse 
direction.  The  crack  direction  is  less  than  45°,  indicating  the  effect  of  triaxiality  on 
cracking.  A  macro-mechanical  model  for  crack  advance  by  combined  fracture  and 
sliding  off  along  two  slip  planes  gives  the  independent  physical  parameters  (cracking 
and  two  shear  directions,  relative  amounts  of  cracking  and  shearing)  in  terms  of  the 
observable  quantities  of  the  macroscopic  fracture  (flank  angles,  flank  lengths,  back 
angle)  for  both  the  symmetric  and  asymmetric  specimens.  A  finite  element  study  of 
the  asymmetric  specimens  gave  a  crack  direction  within  two  degrees  and  a  far  field 
displacement  vector  at  initiation  within  three  degrees  of  that  experimentally  found. 
Stress  and  strain  fields  indicate  the  presence  of  a  Mode  I  component.  Early  growth, 
studied  by  successive  removal  of  the  most  damaged  element,  resulted  in  crack 
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growth  rate  for  the  lower  hardening  case  about  twice  that  of  the  higher  hardening 
one. 


f  INTRODUCTION 
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In  symmetric  singly  grooved  tensile  specimens  the  crack  advances  into  the 
relatively  undamaged  region  between  two  symmetric  shear  zones.  In  the  fully  plastic 
case  these  zones  narrow  into  bands  that  traverse  the  section.  An  asymmetry, 
introduced  through  a  weld  fillet  or  a  harder,  heat-affected  zone  or  a  shoulder  on  one 
side  of  the  crack  (Fig.  1)  suppresses  one  of  the  two  shear  bands  that  would  appear  in 
a  symmetrical  specimen.  In  that  case  the  crack  advances  asymmetrically,  along  the 
remaining  active  slip  band  into  previously  damaged  material.  Thus  one  might  expect 
that  the  ductility  would  be  less  than  that  of  pure  Mode  I  unconstrained  symmetric 
case. 


£ 


Near  the  tip  of  the  growing  crack,  strain  hardening  will  cause  the  deformation 
field  to  fan  out.  For  power  law  creep  or  deformation  theory  plasticity  and  a 
stationary  crack,  the  asymptotic  stress  and  strain  distribution  may  be  found  from 
the  extended  by  Shih’s  [1]  HRR  [2,3]  fields  for  the  general  mixed  mode  case.  Notice, 
however,  that  such  a  superposition  of  stationary  singularities  does  not  take  into 
account  the  hardening  of  the  material  left  behind  the  growing  crack.  Indeed,  the 
stress  and  strain  fields  near  the  tips  of  growing  cracks  in  ductile  materials  are  known 
to  differ  from  the  stress  and  strain  state  around  stationary  cracks  in  the  same 
materials  as  is  shown  from  asymptotic  solutions  [4, 5, 6, 7],  supplemented  through 
finite  element  calculations  [8,9],  Thus,  more  accurately,  a  solution  for  the 
distribution  of  strain  increments  of  a  growing  mixed  mode  crack  should  be  used; 
however  such  a  solution  is  not  yet  available. 
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A  test  with  pure  shear  (Mode  II)  loading  was  carried  out  by  Chant  et  al.  (10]  of 
high  hardening  carbon  manganese  steel  (B.S.  1501-151-430A,  Y.S.=329  MN/m2,  T.S. 
=  490  MN/m2).  Small  specimens  were  subjected  to  both  Mode  II  and  Mode  I 
testing  but  the  ductility,  measured  by  dJ/da,  was  practically  the  same  although  the 
microscopic  features  for  the  pure  sh^ar  specimens  are  different  than  those  observed 
in  the  Mode  I  specimens. 

Representing  ductile  crack  propagation  has  been  in  general  based  on  the 
introduction  of  d(COD)/da  [11,12,13]  and  the  tearing  modulus  T  or  dJ/da  [14,15] 
concepts.  Such  single-parameter  measures  are  incomplete  since  the  triaxiality  and 
the  local  distribution  of  strain  are  affected  by  the  geometry  and  mode  of  loading. 

The  triaxiality  and  strain  distribution  in  turn  strongly  affect  the  cleavage  and  hole 
growth  mechanisms  of  crack  growth. 

The  objective  of  the  current  chapter  is  to  summarize  the  important  findings  of 
the  experimental,  analytical,  and  numerical  work  that  deals  with  the  effect  of  an 
asymmetry  in  crack  propagation  of  unconstrained  fully  plastic  singly  grooved  tensile 
configurations.  First,  approximate  solutions  based  on  the  superposition  of  stationary 
singularities  are  presented.  Next,  test  results  on  symmetric  and  asymmetric 
specimens  of  six  alloys  are  summarized,  along  with  a  method  for  quantifying  and 
representing  the  ductility.  In  addition,  a  macro-mechanical  model  of  crack  growth 
by  combined  fracture  on  one  plane  and  sliding  off  along  two  others,  describes,  for 
this  idealization  of  the  physical  mechanisms,  the  ductile  crack  growth  for  both  the 
asymmetric  and  symmetric  specimens.  To  account  for  the  effect  of  the  finite  width 
of  the  shear  band  and  study  the  stress  and  strain  fields  at  initiation,  a  finite  element 
investigation  of  the  asymmetric  specimens  is  performed.  Early  growth  is  also  studied 
by  successive  removal  of  the  critical  elements. 
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INTEGRATED  STATIONARY  SOLUTIONS 

1.  Pure  Mode  II  approximation. 

A  formulation  for  the  accumulation  of  damage  directly  ahead  of  an  asymmetric 
crack,  based  on  strain  increments  adapted  from  Shih’s  [1]  analysis  was  developed 
[16].  Strain  increments,  following  a  power  law  relation  were  integrated  both 
numerically  and  quasi-steadily.  The  crack  was  assumed  to  follow  the  center  of  a  45° 
shear  band  of  infinitesimal  width  with  the  far  field  displacement,  u,  being  parallel  to 
the  shear  band  (Fig.  2a).  The  critical  fracture  strain  is  determined  from  the  fracture 
criterion  by  McClintock,  Kaplan  and  Berg  [17]. 

The  predicted  displacement  to  crack  initiation  is  found  m  =  0{p),  of  the  order 
of  mean  incusion  spacing  p.  The  crack  growth  per  unit  displacement  was  predicted 
dc/du  =  0[  ln(c/p)  ],  i.e.  to  increase  approximately  as  the  logarithm  of  the  total 
crack  advance  per  inclusion  spacing  p.  The  growth  rate  was  found  to  be  practically 
unaffected  by  strain  hardening.  The  increasing  crack  growth  rate,  associated  with 
the  strain  distribution  flattening  out  in  front  of  the  crack  at  a  decreasing  rate  that 
does  not  reach  a  steady  state,  leads  to  size  effects  in  crack  growth. 

2.  Directional  effects. 

Due  to  the  higher  triaxiality  on  one  side,  there  is  a  tendency  for  fracture  in 
that  direction.  Thus  the  previous  pure  Mode  II  work  was  extended  to  include  sites 
at  several  angles  ahead  of  the  crack.  Far  field  displacement  is  again  assumed  to  take 
place  parallel  to  the  shear  band  (Fig.  2b).  Strain  increments  are  approximated  from 
the  mixed  mode,  power-law  elastic  solution  for  a  stationary  crack  [I]  and  used  with 
the  fracture  criterion  for  hole  growth  in  shear  bands  [17]  to  predict  the  critical 


direction.  The  crack  is  assumed  to  advance  to  advance  in  the  direction  that  requires 
the  minimum  far  field  displacement  to  reach  critical  damage. 

At  initiation,  several  sites  around  the  tip  are  considered.  The  strain  and  hence 
the  required  displacement  for  damage  of  unity  is  found.  The  critical  direction  is  that 
which  minimizes  the  required  displacement.  In  growth,  the  accumulated  damage  due 
to  initiation  and  prior  growth  is  found  and  then  the  required  increment  in  damage 
and  hence  far  field  displacement  is  determined. 

For  a  shear  band  at  45°  the  crack  progresses  at  an  angle  of  21°-30°  from  the 
transverse,  depending  on  the  strain  hardening,  indicating  the  effect  of  higher 
triaxiality.  The  crack  growth  rate  is  about  6-15%  higher  than  with  growth  along  the 
shear  band.  Lower  strain-hardening  results  in  the  final  crack  orientation  being  closer 
to  the  45°  shear  band  and  the  higher  crack  growth  rate. 

EXPERIMENTAL  STUDY 

Tests  were  performed  on  fatigue  precracked  asymmetric  (Fig.  3)  and  symmetric 
(Fig.  4)  specimens  of  six  alloys:  1018  cold  finished,  1018  normalized,  A36  hot  rolled, 
HY80,  HY100  steel,  5086-Hlll  aluminum.  The  low-hardening  alloys  are  the  1018 
CF,  HY80,  HY100  steel  (n~0.10)  and  the  high-hardening  alloys  are  the  A36  HR  and 
1018  normalized  steel  (n~0.23).  In  addition  to  the  load-displacement,  data,  the 
topographies  of  the  fracture  surfaces  were  plotted  using  a  metallurgical  microscope 
with  a  travelling  stage.  A  schematic  of  the  fracture  surface  profile  is  shown  in  Fig.  8. 
These  profiles  allow  determining  the  growth  displacement  vector  7^,  the  total 
displacement  v*j,  and  hence  the  initiation  displacement  7*  =7>1-v^,  as  well  as  the 
geometry  of  the  fracture  (flank  angles,  flank  lengths,  crack  orientation). 


1.  Quantifying  crack  initiation  and  growth. 

Initiation.  As  a  convenient  measure  of  crack  initiation  displacement  from  the 
load-displacement  curves,  define  the  “idealized  initiation  displacement”,  u ^/lQ  as  the 
normalized  extension  between  initial  elastic  loading  and  steepest  unloading  lines  at 
maximum  load  (Fig.  5).  The  normalized  form  is  used  in  the  plots  for  convenience  in 
correlating  crack  growth;  multiplied  by  lQ  it  becomes  analogous  to  the  more  familiar 
COD.  The  tests  gave: 

(u//*o)aSY  —  (ui,//o)sY  ’ 

and 

W/Whigh  n  -  (^K'/Wlow  n  •  (2) 

The  axial  component  of  the  initiation  displacement,  as  measured  from  the  profiles  of 
the  fracture  surfaces,  Uj//0,  has  the  same  behavior  as  the  previously  defined 
“idealized  initiation  displacement”,  u^/fg,  i.e.,  u^/q  is  not  different  between 
asymmetric  and  symmetric  cases.  In  addition,  for  the  high  hardening  alloys  it  is 
about  two  to  four  times  that  of  the  low  hardening  ones  for  both  geometries.  A 
noteworthy  observation  is  the  fair  amount  of  blunting  occuring  in  both  geometries 
(more  blunting  with  higher  n).  An  approximate  relation  can  also  be  observed: 

K'/y  -  (1.5-2.2HVW  ■  (3) 

In  short,  initiation  displacement  is  almost  the  same  same  in  both  asymmetric  and 
symmetric  cases;  strain  hardening  affects  initiation  in  both  symmetric  and 
asymmetric  specimens. 


Growth 
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(i)  For  a  measure  of  crack  growth  resistance,  define  the  crack  ductility,  Dg,  as 
the  minimum  displacement,  duc,  per  unit  ligament  reduction  d /.  The  displacement 
duc  is  associated  with  the  crack  opening  stretch  and  consists  of  the  gauge 
displacement,  du,  and  the  elastic  unloading  duun)  (Fig.  5).  The  ligament  reduction, 
d /,  is  approximated  from  the  relative  load  drop,  d/~(dP/Pmax)/Q.  Notice  that 
thinning  of  the  ligament  from  the  far  side  in  fully  plastic  flow  makes  the  reduction  in 
ligament  rather  than  crack  advance  the  appropriate  measure  of  load  drop.  Thus 


The  crack  ductility  Dg  is  analogous  to  d(COD)/da  and  is  related  to  the  crack 
opening  angle  (COA): 

Dg  ~  COA/cos20.,  (5) 

where  0  is  the  crack  orientation.  It  is  also  the  normalized  maximum  axial 
c 

compliance  allowed  for  stability: 

compliance  allowed  <  Dg/Q/Pmax  .  (6) 

Tests  showed  that  the  crack  ductility  of  the  asymmetric  specimens  vs.  that  of  the 
symmetric  ones  is  primarily  affected  by  strain  hardening.  For  example, 

(Dg)sY/<DgUsY  =  3  37  HY-lOO  steel  (n=0.10) 

=  1.06  A36  HR  steel  (n=0.21) 

In  short,  substantial  reduction  in  crack  growth  ductility  with  asymmetry  occurs  in 
low  hardening  alloys.  High  hardening  alloys  are  almost  as  ductile  in  the  asymmetric 
configuration  as  in  the  symmetric  one. 

(ii)  Other  possible  measures  of  growth  are  related  to  D  .  A  parameter,  T  , 

o 

_  <■) 

analogous  to  tearing  modulus,  T =(K/<T0“)(d J/da ),  can  be  defined.  By  approximating 


J  by  the  non-hardening  limit  (18], 

JASY=(T  S  /v^)uv^>  Jsy=2(T.S./n/3)u  ,  (7) 

we  can  define  T*  in  terms  of  Dg,  the  tensile  strength,  T.S.,  and  the  modulus  of 
elasticity,  E: 

T<asy=D*(E/v^)/T.S.,  T*SY=Dg(2E/N/3)/T.S.  (8) 

Tests  gave: 

(T  ASY^high  n/(T  ASY^low  n  >  3 
(T  SY^high  n  —  SyW  n 

In  conclusion,  strain  hardening  does  not  affect  the  ductility  of  symmetric  specimens; 
it  does  affect  the  ductility  of  asymmetric  specimens. 

(iii)  The  growth  displacement  as  measured  from  the  profiles  of  the  fracture 
surfaces,  ug//Q,  has  the  same  behavior  as  the  crack  ductility  Dg:  for  the  low 
hardening  alloys  it  is  smaller  in  the  asymmetric  configuration  than  the  symmetric  by 
a  factor  of  more  than  three  whereas  in  the  high  hardening  alloys  it  is  reduced  by  a 
factor  of  1.2  at  most. 

(iv)  The  displacement  vector  in  the  asymmetric  specimens  is  more  axial  than 
45°,  suggesting  a  Mode  I  component.  The  angle  from  transverse,  <j>,  is  between  51° 
for  the  1018  CF  steel  and  63°  for  the  1018  normalized  steel. 

(v)  The  crack  direction  in  the  asymmetric  specimens  is  less  than  45°,  indicating 
the  effect  of  triaxiality.  The  angle  0,  is  38°-41°  from  transverse,  larger  values  for  the 
lower  hardening  alloys. 
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In  conclusion,  the  experiments  showed  that  while  the  crack  initiation 
displacements  are  similar,  the  growth  displacement  for  the  low  hardening  alloys  in 
the  asymmetric  case  is  much  less  than  that  of  the  symmetric.  Triaxiality  on  one  side 
of  the  asymmetric  crack  diverts  it  from  45°  to  38°-41°  while  the  far  field 
displacement  vector  is  more  axial  than  45°,  at  51°-63°,  suggesting  a  Mode  I 
component,  even  with  asymmetry. 

Table  1  compares  the  experimental  findings  with  the  predictions  of  the 
integrated  stationary  crack  field  solutions.  The  initiation  displacement  is  an  order  of 
magnitude  larger  than  the  theoretical  one,  apparently  due  to  blunting.  The 
incremental  models  cannot  account  for  the  big  effect  of  strain  hardening  in  crack 
growth;  notice  that  these  models  are  based  on  a  superposition  of  stationary 
singularities  and  thus  do  not  take  the  convection  of  hardened  material  into  account. 
The  size  effects  in  fully  plastic  crack  growth  that  are  predicted  from  the  incremental 
pure  Mode  II  analysis  are  associated  with  the  transient  behavior  (increasing  crack 
growth  per  unit  displacement). 

Finally,  fractographic  observations  show  as  noteworthy  features  in  the 
asymmetric  specimens  the  “shear  type”  fracture,  more  evident  in  the  lower 
hardening  alloys  and  a  larger  amount  of  sliding  off  in  the  lower  flank.  The 
symmetric  specimens,  with  the  larger  ductility,  show  in  turn  the  “normal  type” 
fracture  with  more  equiaxed  voids  than  the  corresponding  asymmetric  specimens. 

SHEAR  BAND  CHARACTERIZATION  OF  CRACK  GROWTH 

To  provide  a  physical  basis  for  interpreting  the  test  data,  a  macromechanical 
model  for  crack  advance  by  sliding  off  and  fracture  was  developed.  The  model 
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assumes  in  the  general  mixed  mode  case  sliding  off  along  two  slip  planes  and  fracture 
on  a  third  and  gives  the  independent  parameters  (shear  and  cracking  directions, 
relative  amounts  of  cracking  and  shearing)  in  terms  of  the  observable  quantities  of 
the  macroscopic  fracture  (flank  angles,  flank  lengths,  back  angle). 

To  describe  the  development  of  deformation,  assume  cycles  of  sliding  off  on  an 
upper  slip  plane  at  0gu  through  a  distance  su,  then  on  a  lower  at  9$l  through  s,  and 
finally  fracture  at  6{  over  a  distance  f  (Fig.  7).  The  limiting  Mode  I  case  with  two 
symmetric  slip  planes  corresponds  to  *su=-*st>  *r=o°,  s  and  the  limiting  Mode 
II  single  slip  plane  case  corresponds  to  sj=0.  Thus  there  are  5  independent  physical 
variables:  the  slip  and  fracture  angles,  ^su’  dsl'  0V  the  cracking  ratio  \=f/su  and  the 
shearing  ratio  £=Sj/su.  Observable  quantities  that  allow  solving  for  the  physical 
variables  are  the  flank  angles,  ffu,  the  flank  lengths  normalized  with  the  initial 
ligament,  /,//0,  and  the  back  angle,  /?u,  defined  as  the  angle  the  deformed  upper 
back  surface  makes  to  the  load  axis.  Closed-form  expressions  are  given  in  chapter 
four. 

Examples  (HY-I00  and  1018  normalized  steel)  are  shown  in  Table  2. 

For  the  asymmetric  specimens  the  shearing  ratio  £  is  found  to  be  about  0.5 
indicating  shearing  in  lower  flank  twice  that  in  upper  flank.  SEM  fractographs  have 
confirmed  that  the  lower  flank  shows  indeed  more  “shear  type’’  fracture  than  the 
upper  one.  The  slip  angle  difference  0sfOsxl  is  a  measure  of  the  spreading  out  of 
deformation  and  is  found  to  be  in  the  high  hardening  alloys  4°-6°  as  opposed  to  l°-2° 
for  the  low  hardening  ones.  The  cracking  ratio  \  is  a  measure  of  the  relative  amount 
of  fracture  and  sliding  off  on  the  upper  surface  and  allows  comparing  with  the 
“apparent  crack  ductility’’,  DAC,  as  observed  fractographically  and  defined  as  the 
sliding  off  to  total  area  including  fracture.  Thus, 
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In  upper  flank  DAC,u  =  S/Af+S/)  =  !/(x/^+l)  •  (9) 

In  lower  flank  DAC  ,  =  su/(f+su)  =  l/(x+l)  •  (10) 

A  comparison  with  values  observed  from  SEM  fractographs  (Table  3),  shows  that  the 
values  from  the  fractographs  are  bigger  by  about  a  factor  of  two.  Considering  the 
idealization  of  the  complex  hole-crack  tip  interaction  and  the  difficulty  in  measuring 
Dac  (from  the  extent  of  void  growth)  in  the  fractographs,  the  agreement  is 
encouraging,  giving  the  right  trend  (low  hardening  alloys  less  ductile  in  the 
asymmetric  configuration  than  the  symmetric  but  high  hardening  alloys  almost 
equally  ductile  in  both  geometries). 


FINITE  ELEMENT  STUDY  OF  THE  ASYMMETRIC  SPECIMENS 


For  further  insight,  a  finite  element  study  of  the  asymmetric  specimens  is 
performed.  This  work  is  needed  to  relax  the  assumption  of  the  far  field  displacement 
being  parallel  to  shear  band  (as  was  presumed  in  the  superposition  of  stationary 
singularities),  to  account  for  the  finite  width  of  the  shear  band  and  to  allow 
describing  the  Mode  I  component  at  initiation.  In  this  finite  element  work  (mesh  is 
shown  in  Fig.  8)  blunting  was  neglected.  Besides  initiation,  early  growth  was  studied 
by  successive  removal  of  elements  reaching  unit  damage.  To  describe  the  Mode  I 
component,  use  the  Mode  I  mixity  parameter  Mp  based  on  stresses  [1]. 


Mp  =  -  tan* 

7T 


>L^ 

r-*0  ffjr,9= 


'=<>  Jl 


=9  , ) 


Results,  compared  with  the  superposition  of  stationary  singularities  and  test 
data  are  shown  in  Table  4.  Notice  the  presence  of  a  large  Mode  I  component  with 
the  far  field  displacement  vector  not  along  the  45°  shear  band  but  at  an  angle  about 
68°  from  the  transverse.  In  addition  stress  and  strain  fields  are  found  consistent 


v.v.v; >  vv  ,\v.  a 
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with  the  solutions  for  the  mixed  mode  extended  HRR  fields.  Displacement  to  crack 
initiation  is  of  the  order  of  the  fracture  process  zone  size.  The  critical  direction  is 
predicted  at  an  angle  of  39°-43°  from  the  transverse,  increasing  for  a  lower  strain 
hardening  exponent.  Finite  element  study  of  early  growth  resulted  in  extension  rate 
for  the  lower  hardening  case  about  half  that  of  the  higher  hardening  one. 

CONCLUSIONS 

1)  Low  hardening  asymmetric  specimens  are  substantially  less  ductile  than  the 
symmetric  ones.  For  the  crack  ductility,  D^,  defined  as  the  minimum  displacement 
per  ligament  reduction, 

(Dg)sY/(Dg)ASY  ^  ^  ^or  l°wer  hardening  alloys 
<  1.2  for  higher  hardening  alloys 

.  Thus,  there  is  a  significant  effect  of  strain  hardening  in  mixed  Mode  I  and  II  fully 
plastic  crack  growth. 

2)  The  initiation  displacement  is  insensitive  to  geometry;  however  it  depends 
on  strain  hardening.  Blunting  of  the  order  of  10 p  or  0.1mm  occurs  in  both 
geometries. 

3)  The  displacement  vector  is  more  axial  than  45°,  at  51°-63°  from  transverse, 
suggesting  a  Mode  I  component  even  with  asymmetry  (where  nonhardening  solutions 
give  pure  shear). 

4)  The  crack  direction  is  less  axial  than  45°,  at  38°-41°  from  transverse  (closer 
to  45°  with  less  hardening);  this  indicates  the  effect  of  triaxiality. 

5)  A  superposition  of  stationary  singularities  gives  practically  no  effect  of 
strain  hardening;  it  overestimates  the  ductility  of  low  hardening  asymmetric 


specimens. 


6)  Finite  element  study  of  the  asymmetric  specimens,  neglecting  blunting, 
predicts  at  initiation  a  critical  direction  at  about  40°,  a  far  field  displacement  at 
about  68°,  gives  HRR  consistent  stress  and  strain  fields  and  describes  the  Mode  I 
component. 

7)  Early  growth,  studied  by  successive  removal  of  critical  elements  shows  an 
effect  of  strain  hardening;  the  crack  growth  rate  for  n=0.12  was  twice  that  of 
n=0.24. 

8)  A  shear  band  model  by  fracture  and  sliding  off  on  two  planes  describes 
mixed  mode  crack  growth;  provides  a  physical  basis  of  interpreting  results. 

SIGNIFICANCE  AND  RECOMMENDATIONS 

Although  much  work  has  been  done  in  the  elastic  and  elastic-plastic  fracture 
mechanics,  less  is  known  for  fracture  under  fully  plastic  conditions.  In  structures, 
fully  plastic  flow  before  fracture  is  desirable  since  large  deformations  help  detect 
impending  fracture  as  well  as  help  stabilize  it  by  load  shifting.  Most  fracture  tests 
use  symmetric  specimens  (e.g.  bend,  compact  tensile  specimens).  An  asymmetric 
configuration,  however,  may  arise  due  to  a  weld  fillet  or  a  harder,  heat-affected  zone 
or  a  shoulder  on  one  side  of  the  crack.  The  effect  of  asymmetry  on  unconstrained 
tensile  specimens  has  been  quantified  and  analyzed  for  several  alloys.  Results 
reported  here  show  that  asymmetric  (mixed  mode  I  and  II)  fully  plastic 
configurations  in  low  hardening  alloys  may  be  less  ductile  than  the  corresponding 
symmetric  singly  grooved  tensile  specimens  by  more  than  a  factor  of  3,  increasing 
the  stiffness  requirements  for  fracture-stable  design.  In  addition  a  standard  way  of 
representing  tests  and  comparing  the  two  geometries  is  suggested.  In  the  fully 
plastic  state,  since  geometry  and  mode  of  loading  can  affect  the  triaxiality  and  the 
local  strain  fields,  the  crack  growth  ductility  will  not  be  a  single  parameter  but  a  set. 


each  referring  to  a  certain  configuration  and  triaxiality. 

Further  work  should  include  studying  the  effect  of  triaxiality  by  performing 
constrained  asymmetric  tests.  For  example,  tensile  testing  on  doubly-grooved 
specimens  with  the  asymmetry  introduced  through  varying  notch  angles  and 
positions;  or  laterally  constrained  singly-grooved  tensile  tests.  Fully  plastic  fracture 
under  high  triaxiality  could  be  studied  by  wedge-splitting  of  a  doubly  grooved 
specimen.  Another  extension  could  involve  testing  part-through  cracks  in  plates 
with  asymmetric  shoulders. 

On  the  analytical  side,  there  is  a  need  for  an  asymptotic  solution  of  mixed 
mode  growing  cracks,  coupled  with  finite  element  solutions  that  connect  far-field  and 
near-field  parameters;  and  possibly  finding  a  rigid-plastic  singularity  for  a  growing 
crack  with  deforming  flanks. 
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Approximate 
Mode  II  soln. 


TABLE  1 

Asymmetric  fully  plastic  cracks 

Superposition  Tests 

of  Shih  soln. 


Initiation  displacement: 

Uj=(0.8-0.6)/>  u.=(0.6-0.4)p 


Uj=(8-39)p 
due  to  blunting 


Crack  direction: 

9  =45°  9  =21°-30° 

c  c 

(assumed)  0C|  with  nj 

Growth  Ductility  (D  =du Jdl): 

6  #  c 

0.190-0.200  0.180-0.170 

no  effect  of  n  little  effect  of  n 


0c=38°-41° 

&c  \  with  n| 

0.215-0.072 
strong  effect  of 


Size  effects  in  fully  plastic  growth: 

1.5in.  dia.  specimens  vs.  0.5in.  dia.  specimens: 

(Dg^5in/(DgH.5in=L2° 


(Dg^0.5in^Dgh.5in 


1.04 
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TABLE  2 

Deformation  of  singly-grooved 
asymmetric  and  symmetric  specimens 


Alloy:  HY-100 

n=0.10 

1018  normalized 

n=0.24 

Observations 

Asy 

Sy 

Asy 

Sy 

Length  ratio, 

0.820 

0.780 

0.750 

0.74C 

Length  ratio,  /j//q 

0.900 

0.780 

0.870 

0.74C 

Flank  angle,  9 

39° 

-14° 

36° 

-12^ 

Flank  angle,  0j 

41° 

14° 

42° 

12° 

Back  angle, 

14° 

(13°) 

13° 

(16°) 

Corresponding  slip  and  fracture 

parameters 

Slip  angle,  0 

53° 

-41° 

52° 

-31° 

Slip  angle,  0su 

54° 

41° 

58° 

31° 

Cracking  angle,  9^ 

37° 

0° 

31° 

0° 

Shearing  ratio,  £ 

0.536 

1.00 

0.445 

1  00 

1.579 


Cracking  ratio,  \ 


2.912 


1.907 


1.518 


■-w* 


Shear  band  fractographs 

model 

HY100  steel 


L-. 

Asy 

DAC,1 

0.26 

0.51 

DAC,u 

0.16 

0.37 

Sy 

dac,i 

0.34 

0.64 

«*• 

DAC,u 

0.34 

0.64 

■ 

1018  normalized  steel 

» 

* » 
pm , 

Asy 

dac,i 

0.40 

0.68 

DAC,u 

0.23 

0.52 

if. 

Sy 

dac,i 

0.39 

0.67 

DAC,u 

0.39 

0.67 

*3. 


TABLE  4 

Comparing  with  the  finite  element  results. 


Finite  Element 


(blunting 

neglected) 


Initiation  displacement: 


Supemosition 


of  Shih  soln. 


Uj=(0.5-0.6)p 


u-=(0.6-0.4)/? 


Mode  I  component  at  initiation,  M^: 
relative  to  initial  crack  direction 


~  0.93  ~  0.50 


relative  to  final  crack  direction 
~  0.71  ~  0.25 


Displacement  vector  at  initiation: 


68l 


Crack  direction: 


0c==39°-43° 
0cf  with  n| 


Gauge  extension  rate 


0.158-0.075 


(Early  growth) 
strong  effect  of  n 


0= 45u 
(assumed) 


0c=21°-3O° 
0cf  with  nj 


(0.180-0.170) 


little  effect  of  n 


ui=(8-39)p 
due  to  blunting 


^~69°-58° 


0c=38°-41° 
0cf  with  nj, 


(o.iss-o.oeo^n 


strong  effect  of 


Second  machining  (after  fatigue  precracking)  for  the 
asymmetric  specimens,  a  is  the  fatigue  precrack 


symmetric  specimens,  a  is  the  fatigue  precrack 


Growth  Zone 
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APPENDIX  A  -  Experimental  Techniques 

In  this  section  the  experimental  techniques  used  (tensile  testing  and 
extensometer  connections,  microscopic  surface  topography  mapping,  fatigue¬ 
precracking)  are  described. 


1.  TENSILE  TESTS  WITH  MTS.  50  METRIC  TON  MACHINE 
Preliminary 

Set  Console  Power  switch  On  (Master  Control  Panel  413). 

Set  Hydraulic  Pump  switch  in  Room  1-014  ON. 

Check  that  Feedback  Selector  (Model  440-32  behind  the  442  Controller  Panel)  is  in 
LOCAL. 

Set  switch  in  PDP-11  behind  the  panel  under  the  disk  drives  to  OUTPUT. 

410  -  Digital  Function  Generator 

Set  to  LOCAL  (to  start  by  pressing  START  -  if  set  to  REMOTE  then  you  must 
start  from  RUN  in  413). 

Select  rates  of  loading  and  unloading,  e.g. 

For  20  mm  stroke  range, 

Rate  1  =  2.4E3  sec.  means  2400  secs  for  20  mm  loading 
Rate  2  =  2.4E2  sec.  means  240  secs  for  20  mm  unloading 
For  50  mm  stroke  range,  same  rates  of  loading  and  unloading  require 
Rate  1  =  6E3  sec., 

Rate  2  =  6E2  sec. 

PRESS  RAMP,  DUAL  SLOPE  for  monotonic  loading  and  unloading. 

Set  BREAKPOINT  to  NORMAL  and  BRKPNT  PERCENT  to  100. 
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442  -  Controller 

Press  STROKE  (for  stroke  control). 

Set  STROKE  range  as  desired  (e.g.  ±20  mm)  by  turning  the  RANGE  knob  behind 
the  panel. 

Set  LOAD  Range  e.g.  to  10K  by  turning  the  RANGE  knob  on  stroke  module  behind 
the  panel;  10K  means  10V  output  correspond  to  10,000  kgf=22,000  lb 
Set  SPAN  1  to  100% 

Zero  load  indication  (see  digital  indicator  channel  1  in  430  panel)  before  inserting 
specimen  by  turning  the  ZERO  knob  on  Load  module  (Model  440.21)  behind  the  442 
panel. 

Press  RETURN  TO  ZERO  (in  410)  to  zero  out  any  pre-existing  function  generator 
signal. 

Set  METER  at  DC  ERROR.  Zero  out  error  by  turning  SET  POINT  right  if  pointer 
is  right;  left  if  pointer  is  left.  If  whole  range  is  not  enough,  use  ZERO  knob  on 
stroke  module  behind  442  panel. 

Press  INTERLOCK  RESET. 


413  -  Master  Control  Panel 


Press  RESET. 

Press  HYDRAULIC  PRESSURE  LOW,  then  HIGH. 

Put  SET  POINT  in  442  to  5.0  -  Ram  will  move.  (Digital  indicator  channel  3  should 
read  0.0).  To  further  move  the  ram  use  ZERO  knob  on  stroke  module  (Model 
440.21)  behind  442  (Controller)  panel. 

Notes:  1.  Ram  moves  DOWN  when  turning  SET  POINT  right, 

2.  to  increase  stiffness  move  ram  up. 

As  you  tighten  the  locknuts  use  SET  POINT  to  relieve  any  compression:  watch 
digital  indicator  channel  1  for  load  (channel  3  is  stroke). 


La.  •» 


TO  STOP:  Press  HYDRAULIC  PRESSURE  LOW  then  HYD  OFF  on  413. 

To  LEAVE:  Set  Hydraulic  Pump  switch  in  Room  1-014  OFF. 

Set  Console  Power  OFF  on  413. 

Intermediate  Manual  Unloading  and  Reloading 
To  UNLOAD:  Press  HOLD  in  410  (holds  the  test),  then  turn  SET  POINT  in  442  left 
to  unload. 

TO  RELOAD:  Turn  SET  POINT  right. 

To  CONTINUE  with  preset  rate:  Press  HOLD  once  more. 

Moving  Crosshead 

Press  HYDRAULIC  LOW,  then  HIGH  on  413. 

Set  switch  behind  the  MTS  to  UNLOCK. 

Move  crosshead  by  turning  the  UP  or  DOWN  handle  as  you  wish. 

Set  switch  behind  the  MTS  to  LOCK. 

Press  RESET  on  413. 

MTS  Plotter 

3  channels,  X,Yl  and  Y2.  Yl  not  working 

e.g.  Using  X  channel  for  stroke,  selecting  5%  of  range/in  with  20  mm  stroke  range 
corresponds  to  1  mm/in  on  the  plot. 

Using  Y2  channel  for  load,  selecting  20%  of  range/in  with  10K  load  range 
corresponds  to  2000  kgf/in  on  the  plot. 


Stiffness  Data  of  the  MTS. 
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2.  EXTENSOMETER  CONNECTIONS  WITH  THE  VISHAY  AMPLIFIER/CONDITIONER 
Wiring  Correspondence 


AXIAL 

Extensometer 

-  Full  Bridge 

Extens 

# 

Vishay 

A 

RED 

D  out  + 

B 

GREEN 

C  exc  - 

C 

YELLOW 

A  out  - 

D 

BLUE 

B  exc  + 

E 

BLACK 

F  ground 

CAi.  A  corresponds  to  0.0275  in.  (0.6985  mm) 
TRANSVERSE  extensometer  -  Half  Bridge 


Extens . 

A 

B 

D 

E 


Vishay 
A  out  + 

C  exc  -  ,  out  - 
B  exc  + 

F  ground 


CAL  A  corresponds  to  0.055  in.  (1.397  mm) 


3.  SURFACE  TOPOGRAPHY  MAPPING  WITH  THE  BAUSCH  AND  LOMB  TRAVELLING 
MICROSCOPE. 


The  apparatus  consists  of  the  microscope,  linear  potentiometers,  stage 
extension  bar  and  assorted  rubber  bands  and  C-clamps  for  fixing  the  potentiometers 
on  the  microscope.  One  potentiometer  is  clamped  to  the  side  of  the  microscope  and 
one  to  the  travelling  stage. 

Rubber  bands  are  used  to  secure  the  stage  extension  bar  to  the  microscope 
knee.  Note  that  some  of  the  rubber  bands  go  around  the  back  of  the  microscope. 
They  serve  two  purposes;  they  keep  the  stage  extension  bar  firmly  against  the 
microscope  and  they  act  to  offset  its  weight. 

Adjust  the  potentiometer  and  the  stage  extension  bar  so  that  the  points  of  the 
potentiometer  are  reasonably  centered  on  the  bottom  of  the  cap  screws  in  the  stage 
extension  bar.  Use  rubber  bands  to  secure  the  potentiometer  ends  on  the  stage 
extension  bar  screws. 

The  electric  circuit  employs  a  V0=3V  battery  and  a  balancing  10-turn 
potentiometer  for  each  linear  pot.  The  two  green  wires  from  the  linear  pot  go  to  the 
battery  terminals  hooked  with  the  balancing  pot  as  in  Fig.  1.  The  blue  wire  from 
the  linear  pot  and  the  remaining  wire  from  the  balancing  pot  go  to  the  plotter. 
Solving  the  circuit  gives 

Rn  R. 

1  2  r,+r2  r,+r4  ° 

Zero  the  output  V1-V2  at  the  starting  point  by  using  the  balancing  10-turn 
potentiometer  (i.e.  adjust  R4  so  that  V1-V2=0). 
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To  plot  the  topography  use  the  lOx  power  on  the  microscope.  The  higher 
magnifications  require  the  objective  to  be  quite  close  to  the  surface  of  the  specimen 
and  one  could  easily  hit  the  objective  on  a  peak  of  the  viewed  surface  when  trying  to 
focus  into  a  valley.  Move  slowly  your  specimen  along  the  horizontal  axis  and  get  the 
corresponding  vertical  coordinate  by  having  the  centerline  of  the  specimen  always  in 
focus.  Notice  that  the  vertical  fine  scale  on  the  microscope  is  10  rev=l  mm. 
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4.  FATIGUE  PRE-CRACKING 

Specimen  precracking  was  done  on  the  SF-1  Fatigue  Machine,  which  is  a  fixed 
frequency  (3600  rpm)  rotating  mass  machine.  It  was  used  with  the  bending  fixtures. 
The  specimens  were  subjected  to  four-point  bending.  A  uniform  bending  moment  M 
is  applied  across  the  span  of  the  specimen,  given  in  terms  of  the  load  moment  arm  R 
(distance  between  rocker  bearings,  R  =  3  or  6  inches)  and  the  total  applied  force  P 


M  =  PR/2  . 


The  nominal  alternating  stress  *7,  can  be  found  by  using  the  moment  of  inertia 

A 

I=bh3/ 12  (b,  h  are  specimen  width,  and  thickness)  from 


M(h/2)  3PR  M 
1  ~~bh2  ”  bh2/6  ' 


Notice  that  <ra  does  not  depend  on  the  specimen  length.  In  terms  of  the  specimen 
length  between  grips,  L,  (L  =  3,  5,  6  or  8  inches)  and  the  modulus  of  elasticity  of  the 
specimen,  E,  the  amplitude  of  the  vibrating  platen  Y  is 
ML  3PR2L 

Y  =  R  —  = - —  .  (3) 

2EI  Ebh3 

The  following  restrictions  apply 

Maximum  applied  force,  P=1,000  lb 

Maximum  amplitude  of  reciprocating  platen  Y=0.4  inches 


As  a  rough  approximation  assume  tr^T.S.  Then  from  (2)  find  the  necessary 
M=(T.S.)bh2/6.  For  a  chosen  R  find  the  required  P  from  (1),  P=2M/R.  Next,  for 
the  chosen  specimen  length  L  check  that  the  resulting  Y  from  (3)  is  less  than  0.4 
inches.  Notice  that  P  should  be  less  than  1,000  lb. 


Indicative  data  for  fatigue  pre-cracking  on  the  SF-1  with  R=3  inches  are: 

1018  CF  steel  (HBN=157)  0.50”  dia.  needs  about  10,000  cycles  with  P=90  lb  to 
grow  0.050”  fatigue  crack. 

A36  steel  (HBN=105)  0.50”  dia.  needs  about  7,000  cycles,  P=80  lb  to  grow  0.050” 
fatigue  crack 

5086-Hlll  aluminum  (HBN=70)  1.50”  dia.  needs  about  12,000  cycles,  P=1,000  lb  to 
grow  0.150”  fatigue  crack. 

To  operate  the  machine  (for  details  see  instruction  manual): 

Turn  the  CONTROLLER  POWER  switch  to  STANDBY. 

Turn  on  main  power  switch  on  wall  behind  machine.  Wait  for  at  least  half  hour  to 
allow  warmup. 

Attach  tuning  weights  to  the  studs  on  either  side  of  the  orange  cage. 

Set  the  oscillating  load  P. 

Check  that  knob  of  variable  transformer  is  at  zero. 

Press  START  button.  Gradually  turn  knob  of  variable  transformer,  increasing  the 
motor  speed  to  the  extreme  100  position.  This  should  take  from  20  to  80  sec. 

Turn  CONTROLLER  POWER  switch  ON. 

To  STOP  press  the  STOP  button  and  turn  the  variable  transformer  back  to  zero. 
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C  This  program  creates  the  node  nos .  ,  element  no*. .  and 
C  the  corresponding  MFC  constraints  for  the  fan. 

C  We  have  maximum  Nr  segments  radially  and  N6  radial  lines 
C  around;  here  Nr* 3 2  and  N6*19. 

C  RQ=10  microns  (one  inclusion  spacing)  3  0.01  ms 
C  s  is  the  size  ratio  in  the  log  circular  mesh=1.155 
C  Angular  spacings  are  (larger  to  smaller)  thet.  thetl.  thet2, 

C  thet3?  here  thet*30,  thetl* IS.  thet2*7.5.  thet3=3.75  deg. 

C  Lis  the  node  no.  with  coords  x(L).  y(L);  M  is  the  element  no. 

C  with  nodes  Nodi  (M) ,  Nod2  (M) . Nod8  (M) . 

C  File  no.  15  contains  the  nodes,  no.  16  contains  the  element  data 
c  and  no.  17  contains  the  MPC  constraints. 

DIMENSION  x (15000) .  y (15000).  Nodi (15000).  phi (100). 

Nod2(15C00).  Nod3(15000),  Nod4 (15000) .  Rc(100), 

,  -  NodS (15000).  Nod6 (15000).  N0d7 (15000),  Nod8 (15000) 

I  Nr  *  32  More;  04r4 

e  Note:  Nr  should  be  a  multiple  of  8  (so  that  the  minimum  no.  r-rc. 

c  of  segments  radially,  corresp.  to  the  largest  angle,  is  Nr/8)  .  ,ir. 
c  Rmax  =  R0  (s*  *Ni — 1)  /  (s-1)  1~ 

I  R0  *  0.01 
1  S  a  1.155 
pi  *  3.14159 
thet  *  30 . *pi/180 . 
thetl  *  IS . *pl/180 . 
thet2  *  7. 5 *p 1/180 . 
thet3  *  3.75*pi/180. 

R30  *  R0*(l+s+s**2+s**3+s**4+s**5+s**6*s**7) 
s30  *  s**8 
N30  *  Nr/a 

R15  *  R0*(l+s-*-s**2+s**3) 

slS  *  s**4 

N15  *  Nr/4 

R75  *  RQ*  (1+s) 

s75  *  s*«2 

N75  *  Nr/2 

1001  FQRMAT(I5,F1Q.7,F10.7) 

c  Dtheta  for  radial  lines  1-N1  is  15  deg.;  N1-N2:  7.5  deg.; 
c  N2-N3:  3.75  deg.;  N3-N4:  7.5  deg.; 
c  N4-N5 :  15  deg.;  N5-N6 :  30  deg.;  N7:  crack  flank. 

'  Nl  »  3 
N2  *  4 
N3  *  8 
N4  *  9 
N5  *  11 
N6  *  19 
N7  *  14 

c  Input  Nodes  for  the  15  deg.  sectors 
DO  90  1*1,  Nl-1 
phia  *(I-l)*thetl 
phi (I)  *  phia 
DO  100  J*l,  N75 

R1  *  R75*  (  s75**J  -  1  )/(s75-l) 

L  *  500*1  +  4*J 
x(L)  *  Rl*cos(phia) 
y (L)  *  Rl* sin (phia) 

WRITE (15, 1001)  L.  x(L) ,  y(L) 

100  CONTINUE 
90  CONTINUE 

DO  911  1*1.  Nl-1 


phia  =  (1-1) *thetl  +  thetl/2. 

DO  922  J=l.  N15 

R1  =  R15*(  sl5**J  -  1  )/  (sl5-l) 

L  =  500*1  +  200  +  8*J 
x(L)  =  Rl*cos(phia) 
y  (L)  =  Rl*sin  (phia) 

WRITE  (15,  1001)  L,  x  (L)  ,  y(L) 

!  CONTINUE 
L  CONTINUE 

Input  Nodes  for  the  7.5  deg.  sector 
DO  82  I=N1.  N2-1 

phia  =(N1— 1) *thetl  +  (I-Nl) *thet2 
phi (I)  =  phia 
DO  101  J=l.  Nr 

R1  =  RO*  (  s**J  -  1  )/(s-l) 

L  =  500*1  ♦  2*J 
x(L)  =  Rl*cos(phia) 
y(L)  =  Rl*sin(phia) 

WRITE  (15, 1001)  L,  x  (L)  ,  y(L) 

L  CONTINUE 

CONTINUE 

DO  912  I=N1,  N2-1 
phia  =  (Nl-1) *thetl  +  (I-Nl) *thet2  +  thet2/2. 
DO  923  J=l,  N75 

R1  =  R75* (  s75**J  -  1  )/(s75-l) 

L  =  500*1  +  200  +4*J 
x  (L)  =  Rl*cos(phia) 
y(L)  =  Rl*sin  (phia) 

WRITE (15,  1001)  L,  x(L),  y(L) 

3  CONTINUE 

l  CONTINUE 

Input  nodes  for  the  3.75  deg.  sectors 
DO  83  I=N2 ,  N3 

phia  =(N1— 1) *thetl  +  (N2-N1) *thet2 
phi  (I)  =  phia 
Rip  =  0. 

DO  132  J=l,  Nr 

R1  =  RO*  (  s**J  -  1  )/(s-l) 

L  =  500*1  +  2*J 
x(L)  =  Rl*cos(phia) 
y(L)  =  Rl*sin  (phia) 

WRITE (IS. 1001)  L,  x (L) ,  y (L) 

L  =  500*1  ♦ (2* J-l) 

R2  =  (Rl+Rlp) /2 . 
x(L)  =  R2*cos(phia) 
y(L)  =  R2*sin(phia) 

WRITE (15,  1001)  L,  x (L) ,  y (L) 

Rip  =  R1 
l  CONTINUE 

CONTINUE 

DO  913  I=N2,  N3-1 

phia  =  (Nl-1) *thetl  +  (N2-N1) *thet2 
+  (I-N2) *thet3  +  thet3/2 . 

DO  924  J=l,  Nr 

R1  =  RO*  (  s**J  -  1  )/(s-l) 

L  =  500*1  +  200  +2*J 
x(L)  =  Rl*cos(phia) 
y(L)  =  Rl*sin(phla) 

WRITE (15.  1001)  L,  x (L) ,  y(L) 
i  CONTINUE 


sectors 


(N2-N1) *thet2  +  (I-N2) *thet3 


E 


913  CONTINUE 

c  Input  nodes  for  the  7.5  deg  sector 
DO  84  I=N3+1,  N4 

phia  =(N1-1) *thetl  ♦  (N2-N1) *thet2 

+  (N3-N2) *thet3  ♦  (I-N3) *thet2 
phi (I)  =  phia 
DO  133  J=l,  Nr 


)/(s-l) 


thet2/2 . 


)/(s75-l) 


R1  -  RG*  (  s**J  -  1  )/(s-l) 

L  =  500*1  +  2* J 
x(L)  =  Rl*cos  (phia) 
y(L)  =  Rl*sin(phia) 

WRITE (15, 1001)  L,  x(L).  y(L) 

3  CONTINUE 

CONTINUE 

DO  914  I=N3,  N4-1 

phia  =  (Nl-1) *thetl  +  (N2-N1) *thet2 

+  (N3-N2) *thet3  +  (I-N3) *thet2  +  thet2/2. 

DO  925  J=l,  N75 

R1  =  R75* (  s75**J  -  1  ) / (s75-l) 

L  =  500*1  +  200  -*■  4* J 
x(L)  =  Rl*cos(phia) 
y(L)  =  Rl*sin(phia) 

WRITE (IS,  1001)  L,  x (L) ,  y(L) 

5  CONTINUE 

1  CONTINUE 

Input  nodes  for  the  15  deg.  sectors 
DO  85  I=N4+1,  N5 

phia  = (Nl-1) *thetl  +  (N2-N1) *thet2 

+  (N3-N2) *thet3  >  (N4-N3) *thet2  +  (I-N4) *thetl 
phi (I)  =  phia 
DO  134  J=l.  N75 

R1  =  R75* (  s75**J  -  1  )/(s75-l) 

L  =  500*1  +  4*J 
x(L)  =  Rl*cos  (phia) 
y(L)  =  Rl*sln(phia) 

WRITE (15, 1001)  L.  x (L) ,  y(L) 

1  CONTI JUE 
CONTINUE 

DO  915  I=N4,  N5-1 

phia  =  (Nl-1) *thetl  +  (N2-N1) *thet2 
-  +  (N3-N2) *thet3  +  (N4-N3) *thet2  +  (I-N4) *thetl  + 

DO  926  J=l,  N15 

R1  =  R15* (  sl5**J  -  1  ) / (sl5-l) 

L  =  500*1  +  200  +  8* J 
x(L)  =  Rl*cos(phia) 
y(L)  =  Rl* sin  (phia) 

WRITE (15.  1001)  L.  x (L) ,  y(L) 

>  CONTINUE 

>  CONTINUE 

Input  nodes  for  the  30  deg.  sectors 
DO  86  I=N5+1,  N6 

phia  = (Nl-1) *thetl  +  (N2-N1) *thet2  +  (N3-N2) *thet3 
♦  (N4-N3) *thet2  +  (N5-N4) *thetl  +  (I-N5) flthet 
phi (I)  =  phia 
DO  135  J=l,  N15 

Rl  =  R15* (  sl5**J  -  1  )/(sl5-l) 

L  =  500*1  +  8*J 
x(L)  =  Rl*cos(phia) 
y(L)  =  Rl*sin(phia) 

WRITE (15, 1001)  L,  x (L) ,  y(L) 


thetl/2 . 


) / (sl5-l) 
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135  CONTINUE 

86  CONTINUE 

DO  916  I=N5,  N6 

phla  =  (Nl-1) * thet I  +  (N2-N1) *thet2 

♦  (N3-N2) *thet3  +  (N4-N3) *thet2 
+  (NS-N4) *thetl  +  (I-N5) *thet  >  thet/2 
DO  927  J=l,  N30 

R1  =  R30* (  s30**J  -  1  )/(s30-l) 

L  =  500*1  *  200  +  16 *J 
x(L)  a  Rl*cos  (phia) 
y(L)  =  Rl*sin(phla) 

WRITE  (15,  1001)  L,  x  (L)  ,  y  (L) 

927  CONTINUE 
916  CONTINUE 
c  Nodes  for  crack  flank 
I=N7 

phia  =  pi 
phi (I)  =  phia 
DO  11  J=l.  N1S 

R1  «  R15* (  slS**J  -  1  )/(sl5-l) 

l  =  soo*i  +  loo  +  a*j 

x(L)  =  Rl*cos  (phia) 
y(L)  =  Rl*sln(phia) 

WRITE (IS. 1001)  L.  x(L) .  y(L) 

11  CONTINUE 

c  Nodes  for  crack  tip 
DO  87  1*1,  N6 
L  =  500*1 

X(L)  =  0. 

y(L)  =  0. 

WRITE  (15.1001)  L.  x(L)  ,  y(L) 

L  =  500*1  +  200 
X(L)  =  0. 
y(L)  a  0. 

WRITE  (15,1001)  L.  x(L)  ,  y(L) 

87  CONTINUE 

L  =  500 *N7  +  100 
X(L)  =  0. 

y(L)  a  0. 

WRITE  (15,1001)  L,  x(L)  ,  y(L) 
c 

C  Elements  -  M  Is  the  element  no. 

1002  FORMAT (15,  8IS) 
c 

N8  a  12 
N9  *  16 
N99  =  17 
N10  a  18 
XC  a  3.81 

phia  =  (N7-N8) *thet 
Rcy  =  3.31/cos (phia) 
yc  a  Rcy*sin  (phia) 

c  Notice:  Fpr  the  shoulder,  the  following  constraint 
c  should  be  fullfllled: 
c  (xcl  ♦  rad  -  xu) **2  +  dy**2  =  rad* *2 
c  with  dyaxcl*tan (thet) 

|  rad  a  1.49 

alpha  *  1.  ♦  tan (thet) • tan (thet) 
beta  *  xc  -  rad 

gamma  a  (xc-rad) * (xc-rad)  -  rad*rad 


delta  =  fceta*beta  —  alpha*gamma 
xcl  =  (  beta  +  SQRT (delta)  )  /  alpha 
dy  =  xcl*tan (thet) 

Nadd  =  2*  (Nr+1 ) 

Nadi  =  2*Nr+l 

c  In  each  sector,  first  determine  boundary  nodes 
c  and  then  Input  elements 
DO  151  I=N1,  N2-1 
Rc(I)  »  xcl/cos  (phi  (I) ) 
arg  =  (  Rc (I) * (s75-l) /R75  )  +  1 
Nri  =  NINT (  LOG (arg) /LOG (s75)  ) 

LI  =  500*1  +  Nadd 
x(Ll)  =  xcl 

y  (LI)  =  Rc(I)*sin(phi(I)) 

WRITE  (15, 1001)  LI.  x  (LI)  ,  y(Ll) 

L2  =  500* (1+1)  +  Nadd 
Rc(I+l)  =  xcl/cos  (phi  (I >1) ) 
x (L2)  =  xcl 

y(L2)  =  Rc(I+l)  *sin(phi(I+l)) 

WRITE  (15, 1001)  L2,  x(L2)  ,  y(L2) 

L3  =  LI  +  200 

x (L3)  =  (  x (LI)  +  x (L2)  )/2. 
y(L3)  =  (  y (LI)  +  y(L2)  )/2. 

WRITE (15, 1001)  L3.  x(L3),  y(L3) 

DO  141  J=l,  Nri-1 
M  =  100*J  +  I 
Nodi  (M)  =  500*1  -*•  4*(J-1) 

Nod2 (M)  =  500*1  +  4*J 
Nod3 (M)  =  500* (1+1)  +  4*  J 
Nod4 (M)  =  500*  (1+1)  +  4*(J-1) 

Nod5 (M)  =  500*1  +  2* (2*J-1) 

Nod6  (M)  =  500*1  -*■  200  +  4*J 
Nod7 (M)  =  500*  (1+1)  +  2*(2*J-1) 

Nod8 (M)  =  500*1  +  200  +  4* (J-l) 

WRITE (16. 1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3 (M) , 

Nod4 (M) .  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 
L41  CONTINUE 
J  =  Nri 

LI  =  500*1  +  4* (J-l) 

L2  =  500*1  +  Nadd 
L3  =  500*  (1+1)  +  Nadd 
L4  =  500* (1+1)  +  4* (J-l) 

L12  =  500*1  +  Nadi 
x (L12)  =  (  x (LI)  +  x (L2)  )  /  2. 

y  (L12)  =  (  y  (LI)  +  y(L2)  )  /  2. 

WRITE (15, 1001)  L12,  x(L12) ,  y(L12) 

L34  =  500* (1+1)  ♦  Nadi 
x(L34)  =  (  x (L3)  +  x (L4)  )  /  2. 

y(L34)  =  (  y(L3)  +  y(L4)  )  /  2. 

WRITE (15, 1001)  L34,  x(L34) ,  y(L34) 

M  =  100*J  +  I 
Nodi (M)  =  LI 
Nod2(M)  =  L2 
Nod3 (M)  =  L3 
Nod4 (M)  =  L4 
Nod5(M)  =  LI  2 

Nod6  (M)  =  500*1  200  Nadd 

Nod7 (M)  =  L34 

Nod8 (M)  =  500*1  +  200  +  4* (J-l) 

WRITE  (16, 1002)  M,  Nodi (M) ,  Nod2 (M) ,  Nod3 (M) . 


3/3 


AD-A162  118  MIXED  MODE  I  AND  II  FULLV  PLASTIC  CRACK  GROWTH  FROM 
SIMULATED  HELD  DEFECTS(U)  MASSACHUSETTS  INST  OF  TECH 
CAMBRIDGE  DEPT  OF  MECHANICAL  ENGIN.  .  G  A  KARDOHATEAS 
UNCLASSIFIED  23  OCT  85  N88B14-82-K-8B25  F/G  13/8 


.a*  k_ 


193 


I 


t 

c 


Nod4  (M) ,  Nod5  (M)  .  Nod6  (M)  ,  Nod7  (M) ,  Nod8  (M) 

151  CONTINUE 

DO  152  I=N2,  N3-1 
Rc (I)  =  xcl/cos (phi (I) ) 
arg  =  (  Rc (I) *  (s-1) /RO  )  +  1 
Nrl  =  NINT (  LOG (arg) /LOG (s)  ) 

LI  *  500*1  ♦  Nadd 
x(Ll)  =  xcl 

y  (LI)  =  Rc(I)*sin(phi(I)) 

WRITE (15.1001)  LI,  x(Ll)  ,  y(Ll) 

L2  =  500* (1*1)  +  Nadd 
Rc(I+l)  *  xcl/cos  (phi  (I +1) ) 
x(L2)  =  xcl 

y(L2)  =  Rc(I+l)  *sln(phi(I+l)) 

WRITE  (15. 1001)  L2,  x(L2)  ,  y(L2) 

L3  =  LI  +  200 

x(L3)  =  (  x(Ll)  +  x(L2)  )/2. 
y(L3)  =  (  y (LI)  ♦  y(L2)  )/2. 

WRITE  (15. 1001)  L3.  x(L3),  y(L3) 

DO  142  J=l,  Nr 1-1 

M  =  100*J  ♦  I 

Nodi (M)  =500*1  +  2*(J-1) 

Nod2(M)  =  500*1  +  2*J 
Nod3(M)  =  500*  (1+1)  ♦  2*J 
Nod4  (M)  =  500* (1+1)  +  2*(J-1) 

Nod5(M)  =  500*1  +  (2*J-1) 

Nod6(M)  =  500*1  +  200  ♦  2*J 
Nod7 (M)  =  500* (1*1)  +  (2*J-1) 

Nod8 (M)  =500*1  ♦  200  +  2* (J-l) 

WRITE  (16. 1002)  M,  Nodl(M),  Nod2  (M)  ,  Nod3(M), 

-  •  Nod4  (M)  .  Nod5  (M)  ,  Nod6  (M)  ,  Nod7  (M)  ,  Nod8  (M) 

142  CONTINUE 
J  =  Nrl 

LI  =  500*1  +  2*  (J-l) 

L2  =  500*1  ♦  Nadd 
L3  =  500* (1+1)  +  Nadd 
L4  =  500* (1+1)  +  2* (J-l) 

LI 2  =500*1  +  Nadi 

x(L12)  =  (  x(Ll)  +  x(L2)  )  /  2. 

y (L12)  =  (  y (LI)  ♦  y(L2)  )  /  2. 

WRITE (15, 1001)  L12,  x(L12)  ,  y(L12) 

L34  =  500*  (1*1)  +  Nadi 
x (L34)  =  (  x(L3)  +  x (L4)  )  /  2. 

y (L34)  =  (  y(L3)  *  y(L4)  )  /  2. 

WRITE  (15, 1001)  L34,  x(L34)  ,  y(L34) 

M  =  100*J  +  I 
Nodi (M)  =  LI 
Nod 2 (M)  =  L2 
Nod3 (M)  =  L3 
Nod4(M)  =  L4 
Nod5 (M)  =  LI 2 

Nod6(M)  =  500*1  +  200  ♦  Nadd 
Nod7 (M)  =  L34 

Nod8 (M)  =  500*1  +  200  +  2* (J-l) 

WRITE (16, 1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3 (M) , 

Nod4 (M) .  NodS (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8  (M) 

152  CONTINUE 

DO  153  I=N3,  N4-1 
Rc(I)  =  xcl/cos  (phi  (I) ) 
arg  =  (  Rc (I) * (s7S-l) /R75  )  +  1 
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Nr I  »  NINT(  LOG (arg) /LOG {s 75)  ) 

LI  *  500*1  ♦  Nadd 
x(Ll)  =  xel 

y  (LI)  *  Rc(I)*sin(phi(I)) 

WRITE (15,1001)  LI,  x(Ll),  y(Ll) 

L2  =  500* (1+1)  +  Nadd 
Rc(I+l)  =  xcl/cos (phi (I+l) ) 
x(L2)  *  xcl 

y(L2)  =  Rc(H-l)*sin(phi(I*l)) 

WRITE  (15, 1001)  L2,  x(L2),  y(L2) 

L3  *  LI  ♦  200 

x(L3)  =  (  x(Ll)  +  x(L2)  )/2. 
y(L3)  =  (  y  (LI)  «•  y(L2)  )/2. 

WRITE (15, 1001)  L3,  x(L3),  y(L3) 

DO  143  J=l,  Nr 1-1 
M  *  100*J  +  I 
Nodi (M)  =  500*1  +  4* (J-l) 

Nod2 (M)  =  500*1  ♦  4*J 
Nod3(M)  *  500* (I+l)  +  4*J 
Nod4 (M)  =  500* (1+1)  +  4* (J-l) 

Nod5(M)  =  500*1  +  2* (2*J-1) 

Nod6 (M)  =  500*1  +  200  +  4*J 
Nod7(M)  =  500* (1+1)  +  2* (2* J-l) 

Nod8 (M)  a  500*1  ♦  200  ♦  4* (J-l) 

WRITE (16. 1002)  M,  Nodl(M)  ,  Nod2  (M) ,  Nod3(M), 

Nod4(M)  ,  Nod5(M),  Nod6  (M)  ,  Nod7  (M)  ,  Nod8(M) 
CONTINUE 
J  a  Nrl 

LI  a  500*1  4*  (J-l) 

L2  =  500*1  +  Nadd 
L3  a  500* (I+l)  ♦  Nadd 
L4  =  500* (I+l)  ♦  4* (J-l) 

LI 2  =  500*1  +  Nadi 
x (L12)  a  (  x(Ll)  ♦  x(L 2)  )  /  2. 

y(Ll2)  *  (  y (LI)  ♦  y(L2)  )  /  2. 

WRITE (15. 1001)  L12,  x(L12) ,  y(L12) 

L34  =  500* (I+l)  ♦  Nadi 
x (L34)  a  (  x(L3)  ♦  x(L4)  )  /  2. 

y (L34)  =  (  y(L3)  *  y(L4)  )  /  2. 

WRITE  (15. 1001)  L34,  x(L34),  y(L34) 

M  a  100*J  ♦  I 
Nodi  (M)  »  LI 
Nod2(M)  a  L2 
Nod3 (M)  a  L3 
Nod4(M)  *  L4 
Nod5(M)  a  L12 

Nod6 (M)  *  500*1  ♦  200  ♦  Nadd 
Nod7  (M)  a  L34 

Nod8(M)  a  500*1  ♦  200  +  4*  (J-l) 

WRITE (16, 1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3(M), 

Nod4 (M) ,  Nod5 (M) .  Nod6 (M) ,  Nod7 (M) ,  NodS (M) 

CONTINUE 

DO  154  I»M5,  N8-1 
Rc(I)  *  yc/sin(phi(I)) 
arg  a  (  Rc(I)  *  (s30-l)/R30  )  ♦  1 
Nrl  »  NINT(  LOG (arg) /LOG (s3C)  ) 

LI  *  500*1  ♦  Nadd 
X  (LI)  *  Rc  (I)  *coa  (phi  (I)  ) 

Y  (LI)  *  yc 

WRITE (15. 1001)  LI.  x(Ll).  y(Ll) 
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L2  =  500* (I+l)  +  Nadd 
Rc(I+l)  =  yc/sin(phi(I+l)) 
x(L2)  =  Rc(I+l)  *cos(phi(I+l)) 
y(L2)  =  yc 

WRITE (15, 1001)  L2,  x(L2).  y(L2) 

L3  =  LI  +  200 

x(L3)  *  (  x(Ll)  +  x (L2)  )/2. 
y(L3)  =  (  y (LI)  +  y(L2)  )/2. 

WRITE (15, 1001)  L3,  x(L3),  y(L3) 

DO  144  J=l,  Nr 1-1 
M  *  100*J  +  I 
Nodi  (M)  =  500*1  +  16*  (J — 1) 

Nod2 (M)  a  500*1  +  16 *J 
Nod3 (M)  =  500* (I+l)  +  16*J 
Nod4(M)  =  500*  (I+l)  +  16*  (J-l) 

NodS (M)  =  500*1  +  8* (2* J-l) 

Nod6(M)  *  500*1  +  200  +  16 *J 
Nod7(M)  =  500* (I+l)  +  8* (2*J-1) 

Nod8(M)  m  500*1  +  200  +  16* (J-l) 

WRITE (16. 1002)  M.  Nodl(M).  Nod2 (M) .  Nod3(M), 

Nod4(M),  NodS(M),  Nod6(M),  Nod7(M),  Nod8(M) 
144  CONTINUE 
J  =  Nrl 

LI  a  500*1  +  16* (J-l) 

L2  =  500*1  +  Nadd 
L3  =  500* (I+l)  +  Nadd 
L4  *  500*  (I+l)  +  16*  (J-l) 

L12  *  500*1  +  Nadi 

X(L12)  =  (  x(Ll)  +  x(L2)  )  /  2. 

y (L12)  a  (  y (LI)  +  y (L2)  )  /  2. 

WRITE (15. 1001)  L12,  x(L12),  y(L12) 

L34  a  500* (I+l)  +  Nadi 
x (L34)  =  (  x(L3)  +  x (L4)  )  /  2. 

y(L34)  =  (  y(L3)  +  y(L4)  )  /  2. 

WRITE (15. 1001)  L34.  x(L34),  y(L34) 

M  a  100*J  +  I 
Nodi (M)  =  LI 
Nod2 (M)  a  L2 
Nod3 (M)  a  L3 
Nod4(M)  a  L4 
Nod5(M)  a  LI  2 

Nod6 (M)  a  500*1  +  200  +  Nadd 
Nod7 (M)  a  L34 

NodS (M)  *  500*1  +  200  +  16* (J-l) 

WRITE (16, 1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3(M), 

Nod4  (M)  ,  NodS  (M)  ,  Nod6(M),  Nod7(M),  Nod8(M) 

154  CONTINUE 

DO  155  IaN8,  N7-1 

Rc(I+l)  =  -xc/cos(phi  (I+l)  ) 
arg  «  (  Rc  (I+l)  •  (s30-l)  /R30  )  +  1 
Nrl  *  NINT (  LOG (arg) /LOG (s30)  ) 

LI  a  500*  (I+l)  +  Nadd 
X(L1)  »  -xc 

Y  (LI)  a  Rc  (I+l)  *sln  (phi  (I+l) ) 

WRITE (15,1001)  LI.  x(Ll),  y(Ll) 

L2  a  500*1  +  Nadd 
Rc(I)  a  -xc/cos  (phi  (I) ) 
x(L2)  a  -xc 

y(L2)  »  Rc  (I) -sin  (phi  (I)) 

WRITE (15. 1001)  L2.  x(L2)  .  y(L2) 
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L3  *  L2  +  200 

x(L3)  a  (  x(H)  x (L2)  )/2. 

y(L3)  »  (  y (1*1)  ♦  y(L2)  )/2. 

WRITE  (15,1001)  L3,  x(L3).  y(L3) 

DO  145  J»l,  Nri-1 
M  *  100*J  +  I 
Nodi (M)  =  S00*I  +  16*(J-1) 

Nod2 (M)  =  500*1  +  16 *J 
Nod3(M)  »  500*  (I+l)  +  16 *J 
Nod4(M)  a  500*  (1+1)  ♦  16*(J-1) 

Nod5(M)  a  500*1  ♦  8* (2*J-1) 

Nod6(M)  a  500*1  ♦  200  +  16 *J 
Nod7 (M)  =  500* (1+1)  +  8* (2*J-1) 

Nbd8 (M)  =  500*1  ♦  200  ♦  16*(J-1) 

WRITE (16,1002)  M,  Nodi (M) ,  Nod2 (M) ,  Nod3(M), 

Nod4 (M) ,  Nod5 (M) .  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 
CONTINUE 
J  =  Nrl 

LI  =  500*1  +  16* ( J— 1) 

L2  =  500*1  ♦  Nadd 
L3  =  500*  (I+l)  ♦  Nadd 
L4  =  500*  (1+1)  +  16*(J-1) 

LI 2  a  500*1  +  Nadi 

x(L12)  =  (  x(Ll)  ♦  x(L2)  )  /  2. 

y (L12)  a  (  y (Ll)  ♦  y(L2)  )  /  2. 

WRITE (15. 1001)  L12,  x(L12).  y(L12) 

L34  =  500* (I+l)  +  Nadi 
x(L34)  a  (  x (L3)  ♦  x(L4)  )  /  2. 

y(L34)  a  (  y(L3)  +  y(L4)  )  /  2. 

WRITE (15, 1001)  L34,  x(L34) .  y(L34) 


M  a  100 *J 
Nodi  (M)  a 
Nod2  (M)  a 
Nbd3(M)  a 
Nod4 (M)  a 
Nod5 (M)  a 
Nod6  (M)  a 
Nod7 (M)  a 
Nod8 (M)  a 


♦  I 

Ll 

L2 

L3 

L4 

L12 

500*1 

L34 

500*1 


Nadd 
16*  (J-l) 


WRITE (16. 1002)  M,  Nodi (M) .  Nod2 (M) ,  Nod3(M), 

a* _  _a  ^  /K  MX  ft*  _  . _ _  _  _ _  '  ' 


Nod4 (M) ,  Nod5 (M) ,  Nod6 (M) . 

CONTINUE 

DO  156  I=N7+1,  N9-1 
Rc(I)  a  -xc/cos  (phi  (I) ) 
arg  a  (  Rc  (I)  *  (s30-l)  /R30  )  +  1 
Nrl  a  NINT(  LOG (arg) /LOG (s30) 
Ll  »  500*1  ♦  Nadd 
X(L1)  a  -xc 

y  (Ll)  a  Rc  (I)  *sln  (phi.(I) ) 

WRITE (15, 1001)  Ll,  x(Ll).  y(Ll) 
L2  a  500* (I+l)  ♦  Nadd 
Rc(I+l)  a  -xc/cos(phi(I+l)) 
x(L2)  a  -xc 

y(L2)  *  Rc  (I+l)  *sln  (phi  (I+l) ) 
WRITE  (IS,  1001)  L2,  x(L2)  ,  y(L2) 
L3  a  Ll  ♦  200 

X(L3)  a  (  x(Ll)  +  x(L2)  )/2. 
y(L3)  a  (  y  (Ll)  +  y(L2)  )/2. 
WRITE (15, 1001)  L3.  x(L3) ,  y(L3) 
DO  146  Jai,  Nr 1-1 


Nod7 (M) ,  Nod8 (M) 
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M  =  100*J  +  I 

Nodi (M)  *  500*1  +  16* (J-l) 

Nod2(M)  =  500*1  +  16 *J 
Nod3(M)  a  500*  (1+1)  ♦  16 *J 
Nod4(M)  =  500*  (1+1)  ♦  16*  (J-l) 

Nod5  (M)  a  500*1  ♦  8* (2*J-1) 

Nod6 (M)  a  500*1  +  200  +  16 *J 
Nod7(M)  a  500*  (1+1)  ♦  8* (2*J-1) 

Nod8 (M)  =  500*1  +  200  +  16* (J-l) 

WRITE  (16.1002)  M.  Nodi  (M)  .  Nod2  (M)  ,  Nod3(M). 

Nod4  (M)  ,  Nod5  (M)  .  Nod6  (M) ,  Nod7  (M)  .  Nod8  (M) 
CONTINUE 
J  =  Nrl 

LI  *  500*1  +  16* (J-l) 

L2  a  500*1  +  Nadd 
L3  a  500* (1+1)  +  Nadd 
L4  a  500* (1+1)  +  16* (J-l) 

L12  =  500*1  +  Nadi 

x(L12)  =  (  x(Ll)  ♦  x (L2)  )  /  2. 

Y (1*12)  =  (  y(Ll)  ♦  y(L2)  )  /  2. 

WRITE  (IS,  1001)  L12.  x(L12),  y(L12) 

L34  =  500* (1+1)  +  Nadi 
x(L34)  =  (  x(L3)  -*■  x (L4)  )  /  2. 

y (L34)  =  (  y(L3)  +  y(L4)  )  /  2. 

WRITE  (15. 1001)  L34,  x(L34)  ,  y(L34) 

M  a  ioo*j  +  I 
Nodi (M)  a  LI 
Nod2 (M)  a  L2 
Nod 3 (M)  a  L3 
Nod4(M)  a  L4 
Nod5 (M)  a  L12 

Nod6 (M)  a  500*1  +  200  ♦  Nadd 
Nod7  (M)  a  L34 

Nod8 (M)  a  500*1  +  200  +  16* (J-l) 

WRITE (16, 1002)  M,  Nodi (M) ,  Nod2 (M) ,  Nod3(M), 

Nod4 (M) .  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 
CONTINUE 

DO  157  I=N9,  N99-1 
Rc(I+l)  a  — yc/sln  (phi  (1+1) ) 
arg  =  (  Rc  (1+1)  *  (s30-l)  /R30  )  +  1 
Nrl  =  NINT(  LOG (arg) /LOG (s30)  ) 

LI  =  500* (I+l)  +  Nadd 
x(Ll)  =  Rc(I+l)*cos(phi(I+l)) 
y  (1*1)  =  -ye 

WRITE (15, 1001)  LI,  x (LI)  ,  y(Ll) 

L2  a  500*1  ♦  Nadd 
Rc(I)  =  -yc/sln  (phi  (I)) 
x(L2)  =  Rc(I)  *cos(phi(I)) 
y  (1*2)  =  -yc 

WRITE  (15, 1001)  L2,  x(L2),  y(L2) 

L3  *  L2  ♦  200 

x(L3)  a  (  x(Ll)  ♦  x(L2)  )/2. 
y (1*3)  =  (  y(Ll)  +  y(L2)  )/2. 

WRITE  (IS,  1001)  L3,  x(L3).  y(L3) 

DO  147  Jai.  Nr 1-1 
M  »  100*J  +  I 
Nodi (M)  a  500*1  ♦  16* (J-l) 

Nod2 (M)  a  500*1  +  16 *J 
Nod3(M)  a  500* (I+l)  ♦  16 *J 
Nod4(M)  a  500*  (1+1)  16*  (J-l) 
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Js-  NodS(M)  =  500*1  +  8*  (2*J-1) 

m  Node (M)  =  500*1  +  200  +  16*J 

H  Nod7(M)  =  500*  (1+1)  ♦  8*(2*J-1) 

iv  Nod8(M)  =  500*1  «■  200  +  16*  (J-l) 

V  WRITE (16, 1002)  M,  Nodi (M) ,  Nod2 (M) ,  Nod3(M), 

V:  -  Nod4(M)  ,  Nod5(M),  Nod6  (M)  ,  Nod7(M).  Nod8(M) 

147  CONTINUE 
&  J  s  Nri 

P  LI  *  500*1  ♦  16* (J-l) 

L2  *  500*1  ♦  Nadd 
L3  =  500* (I+l)  ♦  Nadd 
L4  *  500* (I+l)  +  16* (J-l) 

L12  *  500*1  ♦  Nadi 


x(L12)  *  (  x(Ll)  +  x(L2)  )  /  2. 

y(L12)  =  (  y (LI)  -  y(L2)  )  /  2. 

WRITE (IS. 1001)  L12,  x(L12) ,  y(L12) 

L34  =  500*  (I+l)  +  Nadi 
x (L34)  =  (  x (L3)  x (L4)  )  /  2. 

y (L34)  =  (  y(L3)  ♦  y(L4)  )  /  2. 

WRITE (15, 1001)  L34,  x(L34) ,  y(L34) 

M  =  100*J  +  I 
Nodi (M)  =  LI 
Nod2 (M)  =  L2 
Nod 3 (M)  =  L3 
Nod4 (M)  a  L4 
Nod5(M)  a  L12 

Nod6 (M)  a  500*1  ♦  200  ♦  Nadd 
Nod7 (M)  a  L34 

Nod8 (M)  a  500*1  ♦  200  ♦  16* (J-l) 

WRITE (16.1002)  M,  Nodi (M) ,  Nod2 (M) ,  Nod3 (M) , 

Nod4(M),  Nod5(M).  Nod6(M),  Nod7  (M)  ,  Nod8(M) 

157  CONTINUE 

DO  158  I=N99,  N10-1 
Rc(I)  a  -yc/sln  (phi  (I)  ) 
arg  =  (  Rc  (I)  *  (s30-l)  /R30  )  ♦  i 
Nrl  =  NINT (  LOG (arg) /LOG (s30)  ) 

LI  =  500*1  ♦  Nadd 
x (LI)  =  Rc(I)  *cos(phl(I)) 

/(LI)  =  —ye 

WRITE (15, 1001)  LI.  x (LI) ,  y(Ll) 

L2  *  500* (1+1)  ♦  Nadd 
Rc(I+l)  a  -yc/sin(phi(I+l)) 
x(L2)  a  Rc(I+l)*cos(phi(I+l)) 

y(L2)  =  -yc 

WRITE (15, 1001)  L2,  x(L2),  y(L2) 

L3  »  LI  +  200 

x(L3)  a  (  x(Ll)  +  x(L2)  )/2. 

y(L3)  a  (  y (LI)  +  y(L2)  )/2. 

WRITE (15, 1001)  L3.  x(L3),  y(L3) 

DO  148  J*i,  Nr 1-1 
M  a  100*J  ♦  I 

Nodi (M)  =  500*1  ♦  16* (J-l) 

Nod2(M)  »  500*1  ♦  16*J 

Nod3(M)  a  500* (I+l)  +  16*J 

Nod4(M)  a  500* (I+1)  ♦  16* (J-l) 

Nod5 (M)  a  500*1  ♦  8* (2*J-1) 

Nod6 (M)  a  500*1  ♦  200  ♦  16*J 
Nod7(M)  a  500*  (I*l)  ♦  8*(2*J-1) 

Nod8 (M)  a  500*1  ♦  200  >  16* (J-l) 

WRITE (16. 1002)  M.  Nodl(M).  Nod2 (M) ,  Nod3(M), 
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Nod4(M)  ,  NodS(M),  Nod6  (M) ,  Nod7(M),  Nod8(M) 
148  CONTINUE 
J  =  Nrl 

LI  =  500*1  ♦  16*(J-1) 

L2  =  500*1  +  Nadd 
L3  =  500*  (1+1)  >  Nadd 
L4  =  500*  (1+1)  +  16*(J-1) 

L12  =  500*1  Nadi 
x (L12)  =  (  x(Ll)  +  x(L2)  )  /  2. 

y (L12)  =  (  y (LI)  ♦  y(L2)  )  /  2. 

WRITE (15, 1001)  L12.  x(L12),  y(L12) 

L34  =  500* (1*1)  +  Nadi 
x(L34)  =  (  x(L3)  +  x(L4)  )  /  2. 

y(L34)  =  (  y(L3)  +  y(L4)  )  /  2. 

WRITE (15, 1001)  L34,  x(L34),  y(L34) 

M  =  100*J  +  I 
Nodi (M)  =  LI 
Nod2 (M)  =  L2 
Nod3(M)  =  L3 
Nod4(M)  =  L4 
Nod5  (M)  =  L12 

Nod6 (M)  =500*1  +  200  +  Nadd 
Nod7 (M)  =  L34 

Nod8(M)  =  500*1  +  200  +  16* (J-l) 

WRITE (16,1002)  M.  Nodl(M).  Nod2 (M) ,  Nod3(M), 

Nod4(M).  Nod5  (M)  ,  Nod6  (M)  ,  Nod7(M).  Nod8(M) 
158  CONTINUE 

DO  159  I=N10,  N6-1 
Rc(I+l)  =  xc/cos (phi (1+1) ) 
arg  =  (  Rc(I+l)  *  (s30-l)/R30  )  +  1 
Nrl  =  NINT (  LOG (arg) /LOG (s30)  ) 

LI  =  500*  (1+1)  +  Nadd 
x  (LI)  =  xc 

y  (LI)  =  Hc(X*>l)  *sin  (phi  (1+1) ) 

WRITE (15, 1001)  LI,  x(Ll),  y(Ll) 

L2  =  500*1  +  Nadd 
Rc  (I)  =  xc/cos  (phi  (I) ) 
x(L2)  =  xc 

y(L2)  =  Rc(I)  *sin(phi(I)) 

WRITE (15, 1001)  L2,  x(L2),  y(L2) 

L3  =  L2  ♦  200 

x(L3)  =  (  x(Ll)  +  x(L2)  )/2. 
y(L3)  =  (  y (LI)  +  y(L2)  )/2. 

WRITE (15. 1001)  L3,  x(L3),  y(L3) 

DO  149  J=l.  Nri-1 
M  =  100*J  +  I 
Nodi (M)  =  500*1  +  16* (J-l) 

Nod2(M)  *  500*1  +  16 *J 
Nod3(M)  =  500*  (I+l)  ♦  16* J 
Nod4 (M)  =  500* (I+l)  ♦  16* (J-l) 

Nod5(M)  =  500*1  ♦  8* (2* J-l) 

Nod6 (M)  *  500*1  ♦  200  ♦  16*J 
Nod7(M)  =  500*  (I+l)  +  8*  (2*J-1) 

Nod8(M)  =  S00*I  200  ♦  16*  (J-l) 

WRITE (16, 1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3(M) , 

Nod4 (M) .  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 
149  CONTINUE 
J  =  Nrl 

LI  *  500*1  16*  (J-l) 

L2  ■  500*1  «■  Nadd 
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L3  *  500*  (I •*■!)  ♦  Nadd 
L4  »  500*  (I+1J  ♦  16*(J-1) 

L12  =  500*1  ♦  Nadi 
x (L12)  a  (  x (LI)  +  x (L2)  )  /  2. 

y(L12)  =  (  y (LI)  ♦  y(L2)  )  /  2. 

WRITE (15, 1001)  L12,  x(L12),  y(L12) 

L34  a  500* (I+l)  +  Nadi 
x (L34)  a  (  X(L3)  +  x(L4)  )  /  2. 

y (L34)  =  (  y(L3)  ♦  y(L4)  )  /  2. 

WRITE  (15, 1001)  L34,  x(L34),  y(L34) 

M  =  100*J  ♦  I 
Nodi (M)  =  LI 
Nod2 (M)  a  L2 
Nod3 (M)  =  L3 
Nod4 (M)  =  L4 
Nod5 (M)  =  L12 

Nod6(M)  =  500*1  +  200  +  Nadd 
Nod7  (M)  =  L34 

Nod8 (M)  =  500*1  ♦  200  +  16*(J-1) 

WRITE (16, 1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3 (M) , 

Nod4 (M) ,  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 

159  CONTINUE 
I=N7 

Rc(I)  =  -xc/coa  (phi  (I) ) 
arg  a  (  Rc(I)  *  (s30-l)/R30  )  ♦  1 
Nrl  =  NINT(  LOG  (arg) /LOG  (s30)  ) 

LI  a  500*1  +  100  +  Nadd 
x(Ll)  a  — xc 

y  (LI)  a  Rc(I)*sin(phi(I)) 

WRITE (15, 1001)  LI,  x(Ll),  y(Ll) 

L2  a  500* (1+1)  ♦  Nadd 
Rc(I+l)  a  -xc/cos  (phi  (I+l) ) 

X (L2)  =  -xc 

y(L2)  =  Rc(I-*-l)*sln(phi(I+l)) 

WRITE  (IS,  1001)  L2,  x(L2)  ,  y(L2) 

L3  =  500*1  +  200  +  Nadd 
x(L3)  =  (  x(Ll)  +  x(L2)  )/2. 
y(L3)  =  (  y (LI)  *  y(L2)  )/2. 

WRITE  (15, 1001)  L3.  x(L3),  y(L3) 

DO  161  J=l,  Nr 1-1 
M  a  100*J  +  I 

Nodi  (M)  =  500*1  +  100  +  16*(J-1) 

Nod2 (M)  =  500*1  +  100  +  16*J 
Nod3(M)  a  500*  (1+1)  +  16* J 
Nod4 (M)  =  500*  (I+l)  +  16*(J-1) 

Nod5 (M)  =  500*1  +  100  +  8* (2*J-1) 

Nod6(M)  =  500*1  +  200  +  16 *J 
Nod7 (M)  =  500*  (1+1)  +  8* (2*J-1) 

Nod8 (M)  =  500*1  +  200  +  16* (J-l) 

WRITE (16, 1002)  M,  Nodi (M) ,  Nod2 (M) ,  Nod3(M), 

Nod4 (M) ,  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 
161  CONTINUE 
J  *  Nrl 

LI  a  500*1  ♦  100  ♦  16* (J-l) 

L2  =  500*1  ♦  100  ♦  Nadd 

L3  =  500*  (1*1)  ♦  Nadd 
L4  a  500*  (I+l)  ♦  16*  (J-l) 

L12  *  500*1  ♦  100  ♦  Nadi 
x (L12)  a  (  x(Ll)  ♦  x(L2)  )/  2. 

y (L12)  =  (  y (LI)  ♦  y(L2)  )  /  2. 


WRITE (IS, 1001)  L12,  x(L12),  y(L12) 

L34  =  500* (I+l)  ♦  Nadi 
x (L34)  =  (  x (L3)  +  x(L4)  )  /  2. 

Y (1*34)  =  (  y(L3)  ♦  y(L4)  )  /  2. 

WRITE  (15,1001)  L34,  x(L34),  y(L34) 

M  =  100*J  +  I 
Nodi (M)  =  LI 
Nod2(M)  =  L2 
Nod3 (M)  =  L3 
Nod4(M)  »  L4 
Nod5(M)  =  L12 

Nod6 (M)  =  500*1  +  200  +  Nadd 
Nod7 (M)  =  L34 

Nod8 (M)  =  500*1  +  200  +  16*(J-1) 

WRITE (16, 1002)  M,  Nodi (M) ,  Nod2 (M) ,  Nod3 (M) , 

Nod4 (M) ,  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (H) 

I=N6 

Rc(l)  =  xc/cos  (phi  (1) ) 
arg  =  (  Rc(l)  *  (s30-l)/R30  )  +  1 
Nrl  =  NINT(  LOG  (arg) /LOG  (s30)  ) 

LI  =  500*1  +  Nadd 
x  (Ll)  =  xc 

y  (LI)  =  Rc(l)  *sln(phi(l)) 

WRITE (15. 1001)  Ll,  x(Ll).  y(Ll) 

L2  =  500*1  +  Nadd 
Rc (I)  =  xc/cos (phi (I)) 
x(L2)  =  xc 

y(L2)  *  Rc(I)*sln(phi(I)) 

WRITE  (15. 1001)  L2,  x(L2),  y(L2) 

L3  »  L2  +  200 

x(L3)  =  (  x(Ll)  +  x(L2)  )/2. 

y(L3)  =  (  y (Ll)  ♦  y(L2)  )/2. 

WRITE (15. 1001)  L3,  x(L3)  ,  y(L3) 

DO  162  J=l.  Nr 1-1 
M  =  100*J  +  I 
Nodi (M)  =  500*1  +  16*(J-1) 

Nod2 (M)  *  500*1  ♦  16 *J 
Nod3 (M)  «  500*1  +  16 *J 
Nod4 (M)  *  500*1  +  16*(J-1) 

NodS(M)  =  500*1  ♦  8* (2*J-1) 

Nod6(M)  =  500*1  +  200  +  16*J 
Nod7(M)  =  500*1  +  8* (2*J-1) 

Nod8 (M)  =  500*1  +  200  ♦  16*(J-1) 

WRITE (16, 1002)  M,  Nodi (M) ,  Nod2 (M)  ,  Nod3(M). 

Nod4(M)  ,  Nod5  (M)  ,  Nod6  (M)  ,  Nod7(M),  Nod8(M) 
CONTINUE 
J  =  Nrl 

Ll  =  500*1  ♦  16*(J-1) 

L2  =  500*1  ♦  Nadd 
L3  *  500*1  «•  Nadd 
L4  =  500*1  16*(J-1) 

L12  *  500*1  ♦  Nadi 
x (L12)  *  (  x(Ll)  >  x(L2)  )  /  2. 

y (L12)  =  (  y (Ll)  >  y(L2)  )  /  2. 

WRITE (15, 1001)  L12,  x(L12).  y(Ll2) 

L34  *  500*1  *  Nadi 
x (L34)  =  (  x (L3)  ♦  x (L4)  )  /  2. 

y(L34)  -  (  y(L3)  ♦  y(L4)  )  /  2. 

WRITE (15, 1001)  L34,  x(L34)  ,  y(L34) 

M  =  100*J  +  I 


Nodi (M)  =  LI 
Nod2(M)  =  L2 
Nod3(M)  =  L3 
Nod4(M)  =  L4 
NodS(M)  =  L12 

Nod6 (M)  =  500*1  +  200  +  Nadd 
Nod7(M)  =  L34 

Nodd(M)  =  500*1  +  200  +  16*(J-1) 

WRITE (16, 1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3(M). 

Nod4 (M) ,  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 

I=N5-1 

Rc(I)  =  yc/sln(phi(I) ) 
arg  =  (  Rc (I) * (sl5-l) /R15  )  +  1 
Nrl  =  NINT (  LOG (arg) /LOG (S15)  ) 

LI  =  500*1  +  Nadd 
x(Ll)  *  Rc(I)  *cos(phl(I)) 
y  (LI)  =  yc 

WRITE  (15, 1001)  LI,  x(Ll)  ,  y(Ll) 

Rc(I+l)  =  yc/sin  (phi  (I+l) ) 

L2  =  500* (I+l)  +  Nadd 
x(L2)  =  Rc(I+l)  *cos(phl(I+l)) 
y(L2)  =  yc 

WRITE  (15. 1001)  L2,  x(L2)  .  y(L2) 

L3  =  LI  +  200 

x(L3)  =  (  x (LI)  +  x(L2)  )/2. 
y(L3)  =  (  y (LI)  -  y(L2)  )/2. 

WRITE (15. 1001)  L3,  x(L3)  ;  y(L3) 

DO  163  J=l,  Nr 1-1 
M  =  100*J  +  I 
Nodi (M)  =  500*1  +  8* (J-l) 

Nod2 (M)  *  500*1  +  8*J 
Nod3 (M)  =  500* (I+l)  +  8*J 
Nod4 (M)  =  500* (I+l)  +  8* (J-l) 

Nod5 (M)  *  500*1  +  4* (2*J-1) 

Nod6 (M)  =  500*1  +  200  +  8*J 
Nod7 (M)  =  500* (I+l)  +  4* (2*J-1) 

Nod8 (M)  =  500*1  +  200  +  8* (J-l) 

WRITE (16. 1002)  M,  Nodi (M) ,  Nod2 (M) ,  Nod3(M), 

Nod4 (M) ,  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 
CONTINUE 
J  =  Nrl 

LI  =  500*1  +  8* (J-l) 

L2  =  500*1  +  Nadd 
L3  =  500* (I+l)  +  Nadd 
L4  «  500*  (I+l)  +  8*  (J-l) 

L12  =  500*1  +  Nadi 
x (L12)  =  (  x (LI)  +  x (L2)  )  /  2. 

y (L12)  =  (  y (LI)  +  y(L2)  )  /  2. 

WRITE (15, 1001)  L12.  x(Ll2) ,  y(Ll2) 

L34  =  500* (I+l)  +  Nadi 
x(L34)  =  (  x (L3)  +  x (L4)  )  /  2. 

y(L34)  *  (  y(L3)  +  y(L4)  )  /  2. 

WRITE (15. 1001)  L34,  x(L34) .  y(L34) 

M  »  100*J  +  I 
Nodi (M)  *  LI 
Nod2 (M)  *  L2 
Nod3 (M)  »  L3 
Nod4 (M)  =  L4 
Nod5 (M)  »  L12 

Nod6 (M)  =  500*1  +  200  ♦  Nadd 


Nod7(M)  =  L34 

Nod8 (M)  =  500*1  +  200  ♦  8* (J-l) 

WRITE (16, 1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3(M), 

Nod4 (M) ,  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 

I=N4 

Rc(I)  =  xcl/cos  (phi  (I) ) 
arg  =  (  Rc(I)  *  (sl5-l)/R15  )  +  1 
Nri  =  NINT (  LOG (arg) /LOG (sl5)  ) 

LI  =  500*1  +  Nadd 
x(Ll)  =  xcl 

y  (LI)  =  Rc(I)*sin(phi(I)) 

WRITE  (15,1001)  LI,  x(Ll)  ,  y(Ll) 

L2  =  500* (1+1)  +  8*Nri 
x(L2)  =  Rc(I)  *cos(phi(I+l)) 
y(L2)  =  Rc(I)  *sin03hi(I+l) ) 

WRITE  (15, 1001)  L2.  x(L2).  y(L2) 

L3  =  LI  +  200 

x(L3)  =  (  x  (LI)  +  x(L2)  )/2. 
y(L3)  =  (  y (LI)  +  y(L2)  )/2. 

WRITE  (15, 1001)  L3,  x(L3)  ,  y(L3) 

L4  =  500*1  +  2* (Nr +3) 
x  (L4)  =  xcl 
y(L4)  =  yc 

WRITE  (15, 1001)  L4,  x(L4)  ,  y(L4) 

L5  =  500* (1+1)  +  Nadd 
L6  =  L4  +  1 

x (L6)  =  (  x(L4)  ♦  x(Ll)  )/2. 

y(L6)  =  (  y(L4)  >  y(Ll)  )/2. 

WRITE (15. 1001)  L6.  x(L6).  y(L6) 

L7  =  L4  +  2 

x(L7)  =  (  x (L4)  +  x (L5)  )/2. 

y(L7)  =  (  y(L4)  +  y(L5)  )/2. 

WRITE  (15, 1001)  L7,  x(L7)  ,  y(L7) 

L8  =  L2  ♦  8 

x(L8)  =  (  x(L2)  +  x (L5)  )/2. 

y(L8)  =  (  y(L2)  ♦  y(L5)  )/2. 

WRITE (15, 1001)  L8.  x(L8) ,  y(L8) 

DO  166  J=l,  Nri-1 
M  =  100*J  ♦  I 

Nodi (M)  =  500*1  ♦  8* (J-l) 

Nod2(M)  =  500*1  ♦  8*J 

Nod3(M)  =  500* (I+l)  +  8*J 

Nod4 (M)  =  500* (I>1)  ♦  8* (J-l) 

Nod5 (M)  =  500*1  ♦  4* (2* J-l) 

Nod6 (M)  =  500*1  +  200  ♦  8*J 
Nod7(M)  =  500*  (1  +  1)  ♦  4*  (2*J-1) 

Nod8 (M)  =  500*1  ♦  200  +  8* (J-l) 

WRITE (16, 1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3(M). 

Nod4 (M) ,  Nod5 (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (H) 
CONTINUE 
J  =  Nri 
M  =  100*J  I 
Nodi (M)  =  500*1  +  8* (J-l) 

Nod2 (M)  =500*1  +  Nadd 
Nod 3 (M)  =  500* (I+l)  ♦  8* J 
Nod4  (M)  =  500*  (I>1)  ♦  8*  (J-l) 

Nod5 (M)  =  500*1  ♦  4*(2*J-1) 

Nod6 (M)  =500*1  ♦  200  +  Nadd 
Nod7 (M)  =  500* (I+l)  +  4* (2* J-l) 

Nod8 (M)  =500*1  *  200  +  8* (J-l) 


WRITE (16, 1002)  M,  Nodi (M) .  Nod2 (M) ,  Nod3(M), 

Nod4 (M) ,  NodS (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 
J  =  Nri  +  1 
M  =  100*J  +  I 
Nodi (M)  =  L2 
Nod2 (M)  =  LI 
Nod3 (M)  =  L4 
Nod4 (M)  =  L5 
NodS (M)  =  L3 
Nod6 (M)  =  L6 
Nod7(M)  =  L7 
Nod8 (M)  =  L8 

WRITE (16,1002)  M.  Nodi (M) ,  Nod2 (M) ,  Nod3 (M) , 

Nod4 (M) ,  NodS (M) ,  NodS (M) ,  Nod7 (M) ,  Nod8 (M) 
DO  164  1=1,  Nl-1 

beta  =  (xcl+rad) ‘cos (phi (I+l) )  +  dy'sin  (phi  (I+l) ) 
gamma  =  (xcl+rad) * (xcl+rad)  +  dy*dy  -  rad*rad 
Rc(I+l)  =  beta  -  SQRT (beta ‘beta-gamma ) 
arg  =  (  Rc (I+l) * (slS-1) /R15  )  +  1 
Nri  =  NINT (  LOG (arg) /LOG (slS)  ) 

LI  =  500‘(I+1)  +  Nadd 
x  (LI)  =  Rc(I+l)  ‘cos  (phi  (I+l)) 
y  (Ll)  =  Rc (I+l) ‘sin (phi (I+l)) 

WRITE (15, 1001)  Ll.  x (Ll) .  y(Ll) 

L2  =  500*1  +  Nadd 
phi 3  =  (  phi (I)  +  phi (1+1)  )  /  2. 
beta  =  (xcl+rad)  ‘cos  (phi3)  +  dy*sin(phi3) 
gamma  =  (xcl+rad) * (xcl+rad)  +  dy*dy  -  rad*rad 
Rc3  =  beta  -  SQRT (beta ‘bet a-gamma) 

L3  =  L2  +  200 
x(L3)  =  Rc3*cos  (ph!3) 
y(L3)  =  Rc3*sin  (ph!3) 

WRITE  (IS.  1001)  L3,  x(L3)  .  y(L3) 

DO  165  J=l.  Nr 1-1 
M  =  100*J  ♦  I 
Nodi (M)  =  500*1  +  8* (J-l) 

Nod2 (M)  =  500*1  +  8* J 
Nod3 (M)  =  500* (1+1)  +  8* J 
Nod4  (M)  =  500*  (1  +  1)  +  8*  (J-l) 

NodS (M)  =  500*1  +  4* (2*J-1) 

Nod6 (M)  =  500*1  +  200  +  8*J 
Nod7 (M)  =  500  * (1  +  1)  +  4* (2*J-1) 

Nod8 (M)  =500*1  +  200  +  8* (J-l) 

WRITE (16, 1002)  M.  Nodi (M) ,  Nod 2 (M) ,  Nod 3 (M) , 

Nod4 (M)  .  NodS (M) ,  Nod6 (M) ,  Nod7 (M) ,  Nod8 (M) 
CONTINUE 
J  *  Nri 

Ll  *  500*1  +  8* (J-l) 

L2  =  500*1  +  Nadd 
L3  =  500* (1+1)  ♦  Nadd 
L4  =  500* (I+l)  +  8* (J-l) 

L12  =  500*1  +  Nadi 
x (L12)  »  (  x (Ll)  +  x (L2)  )  /  2. 

y (L12)  =  (  y (Ll)  +  y(L2)  )  /  2. 

WRITE (15, 1001)  L12,  x(L12),  y(L12) 

L34  =  SCO* (I+l)  +  Nadi 
x (L34)  =  (  x (L3)  +  x (L4)  )  /  2. 

y(L34)  =  (  y(L3)  +  y(L4)  )  /  2. 

WRITE (15. 1001)  L34,  x(L34) ,  y (L34) 

M  =  100*J  ♦  I 
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Nodi (M)  =  LI 
Nod2(M)  =  L2 
Nod 3 (M)  =  L3 
Nod4(M)  =  L4 
Nod5 (M)  =  L12 
Nod6 (M)  =  S00*I 
Nod7(M)  =  L34 
Nod8 (M)  =500*1 
WRITE (16, 1002) 
Nod4  (M) , 

CONTINUE 


♦  200  ♦  Nadd 

+  200  +  8*  (J-l) 

M.  Nodi (M) ,  Nod2 (M) .  Nod3(M) , 

Nod5 (M) .  Nod6 (M)  .  Nod7 (M) ,  Nod8 (M) 


c 

c  This  portion  of  the  program  creates  the  MPC  constraints 
2001  FORMAT (515) 

L  =  2 
c 

c  15-30  interface 

c  Nri  referes  to  the  30  i.e.  4  elements 
I  =  11 
Nri  =  4 

DO  111  J=l,  Nr 1-1 

LI  =  500*1  +  4* (4*J-3) 

Lll  =  500*1  +  4* (4*J-1) 

L2  =  500*1  +  16*  (J-l) 

L3  =  500*1  16*J 

L4  =  500*1  +  8* (2*J-1) 

WRITE (17, 2001)  L,  LI,  L2,  L3,  L4 
WRITE (17, 2001)  L.  Lll,  L2,  L3,  L4 
111  CONTINUE 
1  =  1 
Nri  =  4 

DO  12  J=l.  Nr 1-1  . 

LI  =  500*1  +  4* (4*J-3) 

Lll  =  500*1  4*  (4* J-l) 

L2  =  500*1  *  16* (J-l) 

L3  =  500*1  *  16 *J 
L4  =  500*1  8*  (2*J-1) 

WRITE (17,2001)  L,  Ll,  L2,  L3,  L4 
WRITE (17, 2001)  L,  Lll,  L2,  L3,  L4 

12  CONTINUE 
c 

c  7.5-15  interface 

c  Nri  referes  to  the  15  i.e.  7  elements 
1  =  3 
Nri  =  7 

DO  13  J=l.  Nr 1-1 

Ll  =  500*1  +  2* (4*J-3) 

Lll  =  500*1  ♦  2* (4*J-1) 

L2  =  500*1  +  8* (J-l) 

L3  =  500*1  +  8*J 
L4  =  500*1  +  4*(2*J-1) 

WRITE (17,2001)  L.  Ll,  L2,  L3,  L4 
WRITE (17,2001)  L,  Lll.  L2,  L3,  L4 

13  CONTINUE 
1=9 
Nri  =  8 

DO  14  J=l,  Nr 1-1 

Ll  =  500*1  ♦  2* (4*J-3) 

Lll  =  500*1  ♦  2* (4* J-l) 

L2  =  500*1  ♦  8* (J-l) 
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L3  =  500*1  ♦  8*J 
L4  =  500*1  ♦  4* (2*J-1) 

WRITE (17,2001)  L,  LI,  L2.  L3,  L4 
WRITE  (17, 2001)  L,  71.  L2.  L3,  L4 
CONTINUE 


15 


c  3.75-7.5  interface 

c  Nri  referee  to  the  7.5  l.e.  13  elements 
I  =  4 
Nri  =  13 
DO  15  J=l.  Nri-1 

LI  =  500*1  +  (4*J-3) 

Lll  =500*1  +  (4*J-1) 

L2  =  500*1  +  4* ( J — 1) 

L3  =  500*1  ♦  4*J 
L4  =  500*1  2*  (2*J-1) 

WRITE (17, 2001)  L,  LI,  L2.  L3,  L4 
WRITE (17, 2001)  L,  Lll,  L2,  L3.  L4 
CONTINUE 
1=8 
Nri  =  15 
DO  16  J=l,  Nr 1-1 

LI  =  500*1  +  (4*J-3) 

Lll  =  500*1  +  (4*J— 1) 

L2  =  500*1  +  4* (J— 1) 

L3  =  500*1  -►  4*J 
L4  =  500*1  +  2* (2*J— 1) 

WRITE (17. 2001)  L.  LI.  L2.  L3,  L4 
WRITE (17.2001)  L.  Lll,  L2.  L3.  L4 
CONTINUE 
STOP 
END 
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